SPECTRAL PROPERTIES OF A LIMIT-PERIODIC 
SCHRODINGER OPERATOR IN DIMENSION TWO. 



YULIA KARPESHINA AND YOUNG-RAN LEE 



Abstract. We study Schrodinger operator H — —A + V(x) in dimension two, 
V(x) being a limit-periodic potential. We prove that the spectrum of H contains 
a semiaxis and there is a family of generalized eigenfunctions at every point of 
this semiaxis with the following properties. First, the eigenfunctions are close 
to plane waves e 1 ^'^ at the high energy region. Second, the isoenergetic curves 
in the space of momenta k corresponding to these eigenfunctions have a form of 
slightly distorted circles with holes (Cantor type structure). Third, the spectrum 
corresponding to the eigenfunctions (the semiaxis) is absolutely continuous. 

1. Introduction. 

We study the operator 

H = -A + V(x) (1.1) 
in two dimensions, V(x) being a limit-periodic potential: 

oo 

V{x) = Y J V r {x), (1.2) 

r=l 

where {V^}^2 =1 is a family of periodic potentials with doubling periods and decreas- 
ing Loo-norms. Namely, V r has orthogonal periods 2 r di, 2 r ~ 1 c?2 and ||Vr||oo < 
Cexp(—2 vr ) for some rj > rjo > 0. Without loss of generality we assume that (7 = 1, 
di = (di,0), c?2 = (0,d%) and fg V r (x)dx = 0, Q r being the elementary cell of peri- 
ods corresponding to V r (x). We assume that all V r (x) are trigonometric polynomials 
with the lengths growing at most linearly with period. Namely, there exists a positive 
number Rq < oo, such that each potential admits Fourier representation: 



V r {x)= ^2 v r , g expi{2 r+1 q,x), q = 2n(^-,^ 

(•, •) being the canonical dot product in R 2 . 

The one-dimensional analog of (jl.ip . (|1.2p is already thoroughly investigated. It is 
proven in PQ-[I] that the spectrum of the operator H\u = —u" + Vu is generically 
a Cantor type set. It has a positive Lebesgue measure [U [6]. The spectrum is 
absolutely continuous [U [2], [5]- [9]. Generalized eigenfunctions can be represented 



Date: August 30, 2010. 

Supported in part by NSF grant DMS-0800949 (Y.K.) and the National Research Foundation of 
Korea (NRF)-grant 2010-0015575 (Y.-R.L.). 

1 



2 



Y. KARPESHINA AND Y.-R. LEE 



in the form of e lkx u(x), u{x) being limit-periodic El [7]. The case of a complex- 
valued potential is studied in [10]. Integrated density of states is investigated in 
[TT]-[T4]. Spectral properties of Schrodinger operators in l 2 {1) with limit-periodic 
potentials are recently investigated in [15J. It regards such potentials as generated by 
continuous sampling along the orbits of a minimal translation of a Cantor group. It 
is shown that the spectrum is a Cantor set of positive Lebesgue measure and purely 
absolutely continuous for a dense set of sampling functions, and it is a Cantor set of 
zero Lebesgue measure and purely singular continuous for a dense G$ set of sampling 
functions. Properties of eigenfunctions of discrete multidimensional limit-periodic 
Schrodinger operator are studied in [16]. As to the continuum multidimensional case, 
it is proven in [14] that the integrated density of states for (jl.ip is the limit of densities 
of states for periodic operators. 

We concentrate here on properties of the spectrum and eigenfunctions of (jl.ip . (|1.2p 
in the high energy region. We prove the following results for the two-dimensional case. 

(1) The spectrum of the operator (jl.ip . (jl.2p contains a semiaxis. A proof of an 
analogous result by different means can be found in the paper [17j . In [T7], the 
authors consider the operator H = (— A) l +V, 81 > d+3, d ^ l(mod4), d being 
the dimension of the space. This obviously includes our case I = 1, d = 2. 
However, there is an additional rather strong restriction on the potential 
V(x) in [T7], which we don't have here: in [17] all the lattices of periods Q r 
of periodic potentials V r need to contain a nonzero vector 7 in common, i.e., 
V(x) is periodic in a direction 7. 

(2) There are generalized eigenfunctions ^^(k, x), corresponding to the semiaxis, 
which are close to plane waves: for every k in an extensive subset of M 2 , 
there is a solution VP^fe, a?) of the equation -ffvl'oo = ^00^00 which can be 
described by the formula: 

\MM) = e*^ (l + Uoo (k,x)) , (1.3) 



IKoo||= |A ; Koo O(|A:|- 71 ), 71 >0, (1.4) 

where Uoo(k,x) is a limit-periodic function: 

00 

Uoo(k, x) = y~] u r (k, x), (1.5) 

r=l 

u r (k,x) being periodic with periods 2 r ~ 1 di, 2 r ~ x di- The eigenvalue Aoo(fc) 
corresponding to ^oo(k,x) is close to \k\ 2 : 

Aoo(£)= |feV j£| 2 + 0(|£r 2 ), 72 >0. (1.6) 

The "non-resonant" set Qoq of the vectors k, for which (|1 .3j) - (jl.6p hold, is 
an extensive Cantor type set: Goo = C^ =1 G n , where {G n }^Li is a decreasing 
sequence of sets in M 2 . Each Q n has a finite number of holes in each bounded 
region. More and more holes appears when n increases, however holes added 
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at each step are of smaller and smaller size. The set Qoo satisfies the estimate: 
I (S^n Br) 



|Br| 



~R— >oo 

1 + 0(R~^), 73 >0, (1.7) 



where Br is the disk of radius R centered at the origin, | • | is the Lebesgue 
measure in R 2 . 

(3) The set V oc (X), defined as a level (isoenergetic) set for Aoo(£;), 



^oo(A) = {k£ :X 00 (k) = A} 



is proven to be a slightly distorted circle with the infinite number of holes. It 
can be described by the formula: 

^oo(A) = {k:k = ^(A, u)u, v € Boo(A)}, (1.8) 

where £>oo(A) is a subset of the unit circle S\. The set ^oo(A) can be inter- 
preted as the set of possible directions of propagation for the almost plane 
waves (|1.3p . The set £>oo(A) has a Cantor type structure and an asymptotically 
full measure on S\ as A — > oo: 

M£oo(A))= a ^oo2^ + o(a-^ 2 ), (1.9) 

here and below L(-) is the length of a curve. The value Xoo(A, V) in (jl.8j) is the 
"radius" of 2500(A) in a direction V. The function x 00 (A,z?) — A 1 / 2 describes 
the deviation of Poo (A) from the perfect circle of the radius A 1 / 2 . It is proven 
that the deviation is asymptotically small: 

Xoo(A,^)=a^oo \ 1/2 + 0{\~^), 74 >0. (1.10) 
(4) Absolute continuity of the branch of the spectrum (the semiaxis) correspond- 
ing to ^00(^,0;) is proven. 
To prove the results listed above we develop a modification of the Kolmogorov- 
Arnold-Moser (KAM) method. This paper is inspired by [18} LT9| I20j . where the 
method is used for periodic problems. In |18| KAM method is applied to classical 
Hamiltonian systems. In [TU [20] the technique developed in [18] is applied to semi- 
classical approximation for multidimensional periodic Schrodinger operators at high 
energies. 

We consider a sequence of operators 

M„ 

H = -A, # (n) = H + ^ V r, n > 1, M n -> 00 as n 00. 

r=l 

Obviously, \\H-H^\\ -> as n — > 00, where || • || is the norm in the class of bounded 
operators. Clearly, 

M n 

H (n) =H (n-l) + ^ V r . (1.11) 

We consider each operator H^ n \ n > 1, as a perturbation of the previous operator 
ij( n_1 ). Every operator is periodic, however the periods go to infinity as n — > 
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oo. We show that there is a A*, A* = A*(V), such that the semiaxis [A*,oo) is 
contained in the spectra of all operators #K For every operator there is a 

set of eigenfunctions (corresponding to the semiaxis) being close to plane waves: for 
every k in an extensive subset Q n of M 2 , there is a solution ^ n (k, x) of the differential 
equation = A n ^ n , which can be described by the formula: 

V n (k,x) = e i & g > (l + u n (k,x)) , \\un\\ Loom _ = 0(\k\~^), 71 >0, (1.12) 

where u n (k, x) has periods 2 Mn ~ 1 di, 2 M "~ 1 d2^ The corresponding eigenvalue \( n \k) 
is close to \k\ 2 : 

^ {n) (k)= lUHo0 \k\ 2 + 0{\k\-^), 72 >0. (1.13) 
The non-resonant set Q n for which (|1.13p holds, is proven to be extensive in M 2 : 

|g " nBR| = fi ->.l + 0(^). (1.14) 
|-dr| 

Estimates (|1.12p - ()1.14p are uniform in n. The set T> n {\) is defined as the level 
(isoenergetic) set for non-resonant eigenvalue \( n \k): 

V n (\) = {keg n :\^(k)=\}. 

This set is proven to be a slightly distorted circle with a finite number of holes (Fig. 
[Q[2|). The set V n {X) can be described by the formula: 

V n {\) = {k:k = x n {\, v)v, V e B n {\)}, (1.15) 

where B n (X) is a subset of the unit circle Si. The set B n {\) can be interpreted as the 
set of possible directions of propagation for almost plane waves (|1.12p . It is shown 
that {B n (\)}™ = 1 is a decreasing sequence of sets, since on each step more and more 



^Obviously, u n (k,x) is simply related to functions u r (k,x) used in (|1.5p : u n (k,x) 

5Dr="M„_! + l u r(k, x). 
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directions are excluded. Each B n {\) has an asymptotically full measure on 5i as 
A — > oo: 

L (B n (A)) = A ^oo 2vr + O (a^ 2 ) , (1.16) 

the estimate being uniform in n. The set B n has only a finite number of holes, however 
their number is growing with n. More and more holes of a smaller and smaller size 
are added at each step. The value x n (A, u) — A 1 / 2 gives the deviation of V n {\) from 
the perfect circle of the radius A 1 / 2 in the direction v. It is proven that the deviation 
is asymptotically small uniformly in n: 

x n (A, i 7) = A 1 / 2 + 0(A^), gffi*> =o(A-»), 74,75 > 0, (1.17) 

(p being an angle variable V = (cos <p, sin 93). 

On each step more and more points are excluded from the non-resonant sets Q n , 
thus {G n }^Li 1S a decreasing sequence of sets. The set £00 is defined as the limit 
set: Qoq = H^ =1 ^ n . It has the infinite number of holes at each bounded region, but 
nevertheless satisfies the relation (|1.7h . For every k £ and every n, there is a 
generalized eigenfunction of of the type ()1.12|) . It is proven that the sequence of 
^f n (k,x) has a limit in L^M 2 ) as n — > 00, when k £ Q^. The function ^^(k^x) = 
linin^oo *$> n (k,x) is a generalized eigenfunction of H. It can be written in the form 
(|1.3p - (jl.5p . Naturally, the corresponding eigenvalue Aoo(A:) is the limit of \( n \k) as 
n — > 00. 

We consider the limit i3oo(A) of B n {\): 

00 

B 00 (\)=f)B n (\), B„cB„_i. 

ra=l 

This set has a Cantor type structure on the unit circle. It is proven that Boo(X) 
has an asymptotically full measure on the unit circle (see (jl.9D ). We prove that 
the sequence x n (X,u), n = 1,2,...,, describing the isoenergetic curves V n , quickly 
converges as n — > 00. Hence, V 00 (X) can be described as the limit of V n (X) in the 
sense (|1.8p . where Xoo(A, V) = lin^^oo x n (A, V) for every V £ Boo(\). It is shown that 
the derivatives of the functions x n (A, V) (with respect to the angle variable ip on the 
unit circle) have a limit as n — > 00 for every V E ^oo(A). We denote this limit by 
9 *°£ ,ff) . Using (fTTTD we prove that 

9xoo(A, u) 



0(\-r). (1.18) 

Thus, the limit curve Poo (A) has a tangent vector in spite of its Cantor type structure, 
the tangent vector being the limit of the corresponding tangent vectors for T> n (X) as 
n — > 00. The curve Poo(A) looks as a slightly distorted circle with the infinite number 
of holes. 

Absolute continuity of the branch of the spectrum [A*(V),oo), corresponding to 
the functions , J , 00 (fc, x), k G Q^, follows from continuity properties of level curves 
00(A) with respect to A, and from convergence of spectral projections corresponding 
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to *S? n (k,x), k £ Goo, to spectral projections of H in the strong sense and uniformly 
in A, A > A*. 

The limit-periodic operator H = (—A) 1 + V, I > 6, d = 2 is considered in [24[ 125]. 
The results proved in |24[ [25] are analogous to 1-4 on pages [2j [3j The main difficulty 
of the case I = 1 comparing with I > 6 is in starting the recurrent procedure. Step 
1 here is really cumbersome comparing with the case I > 6. This technical difficulty 
is related to the fact that perturbation theory for a periodic operator (—A)' + V per 
is much simpler for I > 1 then for 1 = 1, since Bloch eigenvalues are well-spaced for 
I > 1 at high energies. 

Further steps of approximation procedure are similar to those in [24[ [25] up to some 
technical modifications. The main technical difficulty to overcome is construction of 
non-resonance sets B n (X) for every sufficiently large A, A > A*(V), the last bound 
being uniform in n. The set B n (X) is obtained by deleting a "resonant" part from 
B n -i(X). Definition of B n -\ \ B n , naturally, includes Bloch eigenvalues of i?( n_1 ). 
To describe £> n _i \ B n one needs both non-resonant eigenvalues (|1 . 13[) and resonant 
eigenvalues. No suitable formulae are known for resonant eigenvalues. Absence of 
formulae causes difficulties in estimating the size of B n \B n -\. Let us describe shortly 
how we treat this problem in the second and higher steps of approximation. Indeed, 
on n-th step of approximation we start with the operator H^ 1-1 ^. It has elemen- 
tary cell of periods Q n -i- The corresponding Bloch decomposition is denoted by 
i3 r ( n—1 )(t( n—1 )) ) quasimomentum i( n_1 ) belonging to the dual elementary cell K n -\. 
We take t^ n ~^(if) being equal to x n -i(ip) modulo K n _i, where x n _i(c^) describes 
X> n _i(A) (Fig. 1,2), i? n _i(y?) = x n _i(A,i/)z7, z7 = (cos sin y>). When *€ n _i(ip) is 
a point on 2\_i(A), vector v = (cos 92, sin <p) belongs to B n -i(X). The operator 
H'( n ~ 1 )(t( n— has a simple eigenvalue equal to A, formulas (|1.12p - (|1.17p with n — 1 
instead of n being valid. All other eigenvalues are separated from A by the distance 
much greater than ||W n _i||, see (jl.lip . It is convenient to denote such iy(™ -1 )(£(" -1 )) 
by H( n ~ 1 \iZ n ^i((p)). Next, perturbation W n of operator if( n_1 ) has bigger peri- 
ods than . We assign these bigger periods to H^ n l \ The corresponding 
Bloch decomposition we denote by ii~( n_1 ) (t^), where is quasimomentum in the 
dual elementary cell K n . According to Bloch's theory, for any the spectrum of 
union of the spectra of H^ n ^(fW + b), where b belongs to the 
lattice generated by K n in K n -i, see Fig.4, page[35j We take t^ n \f) being equal 
to x n _i(y) modulo K n . This means that t( n ~ l \tp) and t^ n \ip) satisfy the relation 
= t(»)(p) + 6* for some 6* G P( n ) where fo-Vfa) is K n -i{<p) modulo K n - X . 
It is convenient to denote such operator if( n-1 )(t( n )) by .ff( n-1 )(x n _i(y?)). Obviously, 
H^ n ~ l > {x n -i(ip)) has an eigenvalue equal to A, since i2"( n-1 )(x n _i(^) + &*) does, all 
other eigenvalues of fA n-1 )(x n _i(y>) + 6*) being separated from A. We say that <p is 
resonant if H( n ~ 1 ^ \ % n —\(ip) + bj has an eigenvalue close to A for some b 7^ 6*. Note 
that it happens if and only if the operator 

I + S n {<p) = (H^-V^^) + b) -A-e)(tf (i?n-i(¥>) + 6) +a) _1 (1.19) 
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has an eigenvalue equal to zero for some small e and b ^ b*, the operator Hq being 
the free operator with the same periods as H^ 11-1 ^. Assume for a moment we consider 
a polyharmonic operator H = (—A)' + V, I > 1. Then, S n is in the trace class and 
the determinant of / + S n exists. The resonant set B n -\ \ B n can be described (in 
terms of if) as the set of solutions of the equation Det^I + S n (<p)j = 0. To obtain 

n we take all values of e in a small interval and values of b in a finite set, 
b ^ 0. To estimate the size of B n -\ \ B n we introduce a complex angle variable (p, 
i.e., we extend our considerations to a complex neighborhood of [0, 2tt). We show 
that the determinant (11. 19ft is an analytic function of cp in <Pq, and, by this, reduce 
the problem of estimating the size of B n -\ \ B n to a problem in complex analysis. We 
use Rouche's theorem to count zeros of the determinants and to investigate how far 
the zeros move when e changes or W n is added to H^ 1 ^. It enables us to estimate 
the size of the zero set of the determinants, and, hence, the size of the non-resonance 
set <P n C which is defined as a non-zero set for the determinants. Proving that 
the non-resonance set <P n is sufficiently large, we obtain estimates (j!.16f) for B n , the 
set B n being defined by the real part of <P n . 

To obtain ^> n we delete from more and more holes of smaller and smaller radii 
at each step. Thus, the non-resonance set <S> n C has a structure of Swiss Cheese 
(Fig. HI [6]), pages [HI Ell Deleting resonance set from at each step of the recurrent 
procedure we call a "Swiss Cheese Method" . The essential difference of our method 
from those applied in similar situations before (see e.g. [18j-|21j) is that we construct 
a non-resonance set not only in the whole space of a parameter {k £ M? here), but 
also on all isoenergetic curves V n {\) in the space of parameter, corresponding to suf- 
ficiently large A. Estimates for the size of non-resonance sets on a curve require more 
subtle technical considerations ("Swiss Cheese construction") than those sufficient 
for description of a non-resonant set in the whole space of the parameter. 

When I = 1 (the present case) the determinant of I + S n ((p) (see (|1.19p ) does 
not exist, since S n ((p) is not from a trace class. We approximate S n (<p) by finite 

dimensional operators and consider the solutions of Det(l + Sn(<p)) = 0. The 
accumulation points of these solutions as N — > oo we call the solutions of "Dei" (I + 
S n (ip)) = 0. Swiss Cheese method is applied here with such a modification. 

The requirement for super exponential decay of ||V^|| as r — > oo is essential, since 
it is needed to ensure convergence of the recurrent procedure. At every step we use 
the upper bounds on to prove perturbation formulae for Bloch eigenvalues and 
eigenfunctions when A > A*(V), A* being the same for all steps. It is not important 
that potentials V r have doubling periods, in the sense that the periods of the type 
q r ~ l d\, q r ~ 1 d,2, can be treated in the same way as the doubling. 

The plan of the paper is the following. Section 2 (page [8]) is devoted to the first step 
of the recurrent procedure. Sections 3 (page [34]), 4 (page [50]) and 5 (page [59]) describe 
second, third and n-th steps of the recurrent procedure, respectively. Discussion of 
convergence of the procedure and proofs of the results 1-3, listed above, are in 
Section 6 (page [66]) . Absolute continuity is proved in Section 7 (page [69]) . Proofs 
of geometric lemmas and appendices are in Sections 8 (page [72]) and 9 (page [85]). 
respectively. 
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Further, we denote by c, C absolute constants, by C(||V||) a value depending only 
on ||V||, etc. 

2. The first approximation 

We fix < si < 10~ 4 . Let k > and be large enough so that k Sl > 500-Rch d\, cfe, , . 
We define the first operator by 

Mi 

H W .- -A + Wi, Wi :=^2V r , 2 Ml &k si . (2.1) 

r=l 

Obviously, operator has a periodic potential. We denote its periods by (ai, 0), (0, 02), 
a\ = 2 AIl ~ 1 di, a 2 = 2 Ml ~ 1 d 2 - We employ Bloch theory (see e.g. [22]) for this opera- 
tor, i.e., consider a family of operators H^'(t), t G K\, where K\ is the elementary cell 
of the dual lattice, K\ = [0, 2-na^ 1 ) x [O^vra^ 1 ). Vector t is called quasi-momentum. 
Each operator H^'(t) acts in L 2 (Qi), Q\ = [0,a\] x [0,02]. The operator H^(t) is 
defined by (|2.1I) and quasi-periodic boundary conditions: 

u(ai,x 2 ) = exp(itiai)u(0, x 2 ), u(xi,a 2 ) = exp(it 2 a 2 )u(xi,0), 

u Xl {a u x 2 ) = exp(itiai)u Xl (0,x 2 ), u X2 (xi,a 2 ) = ex.p(it 2 a 2 )u X2 (x 1 , 0). 

Each operator H^(t), t G i^i, has a discrete bounded below spectrum A^\t), 

00 

A«(t) = \J{X^(t)}, \W(t) ->• 00 as n 00. 

n=l 

The spectrum A^ 1 ) of the original operator is the union of the spectra of the 
operators ffW(t) over all t G A^ = U t6Kl AW(i). The functions \£\t) are 
continuous, so has a band structure. Extending all the eigenfunctions of H^\t) 
for all t G K\ by the quasi-periodic boundary conditions to M 2 yields a complete 
system of generalized eigenfunctions of 

Let -ffg 1 ^ be the operator (|2.ip corresponding to V^i = 0. We consider that it has 
periods a±,a 2 and that operators Hq (i) are defined in L 2 (Q±). The eigenfunctions 

of the operator Hq(£) are plane waves satisfying (|2.2p . They are naturally indexed 
by points in Z 2 : 

*$(*,a;) = \Qi\~ 1/2 expHpj(t),x), \Qi\ = a x a 2 , 

where here and below Pj(t) := (27rj 1 /a 1 + ti,2irj 2 /a 2 + ^)- The eigenvalue corre- 
sponding to ^%t,x) is equal to Pj(t), Pj(t) ■= \pj(t)\. 

Next, we introduce an isoenergetic surfaced So(A) of the free operator Hq 1 ^. A 
point t G K\ belongs to So (A) if and only if H^\t) has an eigenvalue equal to A, 
i.e., there exists j G Z 2 such that p 2 j{t) = A. This surface can be obtained as follows: 

the circle of radius k = vA centered at the origin is divided into pieces by the dual 
lattice {Pq(0)} q <=z2 , and then all pieces are translated in a parallel manner into the 

2 We write a(k) « b(k) when the inequalities |6(fc) < a(k) < 2b(fc) hold. 
3 "surface" is a traditional term. In our case, it is a curve. 
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Figure 3. The isoenergetic surface So (A) of the free operator 

cell K\ of the dual lattice. We also can get So (A) by drawing sufficiently many circles 
of radii k centered at the dual lattice {p q (0)} q ^z 2 an d by looking at the figure in 
the cell K\. As the result of either of these two procedures, we obtain a circle of 
radius k "packed into the bag K\" as is shown in the Fig. [3l Note that each piece 
of iSo(A) can be described by an equation = A for a fixed j. If t G So (A), then j 
can be uniquely defined from the last equation, unless t is a point of self-intersection 
of the isoenergetic surface. A point t is a self-intersection of So (A) if and only if 
Pq(t) = p 2 (t) = A for at least one pair of indices q, j, q ^ j. 

Note that any vector x in R 2 can be uniquely represented in the form ok = Pj(t), 
where j G Z 2 and t G K\. Let K,\ be the parallel shift into K\\ 

JCi:R 2 ^K l7 JC 1 p j {t)=t. 

Suppose Q, C M 2 . In order to obtain /Ci$7, it is necessary to partition Q by the lattice 
with nodes at the points p q {0), g G Z 2 and to shift all parts in a parallel manner 
into a single cell. It is obvious that |/CiSl| < for any f2. If fi is a smooth curve, 
then L{K,\£l) < L(Q). For any pair of sets Q,\ and r^, /Cif^i U f^) = /Ci$7i U /C1Q2- 
Obviously, K,\Sk = S${\) and ^(^(A)) = L(Sk) = 2-7rk, k = \/A, Sk being the circle 
of radius k centered at the origin. 

The operator H^\t), t G K\, has the following matrix representation in the basis 
of plane waves *&j(t,x), j G Z 2 : 

H {1 \t) mg = p 2 m {t)5 mq + w m - q , m,q€Z 2 . 
Here and below, 5 mq is the Kronecker symbol, w m - q are Fourier coefficients of W\\ 

W Q = T7TT / Wi(x) exp -i(p q , x)dx, w q = W- q . 
\Ql\ JQi 
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Using assumptions on potentials V r , we easily obtain: 

wq = 0, \w q \ < ||Wi||, w q = 0, when p q > Rq. (2-3) 

The matrix H^(t) mq describes an operator in /2(Z 2 ) unitary equivalent to H^\t) 
in L,2(Qi). From now on, we denote the operator in ^(Z 2 ) also by H^\t). Since 
the canonical basis in ^(Z 2 ) does not depend on t and all dependence on t is in the 
matrix, the matrix H^ l \t) mq and, hence, the operator H^ l \t) : ^(Z 2 ) —> ^(Z 2 ) can 
be analytically extended in t from K\ to C 2 . From now on, when we refer to H^(t) 
for t £ C 2 , we mean the operator in ^(Z 2 ). 

2.1. Perturbation theory. We formulate the main results of the perturbation the- 
ory considering H^ l \t) as a perturbed operator of the free operator i.e., we 
construct the perturbation series for eigenvalues and spectral projections for t in a 
small neighborhood of a non-resonant set xi- 

Lemma 2.1 (Geometric Lemma). For arbitrarily small positive 5 < si and j3 : 
4si < 2/3 < 1 — 15si — 86 and sufficiently large A, A > Ao(/3, 8%, 5), there exists a 
non-resonant set xi(A, f3, si, 5) C Sb(A) satisfying 

(i) For any t G Xi(^j P: s ii $)> there exists a unique j G Z 2 such that Pj(t) = k, 
k := vA- The following inequalities hold: 

4 mm . \p 2 (t) -p)(t)\\p 2 +q {t) -p 2 (t)\ > k 2 \ (2.4) 

0<p q <k s i 



2 min \p 2 +q {t) - p 2 {t)\ > k 1 -^- 5 . (2.5) 

(ii) For any t in the complex (k 2l3 ~ 2 ~ Sl ~ 2S )-neighborhood of the non-resonance set 
in C 2 , there exists a unique j G Z 2 such that 

\p 2 (t) - k 2 \ < Zk 2 ?- 1 -^- 25 . (2.6) 

Estimates (H3D and ([23)1 hold. 

(iii) The non-resonance set xi(A, f3, s±, 5) has an asymptotically full measure on 
So (A) in the sense of 

L(S (\)\xiM,s lt 6)) = Q{k ^ (27) 

L (60(A)) A->oo 

Corollary 2.2. If t belongs to the complex (k 2/3 ~ 2 ~ Sl ~ 2S )-neighborhood of the non- 
resonant set xi(A, j3, s\, 5) in C 2 , then for any z G C lying on the circle 

d = {z : \z - k 2 \ = k 2 ^ 1 - 31 - 5 } (2.8) 
and for any i and q in Z 2 with < p q < k Sl 

2\p 2 (t) -z\> k 2 ^ 1 -^- 5 (2.9) 

16|p?(t) - z\\p 2 +q (t) -z\> k 2 ^- 2S (2.10) 
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The set xi(A, /3, si, 5) is defined below. The lemma and the corollary are proven 
in Section EU Section 18.1.11 is introductory. Properties (i), (ii) and Corollary 12.21 
are proven in Section 18.1.21 (Corollary I8.5p . Section 18.1.31 is devoted to the proof 
of (iii). Here we just note the following. An analog of the lemma and corollary is 
proven earlier in [23] . However, in this work we need more detailed description of 
Xi(A, (3, si,8) than that in [23J. Moreover, we need a complex analog of xi(A, /3, s\,5). 

We define xi(A, f3, si, 5) by the formula: 

Xi(A,/3,si,(5) = K,iV (X) nonres , (2.11) 
where T>Q(X) nonres is a subset of the sphere 

Do{\) nonres = {kv, u = (cosy?, sin ip), ip 6 Si}, (2.12) 
the set 0i being defined as the real part of a set <P\ C C: 

Gi = <f>a n [0,2vr), (2.13) 

= # \0 W , (2.14) 
where <Pq is the complex k~ s neighborhood of [0, 2ir), 

^ = {^eC:%e [0, 2vr), \$s<p\ < AT 5 }, (2.15) 
is the union of discs 0^: 

= U meZ 2 i <p m <4fc, ± Cm- ( 2 - 16 ) 

Discs O m are centered at the points ip^, which are solutions in C (two for each m) 
of the equations: 

\k(tp) +Pm\l = k 2 , k((p) = k(cos ip, sin ip), (2.17) 

here and below \x\1 = x\ + x\ for any x = {x\,xz), X\,X2 G C. The radii r m of the 
discs are chosen to ensure the properties of the set Xl(A, P, si, 6) to be valid. They 
are given by the formula (|8.4p . 

Let Ej(t) be the spectral projection of the free operator corresponding to the eigen- 
value Pj(t), (Ej) rm = Sj r Sj m . In the (/c 2/3 ~ 2 ~ Sl ~ 25 )-neighborhood of xi(A, /3, si, S), 
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we define functions gj- (k, t) and operator-valued functions Gf- (k, t), r = 1, 2, . . . as 
follows. 

9?Hk,t) = ^Tr I - z^WtYdz, (2.18) 

G?\k,t) = I ((H< 1 \t)-z)- 1 W 1 r(H!P{t)-z)- 1 dz. (2.19) 

To find gr l \k,t) and G^\k,t), it is necessary to compute the residues of a rational 
function of a simple structure, whose numerator does not depend on z, while the 
denominator is a product of factors of the type (p 2 (t) — z). For all t in the non- 
resonance set the integrand has a single pole within C\ at the point z = k 2 = p?(i). 

By computing the residue at this point, we obtain explicit expressions for g^\k,t) 
and G^ 1 \k,t). For example, g^\k,t) = 0, 

g ( i\k,t)= Y, \ Wq \ 2 (v){t)-p) +q (t)Y l 



8 Jr^^ (<) -^ (t))( ^ (t) -^ (t)): 



(2.20) 



G i 1)(M)rm " p j(t)-pl(t) 6rj + pj^ r - p 2(t) 6mi > G i ] ( k ^n-°- ( 2 - 21 ) 

For technical reasons, we introduce parameter a in front of the potential W\, Ha^ = 
(—A)' + aW±, < a < 1. We denote the operator i/^ 1 ^ with a = 1 simply by ijW. 

Theorem 2.3. Suppose t belongs to the (k 2l3 ~ 2 ~ Sl ~ 2S ) -neighborhood in K\ of the 
non-resonant set xi(A, /3, si, <5) ; < 5 < si, 8si + 65 < 2j3 < 1 — 15si — 85. Then 
for sufficiently large k, k > fco(|| W\ ||, s\, /3, 5) and for all a, —1 < a < 1, i/iere 
exists a single eigenvalue of the operator Ha\t) in the interval E\(k,5) := (k 2 — 

£,2/3-l-si-<5^2 _|_ fc2f3-l-s 1 -8y j£ -g gj^ ven fry tfr e ser { es 

oo 

\f ) (a,t)=p 2 1 (t) + J2^ r gi 1) (k,t), (2.22) 

r=2 

converging absolutely, where the index j is defined as in Lemma \2.l\ The spectral 
projection corresponding to X^p (a, t) is given by the series 



j 

r=l 

which converges in the trace class Si. 



Ef\a,t) = Ej+Y^cTGMfot), (2.23) 



Moreover, g^{k,t) and Gr(k,t) satisfy the estimates: 

\g?Hk,t)\ < A ; 2 ^ 1 - Sl - 5 ( C ||^ 1 ||A;-/ 3+4si+25 ) r , (2.24) 

\\GW(k,t)\\ 1 <r{ck-^ +2S ) r , (2.25) 
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|| • ||i being the norm in the trace class Si. In addition, 

\g^\k,t)\ < c||Wi|| 2 £;~ 2+10si+45 , \g^(k,t)\ < c|| ^ || 3 fc- 2+10si+4 ' 5 , (2.26) 

g { 2 ] {k,t) > 0, when \\Wi\\ £ 0, (2.27) 

Operators G*j? (k,t) are finite dimensional: 

G^(k,t)u = when \i - j\ > rR or \l - j\ > rR . (2.28) 

The series \2. 22\) and \2. 23\) converge uniformly with respect to a in the complex disk 
\a\ < 1. 

Corollary 2.4. The perturbed eigenvalue and its spectral projection satisfy 

\\f\a,t)-pj(t)\ < a 2 C(||Wi||)A;- 1 - 2 ^ +15 * l+lw , (2.29) 

\\Ef\a,t) -Ej\\ x < c|a|||^||Ar /3+4si + 35 . (2.30) 

Proof. The proof of the theorem is based on expanding the resolvent 
(Ha\t) — z)^ 1 in a perturbation series for z belonging to the contour C\ about the 
unperturbed eigenvalue p 2 (t). It is completely analogous to the proofs of Theorems 



2.1 and 3.1 in |23j . Indeed, it is obvious that 

(ff«(i) - z)- 1 = (H^(t) - z)- l '\l - aA.r^H^it) - z)- l /\ (2.31) 

A x = A x {z,t) := -(H^(t) - zy^WxiH^it) - z)' 1 / 2 . (2.32) 
It follows from Corollary 12.21 estimates (I2.3P and i?o < k Sl , that 

\\(H^(t) - z)- 1 ]] < 2£r 2 ^+ 1+Sl + 5 , z e Ci, (2.33) 

Pill < 16||Wi||AT /,+4 ' 1+M . (2.34) 



Thus, H-Aill < 1 for sufficiently large k > feo(||Wi||, f3, si, 6). Expanding (I—aA\) 1 
in powers of aA±, we obtain 

oo 

(^w-^^-^w-^-^E^^w-^)" 172 ^^^)-^)" 172 - ( 2 - 35 ) 

r=l 

Note that (Hq (t) — z)^ 1 € S2. Taking into account estimates (|2.34[) . we see that 
the series (|2.35p converges in the class Si uniformly with respect to a in the whole 
complex disk |a| < 1. For real as we substitute the series into the following formula 
for a spectral projection 

Ef(a,t) = -^-J (H^^-z^dz. (2.36) 

Integrating termwise, we arrive at (|2.23p . Next, we prove estimate (|2.25p for Gr (k, t). 
Indeed, it is easy to see that 

G^CM) = ~ f (H^\t) - z)-V 2 A\{HW(t) - zyWdz. (2.37) 
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We introduce the operator Aq = (I — Ej)A\{I — Ej). It is obvious that \\Aq\\ < \\Ai\\. 
In addition 

/ - z)-^ 2 A r (H^\t) - z)- l ' 2 dz = 0, (2.38) 

since the integrand is holomorphic inside the circle. Thus, 

GW(M) = / (^o (1) (*) " ^)" 1/2 K " A r )(H^(t) - x)~Wdz. (2.39) 

Since A 1 - A = EjAi(I - Ej) + (I - Ej)AiEj, EjAiEj = and is a one- 
dimensional projection, we get 

< (2.40) 

Since 

\\A\ - 45 ||i < r||ii - Aollipif- 1 , (2.41) 
we obtain from relation (|2.34p that 

\\A\ - Al\\x < r (c||W 1 ||Ar /3+4si+25 Y . (2.42) 

Noting that the length of C\ is equal to 2irk 2 P~ 1 ~ Sl ~ & , we obtain from formula (12.391) : 

||Gf«(fc,t)||i < k^- 1 - s ^ 5 \\(H^ 1 \t) - zr^fU^ - A^U (2.43) 

Using inequalities (|2.33p and (|2.42p . we get (|2.25p . Convergence of the series in the 
complex disk |a| < 1 easily follows. Note that (A\)u = if \i — l\ > Ro- Now it is 
not difficult to check (|2.28p . for details see |23j . 

We show now that operator Ha \ — 1 < a < 1, has a single eigenvalue (n a = 1) 
in the interval e(k, S). It is easy to see that n a is determined from the formula 
n a — n o = Tr(£'j 1 ' 1 — Ej). Considering formula (12.23p . we obtain 

oo 

|n a -T*|<£[|G«(M)[|i = °(l). 

r=l 

Since n a and no are integers, and no = 1 by the hypothesis of the theorem, the 
operator Ha\t) for all a, — 1 < a < 1, has a single eigenvalue in e(k,5). 
Further, we use the well-known formula: 



^^ = -J-Tr/ Wl{ H^ { t)-z)^ d z 



da 2iri Jq 



(2.44) 



Ci 



Using formula (pT35D and considering that Tr( W\ (H^ (i) - z)' 1 ) = 0, we obtain: 



da 



Y^ra^gV^t). (2.45) 



r=2 

Integrating the last relation with respect to a and noting that A^ 1 ) (0, t) = p 2 (t), we get 

formula (|2.22p . To prove estimate (|2.24p we note that gr = — (27rir)~ lr Tr § c A[dz, 
and, therefore, 

|#(M) l^r- 1 ^- 1 - 1 "'^!-^!!!, (2.46) 
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Using estimate (|2.42p . we obtain inequality (I2.24f) . Estimates (|2.26j) are better than 
those provided for gi, 53 by the general estimate (|2.24p . We obtain estimate (|2.26p 
for g2(k,t) using the second part of formula (I2.20|) and the estimate (|2.5p . Estimate 
(|2.27p easily follows from the same formulas. The estimate for g%(k, t) can be obtained 
in the analogous way. For details, see [23], Theorem 3.2. | 

The series (|2.22p and (|2.23|) can be extended as holomorphic functions of t in a 
complex neighborhood of xi(fc, ft, s\, 5). They can be differentiated with respect to 
t any number of times with retaining their asymptotic character. Before presenting 
exact statements, we introduce the following notations: 

Q\m\ 

m = (1711,1712), \m\ := mi + m.2, ml := milmj!. 

Lemma 2.5. Coefficients gr(k,t), and G^ (k,t) can be continued as holomor- 
phic functions of two variables from the real (k 2 P~ 2 ~ Sl ~ 2S ) -neighborhood of the non- 
resonance set xi(k, f3,si,5) to its complex (k 2/3 ~ 2 ~ Sl ~ 2S )-neighborhood. In this com- 
plex neighborhood the following estimates hold. 

\T(m)gV>(k,t)\ < m!^- 1 - Sl - <5 ( C ||VF 1 ||A : ^+ 4si+25 ) r A ; l m l( 2 - 2 ' 3 + Sl+2<5 ), (2.47) 
||r(m)G«(M)ll < m!r(c||^ 1 ||A : -' 3+4si+25 ) r A ; H(2-2^+ Sl+ 25)_ (2>4g) 

The coefficients g r (k,t), r < fc^i?^ 1 can be continued to the complex (k~ 3si ~ 25 )- 
neighborhood of Xi(k, ft, s%, 5) . They obey the estimates: 

\T(m)g ( 2 1) (k,t)\ < cm!||W l || 2 A;- 2+10si+46 fcl m l( 3si+25 ), (2.49) 

\T(m)g^ {k,t)\ < cm!||M/i|| 3 fc- 2+10si+4,5 A;l m l( 3si+2 ' 5 ). (2.50) 
\T{m)g^{k,t)\ < m!A; 2 ^- 1 - sl - 5 (c||Wi||A;-^ +4si+2<5 ) r A;l w K3si+25) ) 4 < r < k Sl R \ 

(2.51) 

Proof. Coefficients g? (k, t), and G> (k, t) can be extended analytically to the (k 2l3 ~ 2 ~ s - 
neighborhood of xi(fc, ft, si, S), since estimates (|2.9p . (|2.10p are stable in this neigh- 
borhood. The coefficients g r (k,t), r < fc Sl i?Q 1 , can be extended to a bigger neighbor- 
hood, since only (|2.5p is required to estimate them. The estimates for the derivatives 
are obtained by means of Cauchy integrals. For details, see Theorem 3.3 and Corol- 
lary 3.3 in [23]. ■ 

Lemma 12.51 implies the following theorem. 

Theorem 2.6. The series (|2.22p and ([2.23P can be continued as holomorphic func- 
tions of two variables from the real (k 2l3 ~ 2 ~ Sl ~ 2S ) -neighborhood of the non-resonance 
set Xi(k, ft, s\, 5) to its complex (k 2l3 ~ 2 ~ Sl ~ 2S )-neighborhood. The following estimates 
hold in the complex neighborhood: 

\T{m){\f\a,t)-p 2 {t))\ < m!C '( Wl ) a 2 A .-l-2^ + 15, 1 +115+|n l |(2-2^+ S1 +25) ) (2>52) 

\\T(m)(E ( j 1) {a,t) - Ej)\\ < cm!||VFi||afc- /3+4si+3<5+|m|(2 - 2/3+Sl+2,5)) . (2.53) 
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There are stronger estimates for \m\ = 1,2: 

VX ( p(a,t) - 2pj{t) < CiW^a 2 ^ 2 ^ 183 ^ 135 , (2.54) 

T(m) (xf\a,t)-^\\ < C'(Wi)a 2 A;- 1 -^ +21si+15 ' 5 , if\m\ = 2. (2.55) 

The next lemma is used in the second step of approximation, where the operator 
(t) plays a role of the initial (unperturbed) operator. 

Lemma 2.7. For any z on the circle C\ given in (12. 8h and t in the (/j 2 / 3 -2- £, i- 2<5 )- 
neighborhood ofxi(k,5) in K±, 

\\(H {1 \t) - 1 1| < 2k~ w+l+Sl+s . (2.56) 

Proof. The estimate follows from fl2~3il . (l2~33D and (|2~34]l . | 

2.2. Non-resonant part of the Isoenergetic Set of H^K In this subsection we 
choose (3 to have the biggest possible value (3 = /3q, 2(3q = 1 — 15si — 95. Let Si (A) 
be the isoenergetic surface of the perturbed operator Ha\ i.e., 

Si (A) = {t G K x : 3n G N s.t. A^cM) = A}, (2.57) 

where {A« (a, t)}^ =1 is the complete set of eigenvalues of Ha\t). We construct 
a "non-resonance" subset x*(A) of Si (A), which corresponds to the non-resonance 

eigenvalues X^\a,t) given by the perturbation series. By Lemma 12.11 for every 
t belonging to the non-resonant set xi(A, /3, si, 5), there is a single j G Z 2 such 
that Pj(i) = k, k = \/A. This means that formula (|2.1ip establishes a one-to-one 
correspondence between xi(A, /3, si, 5) and T>o(X) nonres . Let x G TJo(X) nonres and 
j G Z 2 , t G xi(^) s i> are defined by the relation x = Pj(t)a According to 
Theorem l2.3l for sufficiently large k, there exists an eigenvalue of the operator Ha\t), 
t = JCix, < a < 1, given by (12.22H . It is convenient here to denote AJ- (a, i) by 

A^^(a, x); we can do this since there is a one-to-one correspondence between x and 
the pair (t, j). We rewrite (|2.22p in the form 

oo 

A«(a,x) = x 2 + /i(a,x), x=|x|, /i (a, x) = £ oTg « (x) , (2.58) 

r=2 

<7r being defined by (|2.18p with j and t such that Pj(t) = x. By Theorem 12.61 
/i(a, x) satisfies the following estimates when 2/3 = 1 — 15si — 95: 

\fi(a, x)| < a 2 C'(Wi)*r 2+30si+20 ' 5 , (2.59) 
|T(ro)/i(a, x)| < a 2 c(Wl)A .-2+30, 1 +205 + | m |(i+i6 Sl+ ii5)_ (2 _ 6Q) 

Estimates (ET54"D . (|235|) yield: 

V/i (a, H) = (V 2+33si + 225 ) ) (2 gl) 



Usually the vector pj(t) is denoted by k, the corresponding plane wave being e {k ' x) . We use the 
less common notation iJ, since we already have other k's in the text. 
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VA« (a, x) = 2x + O (V 2+33si + 225 ) , (2.62) 
T(m)fi(a, k) = (V 2+36si+245 ) , T(m)\^\a, x) = 0(1), if |m| = 2. (2.63) 



By Theorem 12. 6\ the series (I2.58D converges in the (k 1 16si 115 )-neighborhood of 
V (\) non res, the estimate flg^gp - flglEgl hold. 

Let us recall that operator H^\t) is defined in Z 2 , i.e., just by the matrix depending 
on i. We define H^^H) by substituting x instead of Pj(t) into H^(t). For any 
i £ Xi(k, /3q, si, 5) matrices H^(t) and x) are the same, up to the shift of indices 
by j. Further we consider (x) for x = xv, when (p is in the complex 2A;~ 2 ~ 16si ~ 11<5 
neighborhood of <P\ and x is in the complex 2/c _1_16si_11<5 neighborhood of k. Using 
Lemma [8.91 we easily obtain the estimates analogous to (|2.33p . (|2.34[) . which provide 
the convergence of the series for the resolvent. Operator H^ l \x) is not self-adjoint 
for complex <p and we do not try to make spectral analysis of it. However, we notice 
that the series (|2.58p converges when cp is in the complex 2/c~ 2_16si_11<5 neighborhood 
of <Pi and x is in the complex 2/c~ 1_16si_11<5 neighborhood of k. Lemma 12.71 admits 
the following generalization. 

Lemma 2.8. For any z on the circle C\ given in (|2.8p and x = xz7 ; (p being in 
the complex 2/c~ 2 ~ 16si ~ 11<5 neighborhood of<Pi, x being in the complex 2/c -1 ~ 16si ~ 11<5 
neighborhood of k, 

||(tf (1) (x) - 1 1| < 2£T 2/m+Sl+5 . (2.64) 

Let -Si (A) be a set of unit vectors corresponding to T>o(\) nonres : 

Bi(X) = {v e S\ : ki> e T> (X) nonres } = {i? £ Si : V = (cos (p, sin >p), <p £ ©i}. 

It is easy to see that £>i(A) is a unit circle with holes, centered at the origin. Since 
formulas (l23Kl) - (|2ll3l ) hold in the (A;- 1 ^ 16si ~ 11 ' 5 )-neighborhood of V (X) nonres , they 
hold for any xz7 such that v £ B\(\), \x-k\ < /j- 1 - 16 *!- 115 . We define Pi (A) as the 
level set of the function \^\a, x) in this neighborhood: 

Pi (A) := {i? = xz7 : v £ 0i(A), \x - k\ < k~ l - l(isi - lis , A (1) (a, x) = A}. (2.65) 

We prove in Lemma [2.101 that Pi (A) is a distorted circle with holes, which is close to 
the circle of radius k; see Fig. [TJ First, we prove that the equation \^(a, x) = A is 
solvable with respect to x = |x| for any u = ^ £ Z?i(A). 

Lemma 2.9. For every V £ £>i(A) and every a, < a < 1, and sufficiently large X, 
there is a unique x\ = xi(A, z?) in the interval 

h : = [ k - jfe-1-16.1-11^ k + ^-1-16^-11^ k 2 = Aj 

such that 

\M(a,x 1 v) = \. (2.66) 
Furthermore, |xi - fc| < C(Wi)k- 3+30si+20S . 

Proof. Formula (|2.62p yields 

-2k, (2.67) 

ax 
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when |x — k\ < k 1 16si 11S . Now the lemma easily follows from (I2.59|) . For details 
see Lemma 2.10 in [24]. ■ 



Let us introduce new notations. Let Si be the (k 2 16si 11<5 )-neighborhood of 3>i. 
Note that (f238]) - ([2T63]) hold for cp in Si and even its (£r 2 " 16si - 11,5 )-neighborhood 
when x = xz7, |x— &| < k~ l ~ 16si ~ lis . Let 

xi (<£>):= xi( A, z?), = (cos sin <p), x\ (</?) = >ci(ip)u, hi{tp) = xi(ip) — k, (2.68) 

xi(A, i?) being defined by Lemma 12.91 

Lemma 2.10. (1) For sufficiently large A, the set 2?i(A) is a distorted circle with 
holes; it can be described by the formula 

2?! (A) = {x G M 2 : x = xi(tp)P, ^eBi(A)}, (2.69) 

where xi(<p) = k + hi{ip) and hi((f) obeys the inequalities 



< C{Wi)k- 1+33si+22S . (2.70) 



dip 

(2) The total length of B\{\) satisfies the estimate 

L(Bi) = 2vr(l + 0(k- 5 )). (2.71) 

(3) The function h\ (ip) can be extended as a holomorphic function of (p from O2 
to Si. Estimates (I2.70P hold in Si too. 

(4) The curve £>i(A) has a length which is asymptotically close to that of the whole 
circle in the sense that 

L(Vi(X)) = 2irk(l + 0{k- 5 )), A = k 2 . (2.72) 

A— >oo 



Proof. Here are the main points of the proof, for details see Lemma 2.11 in 

(1) Inequalities (f!T70|) easily follow from (j238|) - (pT60|) and the definition of 
T>i(X), Implicit function theorem being applied. 

(2) By definition, £>i(A) is the set of directions corresponding to T>o(\) nonres , the 
latter set being a subset of the sphere of radius k. Formula (|2.1ip estab- 
lishes a one-to-one correspondence between Xi(^) A S) and T>o(\) nonres , 
their lengths being equal. Considering (12. 7|) . we obtain L (T>o(\) nonres ) = 
2irk(l + 0(k- 5 )). Hence, (j277T|) holds. 

(3) The series (|2.58p converges for x = xz?, when cp in Si and even its (&- 2 - 16s i- 11<5 ) 
neighborhood and x G to, to being the complex (/c _1_16si_11<5 )-neighborhood 
of k, to = {x G C : |x - fe| < fe-i-ieai-Htf^ Function AW(a,x) is an- 
alytic in x and ip. The estimate (|2.59|) holds. Thus, \\^(a, ku) — x 2 | < 
C(Wi)k~ 2+30si+20S , and |x 2 - /c 2 | ps jfe-i6«i-iW when x G do;. Applying 
Rouche's theorem in w, we obtain that the equations A^^(xz^) = A; 2 and 
x 2 = k 2 have the same number of solutions x in to for every ip £ Si. Obvi- 
ously, x 2 = k 2 has just one solution. Therefore, X^(a, >cv) = k 2 also has just 
one solution xi(<^). Using (|2,67p and Implicit function Theorem, we obtain 
that xi(y?) is an analytic function of ip in and estimates (I2.70P hold. 

(4) Estimate (l2T72jl follows from rf2770j) and (l2?7ll . 
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Next, we define the non-resonance subset X*(A) of the isoenergetic set S\(X) as the 
parallel shift of Pi (A) into K±: 

X * 1 (A):=/C 1 P 1 (A). (2.73) 

Lemma 2.11. The set x*(A) belongs to the (k~ 3+30si+20S ) -neighborhood o/xi(A) 

in K\. If t G X*(A); then the operator Ha if) has a simple eigenvalue \„ (a,i), 
n G N, equal to X, no other eigenvalues being in the interval £i(k,5), £\{k,5) := 
(k 2 — k~ 1Qsi ~ lls , k 2 + k~ 1Qsi ~ lls ) . This eigenvalue is given by the perturbation series 
(|2.22ft . where j is uniquely defined by t from the relation p 2 (t) G £\{k, 5). 

Proof. By LemmaEJin T>i(X) is in the (/c- 3+30si+205 )-neighborhood of P (A). Since 
Xi(A) = tCiVo(X) and x*(A) = K\D\(X), we immediately obtain that X*(A) is in the 
(cfc~ 3+30si+20l5 )-neighborhood of xi(A). The size of this neighborhood is less than 
k-i-ltoi-iw here _i _ 16si _ u $ = 2 /3 - 2 - si - 2<5 for 2/3 = 1 - 15si - 95. Hence, 

Theorem 12.31 holds for any t G X*(A): there is a single eigenvalue of Ha\t) in the 
interval £i(k, 6). Since Xi(A) C Si (A), this eigenvalue is equal to A. By the theorem, 
the eigenvalue is given by the series (|2.22j) . j being uniquely defined by t from the 
relation p 2 (t) G £i(k, 5). | 

Lemma 2.12. Formula (|2.73p establishes a one-to-one correspondence betweenXi(X) 
and Pi (A). 

Remark 1. From the geometric point of view, this means that Xi(A) does not have 
self-intersections. 

Proof. Suppose there is a pair xi,i,><i,2 G Pi(A), such that /Cii?i,i = 
= t, t G Xi(A)- By the definition (|2.65p of Pi (A), we have A^^(a, i?i,i) = 
2) = A, i.e., the eigenvalue A of Ha\t) is not simple. This contradicts the 
previous lemma. | 



2.3. Preparation for the Next Approximation. In the next steps of approxi- 
mations we will need estimates for the resolvent ^ilW (x.\{Lp) + 6) 1 — k 2 ^j , b G 
^1 \ {0}' </> G ^1. Let 60 be the distance of the point b to the nearest corner of K%: 

bn = min \b — 2irm/a\. (2.74) 

m=(0,0),(Q,l),(l,0),(l,l) 

We will consider separately two cases: bo > ^-i- 16 ^- 125 and < 60 < A; -1-1651-125 . 
Our goal is to prove Lemmas 12.291 and 12.341 which give the estimates for the resol- 
vent on a subset of $1 for the cases b > /c _1 ~ 16si ~ 12<5 and < b Q < fc-i-iesi-l^ 
respectively. In this subsection we choose f3 to be relatively small, lOOsi < (3 < 
1/12 - 28si - U5. We denote such by fa. 
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2.3.1. The case b > fe-i-i6«i-i2<y_ The goal is to cons t ruc t a set O^'ib) such that 
Lemma 12.291 holds. The general scheme of considerations is similar to that for the 
case H = (—A)' + V, I > 6, see [23], Section 3.5.1. First, we describe shortly the 
scheme in [24J. Further, we explain what changes have to be made to adjust it for 
the present case. Indeed, in [23] we constructed a set 0(b), which consists of small 
discs centered at the poles of the resolvent 

H^(k(ip)+b)-k 2 y\ (2.75) 

However, the size of discs was chosen sufficiently large to ensure that 

-l 

< fi>0, when <p g 0{b). 



H^ik^ + b) 



(2.76) 



It was shown that the number of poles of (|2.75p in (Pq does not exceed cok 2+2si . Using 
perturbation arguments we showed that the estimate analogous to (12.75D holds for 
the perturbed resolvent 

(2.77) 



in <Z>i \ 0(b): 



< £T M , fi > 0. 



(2.78) 



Next, suppose O c (b) is a connected component of 0(b), which is entirely in <Pi: 
O c (b) C $\. Then (I2.77P is an analytic function in O c (b). Using determinants and 
Rouche's theorem, we showed that (12.751) and (]2.77|) have the same number of poles 
inside each connected component of 0(b). Therefore, the total number of poles of 
(|2.77|) in such components does not exceed cok 2+2si . Next, we replaced the discs in 
0(b) by much smaller discs around the poles of (|2.77j) . By doing this, we obtained 
a much smaller set O s (b), O s (b) C 0(b). Using analyticity of (j2.77j) inside 0(b) and 
the maximum principle, we arrived at the estimate: 

./ 



11^(2^) + b)-k' 



< ck~ 



Rb 



J = cok 



2+2si 



(2.79) 



where Rf, is the size of the biggest component of 0(b) and is the radius of the 
discs constituting O s (b). Thus, we shrank the set 0(b), but still obtained a bound 
for the norm of the resolvent, which is essential for the next step of approximation in 
[24|. It may happen, however, that O c (b) n <P\ / 0, but O c (b) <£_ ®\. In this case we 
showed that the radii of the discs in O c (b) were sufficiently small to make sure that 
the function xi(<p) can be analytically extended from to the interior each O c (b) 
when O c (b) fl^i ^ even if O c (b) <f_ <P\. For / > 6 this condition on radii does not 
contradict to (|2.76p . 

In the present case I = 1 the plan is basically the same. However, we have essential 
technical complications. The set 0(b), constructed for I > 6, is now too small to 
provide convergence of perturbation series for (12.771) when ip 0(b). Therefore, we 
need to construct a bigger set O^ (b), such that the series converges when <p $ O^ (b). 
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It can be constructed by analogy with (|2.16p . Let Oc (b) be a connected component 
of C>W(6). if { c\b) C <2>i, then considering as in [24], we obtain estimates similar 
to (127511 . A difficult case is oP(b) D #i ^ 0, but 0^ 1} (6) £ <f>i- On one hand, 
the function xi(ip), included in (12.160 . is defined on the set which has holes 
as small as e = fc~ 2 ~ 12si ~ 9<5 . These holes are essential in the construction of <Pi. 
Hence, the function H\{ip) can be analytically extended no further than cs, c < 1, 
neighborhood of <P\. Thus, to extend >?i(</?) analytically into a connected component 
0^\b) of a set such that 0^\b) <£. we have to make sure that the size 

of 0^\b) is smaller than e. On the other hand, deleting only such a small set from 
<Pi cannot provide convergence of perturbations series for (pTfTjl in <2>i \ £>M(6). The 
properties of the smallest set 0^\b), which we can construct to provide convergence 
of the perturbation series, is given in Lemma l2.161 Its connected component still can 
have a size up to / c ~ 11 / 6 +/ 3 i+ 12s i+ 4(5 ) a t least this is the strongest estimate we can 
prove. Obviously, the size of such a component is still greater then the size e of the 
neighborhood of where >ci(<p) is holomorphic. Thus, we can not construct the set 
£>M(6) which fits into e-neighborhood of <P\, even though it is much smaller then 
(|2.16p . To overcome this difficulty we need additional considerations. We represent 

ok\ in the form i?i(y?) = i?*(y?) + ( k~ ck3/A ), where >?*(</?) is an analytic function in a 



set ^>o larger then , 5>o 3 ^i > see Lemma 12.211 The smallest hole in <£>o has the size 
/c~ 7 / 4-12si-9<5 . Since the size of each connected component of 0^\b) is smaller than 
fc -n/6+/3i+i2s 1 +4<5 ; ig holomorphic in every oP(b), such that (b) n <2>i ± 0. 

We prove an estimate for the norm of 

H^iM^ + ^-k 2 )' 1 (2.80) 



outside { c\b). Since, the set 0^'(b) is rather small, the proof of an estimate for 
the norm is quite technical (Lemmas 12.161 12.201 12.231 12.25j) . Further we construct a 

set Os 1 " l (b) taking small discs around the poles of (|2.80p . Using the analyticity of the 
resolvent and the maximum principle, we obtain (Lemma I2.25P an estimate of the 
type (12791) for (123D1) : 

+ b)-k 2 y 1 \\ <ck»* (^) J , J = c k 2+2s \, M *>0. (2.81) 

The size of each connected component of 0^\b) is much less than fc- 2 - 125 !- 9 " 5 . 
Therefore, the function x\{^p) can be analytically extended into interior of every 
connected component of (b) intersecting with <&i. Using the estimate i?i(ty?) = 
x*(</?) + O (k~ ck3/4 ^ in the e neighborhood of 3>i, we prove that (|2.77|) obeys an 

estimate similar to (|2.8ip , see Lemma 12.291 This lemma is used in the second and 
step of approximation. 

One more technical difficulty is related to the fact that the resolvent (|2.75p is not 
from the trace class when 1 = 1, while it is from the trace class when I > 1 . This means 
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that the determinant of the operator ( H^> (k(ip) + b) — k 2 J yHfi ' (k(tp) +6) — k' 

which we considered in |24] as a complex function of ip for I > 6, is not defined for the 
present case 1 = 1. To overcome this difficulty, we take a family of finite dimensional 
projections P/v, P/v — > s I as N goes to infinity. We consider a finite dimensional oper- 
ator H$ = P N HV>P N and the determinant of (H§\k(<p) + b) - k 2 ) (H^ ] N {k{ip) +b)-k' 



We prove all necessary results for HQ' and, then, the analogous results for by 
sending N to infinity. 

Definition 2.13. Let 

&o = \ (u, : o< Pv <4kn ,± Of) , (2.82) 

Obviously, <P\ C <£' C <^o- Note that the circles 0^ deleted from <P to obtain <P' 
are relatively large. Their radius is larger that jk~ 4si ~ s , see (|8.4p . These circles are 
essentially bigger that other 0^ constituting CK 1 ). Properties of are stable in its 
£-4si-2i5 neigh^orhoocl. We denote such neighborhood of <P' Q by <^ . 

Definition 2.14. The set C^ 1 **^) is defined by the formula: 

O^(b) = U 0<Pm<ik ,± Ot(b), (2.83) 

where 0^(6) are discs in complex plane centered at ip^(b), which are zeros of \k(<p) + 
Pm{b)\* — k 2 . The radii of the discs are given by Definition! 



Definition 2.15. The total size of 0W(6) is the sum of the sizes of its connected 
components. 

Lemma 2.16. Let lOOsi < ft < 1/12 - 28si - 145. 

(1) If<p G # \ C (1) (&) and m is such that \\k((p) + p m (b)\l - k 2 \ < k^ 1 , then 

min ||%) +p m+q (b)\l ~ k 2 \ > k?\ (2.84) 

|fcC^)-Kjf5W,C^) 12 — >t 2 1 1 1 C^)l2 — ^ 2 1 1 I^C^) C^)l2 — 1 > l^ 2 ^ (2.85) 

when < p qi ,p q2 < k Sl . This property is preserved in the k~4+2Pi-2si-8 
neighborhood ofS' \ 0^(b). 

(2) The size of each connected component 0^\b) ofO^(b) is less than k~^ , 7 = 
11/6 — ft — 12si —45. Each component contains no more than c\k 2 ^ 3+Sl , c\ = 
ci(di,d 2 ) discs. The total size ofO {l \b) does not exceed 27rci&- 5 / 6+/3l+12si+4<5 . 
The set 0^(b) contains less that cok 2+2si discs. 

Corollary 2.17. For every ip in the k~ 4:+2 P 1 ~ 2si ~ s neighborhood of &' Q \ 0^(b) and 
for every m £ Z 2 

\\k{<p)+p m {b)\ 2 -k 2 \ > k - 2+2 ^- 2s \ (2.86) 
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Proof of Corollary \2.ll\ Obviously, 



If 



\k(tp) +p m+qi2 

\k(ip) +p m (b 



< 



\k(lf) +Prn{b)U ~ k 2 



+ 2 



(k(<p) +p m (b),p qi 



< 1, then, clearly, 



\k(<p) +p m + qi , 2 (b)\ ~ k< 



the corollary follows from (|2.85[) and the last estimate. 



< ck 1+n . Now 



Corollary 2.18. If ip in the fc-4+2/3i-2si-<5 neighborhood of S' \ (b) , then the 
following estimates holds: 

-i 



H^{k{ip)+b)~ l -k 2 



< ck 



2-2/3i+2si 



(2.87) 



The resolvent has no more than cik 2 / 3+Sl poles in each connected component of 
0^'(b). The total number of poles in <Pq is less than cok 2+2si . 

Corollary 2.19. // { c\b) n <P' ± 0, then { c\b) C S' . 

This corollary follows from the statement that the size of Oc (b) does not exceed 

and the definition of <P' as the k~ isi ~ 2& neighborhood of 
The lemma is proven is Section [8T2l Note that the total size fc- 5 / 6 +£i+ 1 2si+4<5 f the 
set OW(b)n#' is small comparing even with the smallest of the circle O q , p q < k Sl . 



Lemma 2.20. If ip G S' \ OW(6), then 



[H(k((p) + b)-k 



2\-l 



<ck 4 



(2i 



The only possible singularities of the resolvent {H (k{ip) + b) — A; 2 ) 1 in each connected 

Oc(b) component of 0^(b), such that 0^\b) fl $q ^ 0, are poles. The number of 
poles (counting multiplicity) inside each component does not exceed c\k 2 ^ +Sl . The 
total number of such poles is less then cok 2+2si . 



Remark 2. In Lemma \2.20\ we could prove convergence of perturbation series and 
the estimate for the resolvent analogous to \2.81^ by somewhat longer considerations. 
However, 12.88\) is good enough for our purposes. The main reason that convergence 
holds, in spite of weaker conditions on ip are formulas of the type (|2.20p . 

Proof. Let <p G S' \ 0^\b). We define the set J%) C I? as follows: 



ft = |m G Z 2 : \ \k((p) +p r< 



\ 2 -k 2 \< 



1 if and only if 



Let Pq be the diagonal projection, corresponding to f2: Po mm 

m G n. Let Pi = J - Po- 
rt follows from Lemma 12.161 that PqW\Pq = 0. Indeed, suppose it is not so. Then, 

there is a pair m,m + q G f2 such that w q ^ 0. Function W\ is a trigonometric 

polynomial, see (|2.3p . Hence < p q < k Sl . This contradicts to (|2.84[) . 
Next, we consider Hilbert identity: 



( ff (i) _ A; 2 ) -1 = - A; 2 )- 1 - (HP - k 2 )-^^ - k 



(i) 
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abbreviations = P"W (k(<p) + bj , Pq 1 ^ = yv((f) + bj being used. Applying 
Pi from two sides and solving for Pi(H^ — k 2 )~ 1 P±, we obtain: 

P^pW -k 2 Y l P 1 = {I + E)- 1 I\{H$ ) -k 2 Y l (l-WiPo(HW -fc 2 ) -1 Pi) , (2.89) 



where £ : Pit 2 -> P^ 2 , £ = Pi(P (1) - k^W^. By the definition of P u 

\\P\{H^p — fc 2 ) _1 || < k~^, (2.90) 

\\£\\ < \\W x \\k-^. (2.91) 
Next, we apply Pq from the left and Pi from the right to the Hilbert identity. Con- 
sidering that PqWxPq = and Pc^P^ - A; 2 ) _1 Pi = 0, we obtain: 

P (P (1) - k 2 y 1 P 1 = -P (H^ - k^WxP^H^ - k 2 y x P x . (2.92) 
Substituting (j!T89"]) into (pT92D gives: 

P (P {1) - k 2 )~ 1 Pi = PP (P« - k 2 )~ 1 P 1 - C, (2.93) 

B = P QPi(I + £)- l P x QP 0y C = P QPi(/ + £ )~ 1 Pi(^ 1) - fc 2 ) -1 , 

Q=(p( 1) -fc 2 )"V 1 . 

It follows from (12T571) that 

IIQH < CiW^k 2 - 2 ^ 2 ^. (2.94) 
Using also (pT90D and (f2T9ll we get: 

\\C\\<C{W l )k 2 - ? ^ +2s \ (2.95) 

Let us prove that 

||P|| < C(VFi)£T 2/3l+8si . (2.96) 

Using (/ + f)- 1 = / + Er=i(- £ Y + (-«f) JV+1 U + ^ = Wi\ we expand P 

into sum oi N + 2 operators: 

2V+1 

B = J2 B r, B = PqQPiQPq, B r = {-l) r P QP 1 (QP 1 ) r QP , r = 1, ...,N, 

r=0 

B N+1 = (-l) Ar+1 P QP 1 (QP 1 ) Ar+1 (/ + £:)- 1 PiQPo. 
Let us estimate ||P r ||, r = 0, ...,N. First, we prove that Bq is diagonal and 

\\B \\ < k- 2l3l+6si . (2.97) 

Indeed, suppose it is not diagonal. Then there is a pair m,m + q £ Q, q ^ 0, such 
that Po mim +g = 0. Considering that W\ is a trigonometric polynomial, we obtain 
that p q < 2Rq < k Sl . Now we easily arrive to contradiction with (|2.84p . Lemma [2. 161 
Thus, 1 1 Po| | = sup mgn \ B 0mm \, 

B = KP 
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Considering that w q = W- q , we easily show that the right-hand side is equal to 

\ 2 -k 2 +p 2 



w q\ ( \H<P) +Pn 



q ^ [\ k (V>) +Pm(b)\ 2 ~ k 2 ) (\k(<p)+p m+q (b)\ 2 - k 2 )(\k{p)+p m _ q {b)\ 2 - k 2 ) 
0<p q < R 

Using ([2^841 and (JXHSJ), we obtain B 0mm = 0(£r 2/3l+6si ). Therefore, (l2T9Tj) holds. 
Similar argument yields that B\ is diagonal and ||-Bi|| = sup mg Q |-Bimm|> 



Bi r 



£ 

qi,q2 € Z 2 , 
< |p, lia | < R , 

qi + qi ^ 



flfcfc) +p m (6)| 2 - fc 2 ) +p m+qi (b)\ 2 - k 2 ){\k(<p) +p m+qi+q2 i 



(2.98) 



Applying (|2.85|) in Lemma 12.161 we obtain 

pill < ||Wi|| 3 ^ 2/3l+8si . 

Similar estimates hold for all B r , r = 2, ...N: 

\\B r \\ < {c\\Wi\\) r+2 k-^- 4si){r+1) . (2.99) 

To estimate Bn+i, we n °te that ||PiQ|| < HWiHA; - ^ 1 . Considering that ||PoQ|| < 
c £,2-2ft+2si an( j i s sufficiently large, we arrive to ||Sjv+i|| < c/c~ 2 ^ 1+8si . 
(f2796|) is proven. Next, considering (pT93j) and ([2795]) . (j2lJ6|) . we obtain: 



Hence, 



k 2 )~ l P l 



Using the last estimate in (|2.89p , we get 

Pl(hV> - k 2 y 1 p 1 



< ck 



2-3/3i+2si 



< ck 



2-4/3i+2si 



(2.100) 



(2.101) 



Taking into account that Pi(ijW — k 2 ) 1 Pq corresponding to a <p is the adjoint of 
Po(H^ — k 2 )~ 1 P\ corresponding to (p and the set 0^\b) is symmetric with respect 
to real axis, we obtain: 



P^H^ 



k'r'Po 



< ck 



2-3/31+2S1 



(2.102) 



Applying P to both parts of Hilbert equations, using Po^i-Po = and (|2.87j) . (|2.102p . 
we obtain: 



Po(H 



(i) 



k^P, 



<ck q 



<ck q 



(2.103) 

Combining (|2~T00l) - (12TT031) . we obtain (12T88D . Note that (I2T88D holds on the bound- 
ary dO^ib) of each connected component Oc(b), such that 0^\b) fl #q 7^ 0. It is 
true, since the size of oi 1 ' '(&) is less than cA;^ 7 , i.e. much smaller than /c -4si-2<5 . 

It remains to show that all singularities of (H^ — k 2 )~ l inside each connected 
component O^c ' '(&), such that 0^\b) n <£' 7^ 0, are poles and the number of poles, 
counting multiplicity, does not exceed c\k 2 / i+Sl . We follow here the approach devel- 
oped in [24] for (— A) 1 + V, I > 6, The plan in [24] is the following. We consider the 



k 2 ) 
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operator I + A := (H^ — k 2 )(H^ — k 2 ) 1 and , use a well-known relation (see e.g. 
[52]): 

det(J + A)-1 < p||iell A " 1+2 , Ae S],. (2.104) 

on the boundary of O^ib). Using Rouche's theorem, we arrive at the conclusion 
that I + A has the same number of zeros and poles inside 0^\b). Therefore, the 
resolvents (H^ — k 2 )~ l and (H^ — k 2 )~ l have the same number of poles. Using 
Corollary 12.181 we obtain that the number of poles of (H^ — k 2 )^ 1 does not exceed 
c\k 2 ^ +Sl . However, there is a technical obstacle on the way of this proof in the 
present case 1 = 1: the determinant of (H^ — k 2 )(H^ — k 2 )~ l is not defined, since 
the resolvent of Hq is not from Si. To overcome this obstacle, we introduce a family 
of expanding diagonal finite dimensional projections P/v : Pn ~^ I- We consider a 
finite dimensional analog of with a multiplier a in front of W\\ 

H^ N = H^ ) P N + aP N W 1 P N , 0<a<l. (2.105) 

First, we show that all (H^N — k 2 )~^, < a < 1, have the same number of poles 
inside each 0^\b). Indeed, let a, ao G [0,1]. We introduce operator A by the formula: 
/ + («_ ao )A = (H^ N -k 2 )(H^ N -k 2 )-\ Obviously, A = P N W X P N {H^ N - k 2 )' 1 . 
Applying formula (12.1041) to I + (a — ao)A, we obtain: 

det(H^-k 2 )(H^ N -k 2 )- l -l\ < \a-a \\\A\\ 1 e^ +2 . (2.106) 

Clearly, \\A\\ < \\Wi\\\\(H^l N - k 2 )- 1 ]]. Considering as in the proof of $2Mb . we 

obtain \\(H^ N (<p) - k 2 )- 1 \\ < ck 4 , when (p is on the boundary of Oc^\b). Since A 
is an N dimensional operator: ||^4||i < ||j4||JV. Combining the last three inequalities, 
we get: \\A\\i < ck A N. Using this estimate in (|2.106p . we obtain: 

det(ff£kfoO - k 2 )(H^ N &) - k 2 )- 1 - l| < c\a - a Q \Nk\ (2.107) 

when tp is on the boundary of { c\b). If | a — o.q\ is sufficiently small, then the right- 
hand side of (|2.107j) is less than 1. By Rouche's theorem det(i/£] v -/c 2 )(i^ 1 o ) Ar -/c 2 )- 1 
has the same number of poles and zeros inside Q { c\b). Therefore, det(^ 1 j v - k 2 ) 
and det(-ff^jy — k 2 ) have the same number of zeros. Covering [0, 1] by small intervals 
and using a finite step induction, we obtain that det(H^ N — k 2 ) has the same number 

^0,N 



of zeros as det(H^ N — k 2 ) for any a G [0, 1]. 

Obviously the zeros of det(H^ N - k 2 ) in Ojpfy stay the same after N surpasses 
a certain number. They are solutions of the equations \k(ip) + p m \ 2 = k 2 in 0^\b). 
Hence, the number of zeros of det(H[ 1] N - k 2 ) is the same for all sufficiently large N . 
It is equal to the number of discs in (b) .We denote this number by M. By Lemma 



SPECTRAL PROPERTIES OF A LIMIT-PERIODIC SCHRODINGER OPERATOR 



27 



[2361 M < Cl k 2 / 3+s K We denote the zeros of det^f^-A; 2 ) as tp\ 



AN) 



(N) 



~<Pm 



ity being taking into account. Obviously, the function n^Li ~ fn 

is holomorphic in O^ib). For a fixed n the sequence {tpn }w— i nas an accumulation 

point. We choose a subsequence iVj such that each {p^ 1 ' 1 converges. With slight 

abuse of notations we drop the index i and consider that each {v J n V ^}^_ 1 
<p n . Considering as in the proof of (|2.88j) . we obtain: 



m(l) 



multiplic- 

k 2 )- 1 



has a limit 



< ck 4 , ip G dO^Xb) 



Therefore, 



M 
n=l 



< cr M fc 4 , 



(2.108) 



(2.109) 



when 99 G dO^(b), r being the size of 0^(6), r < k~ u ^ +l3l+12si+4S . By the 
maximum principle, the above estimate holds inside (b) too. Suppose \ip — <p n \ > 
2e for all n and some e > 0. Then \p — (p nt jsf\ > £ for all n and sufficiently large N. 
Using (|2.109p . we obtain: 

(H^iip) - k 2 )- 1 <cr M k 4 £- M . (2.110) 

Thus, all resolvents are bounded uniformly in N, when <p 7^ ip n , n = 1, ...M. Now, it 
is easy to show now that (H^p N (<p) — k 2 )~ l tends to (Hf\ip) — /c 2 ) -1 in the class of 
bounded operators when <p 7^ <p n , n = 1, ...M. Taking the limit in (|2. 109j) . we obtain: 

M 

(2.111) 



Y[(p-p n )(H^\p)-k 2 )-' 



n=l 



< cr M k\ 



This means that the only possible singularities of of {H^\p) — At 2 ) -1 are poles at 
the points ip = cp n , n = 1, ...M. The number of poles, counting multiplicity does not 
exceed M = c±k 2 / 3+Sl . | 



Let us recall that function x\ (p), defined by Lemma 12.91 is holomorphic in 4>i, 
see Lemma 12.101 Now we prove existence of a function (ip), which is a good 
approximation of o<\{ip) in <Pi and defined on a larger set Let 



u iez 2 ,o< Pi <k 3 /^ - 

The size of the discs constituting O^ 1 ' is at least k~~' '< /4 - 16s i- 9<5 . Let be the 
^-7/4-i6si-io<5_ ne jg] 1 j ;)0r ] 100( j Q £ ^ Since the size of the neighborhood is much 

smaller than the size of the discs, and are the sets of the same type. The 
following relations hold: , <P\ C <P* C <P' C li C l» C S' C @q. 

Lemma 2.21. Let 100si < /?i < 1/12 - 28si - 145. T/iere is a function x*(<p) 
holomorphic in which satisfies the estimates analogous to { 2. 70[ ) on this set: 



k\ < ciw^k-^ 30 ^ 205 , 



dip 



< C(W!)k 



-l+33si+22<5 



(2.112) 



28 Y. KARPESHINA AND Y.-R. LEE 

It satisfies the following estimate when ip £ <Pi : 

\K 1 {<p)-x u {<p)\<ck- l i k ' i/i (2.113) 

Corollary 2.22. Function x* is holomorphic and obeys A2.112\) inside each connected 
component { c\b) ofO^ib), such that oP(b)n$* / 0. 

Proof of the corollary. By Lemma 12.161 (b) is in the cfc _7 -neighborhood of . 
Considering that -7/4 - 16si - 105 > -7, we obtain 0^\b) C 

Proof. Let us consider the function: 

|fc 3 / 4 J+i 

\?\a,x) = x 2 + f u (a,x), x=\*\, f U {a, 2)= (2.114) 

r=2 

Obviously /i*(a!, i?) is a finite sum for the series /i(a, x), see (|2.58p . It satisfies the 
estimates analogous to (|2.59p . when ip £ <P\ 

|/i*(a,x)| < 2a 2 k- 2+30 ^ +20S , \T(m)f u (a, x)| < 2a 2 k -^o S1 +20S+\m\(l+i6s 1+ lU) _ 

(2.115) 

By the construction of <P*, the estimates f|2.4j) and (|2.5p hold for all tp G ^* and 
Pi < A; 3 / 4 . Now it is easy to see that all coefficients #r (x), 2 < r < fc 3 / 4 , are 
holomorphic functions in and even to its 2/c -7//4-16si ~ 10<5 neighborhood. The 
estimates (l2~24l . (|2~26]h (p~4Tjh (p^9]) - (l23ll hold when r < A; 3 / 4 . Therefore, the 
finite sum can be analytically extended to 2/c~ 7 / 4_16,Sl ~ 10<5 neighborhood of <P*, the 
estimates (|2.115j) being valid. Estimating the tail of the series, we get 

\h-fu\ 

Solving the equation A* (a, x) = k 2 for x, we obtain that there is a function x*(c/?) 
defined in the such that (|2,113p holds. This function obviously obeys estimates 
(|2~TT2|) . ■ 

Lemma 2.23. (1) Let lOOsi < ft < 1/12 - 28s x - 145. If ip £ S* \ 0^(b) , and 
m is such that ||x*(</?) + Pm(b)\ 2 — k 2 \ < \k^ x , then 

min \\^(<p) +p m+q {b)\l -k 2 \ > -kP 1 , (2.116) 
0<p q <k s i 1 

\\2*{v)+P m {b)\l-k 2 \ \\x*{y)+p m+qi {b)\l-k 2 \ \\2*{}p)+p m + q2 (b)\l-k 2 \ > Y^ fc2ft 

(2.117) 

when < p qi ,p q2 < k si . This property is preserved in the I/ c - 4 + 2 / 3 i- 2s i-< 5 
neighborhood of&* \ 0^\b). 
(2) Each equation |x*(y?) +p m (^)|* = A; 2 /ias the same number of solutions as the 
"unperturbed" equation \k(ip) + p m (b)\ 2 = k 2 inside every O^ib), Oc \b) H 
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Corollary 2.24. The number of points if satisfying one of the equations | + 

p m (b)\l = k 2 , m G Z 2 , in { c\b), Op '(&) n <P* ^ 0, does not exceed Cl k 2 / 3+Sl . The 
total number of such points is less then cok 2+2si . 

The corollary follows from the above lemma and the last statement of Lemma [2. 161 

Proof. Let tp G \ 0^(b). Noticing that \ £>W(6) C <£g \ we obtain that 

(/? satisfies the conditions of Lemma 12,161 Considering the first inequality in (|2,112p , 
we conclude: 

\\M<P) +Mb)\l ~ +Pm(b)\l\ = O(k- 2+30s ^ 205 ) for any m G Z 2 . (2.118) 

Taking into account that 30si + 205 < 2/3\ — 2si and (|2.86p . we get 



+p m (b)\l ~ \k(<P) +Pm(b)\l\ < \\\k<P)+Prn$)\l ~ k 2 



for any m G Z 2 . 



(2.119) 

Therefore, 2\\x*(ip) + p m (b)\l - k 2 \ > \\Je(tp) + p m (b)\ 2 - k 2 \ for all m G Z 2 . Now, the 
first statement of the lemma easily follows from the first statement of Lemma 12.161 
In particular, (|2.119p holds on the boundary of each 0^\b), { c\b) n 0* ^ 0, 
see Corollary 12.221 By Rouche's theorem, each equation \£*((p) + p m (b)\1 = k 2 has 
the same number of solutions inside every 0^\b) as the "unperturbed" equation 
\k{^)+p m {b)\l = k 2 . U 

Lemma 2.25. Let lOOsi < fa < 1/12 - 28si - 145. //</> G #* \ £> W (&) 3 tten 

(i^(^C^) + 6) - A: 2 ) -1 <cfc 4 . (2.120) 
T/zis estimate is stable in the \ fc-4+2#.-2si-a neighhorhood f ^ \ £>(!)(&). The 
only possible singularities of the resolvent (if (i?* ((/j) + b) — A; 2 ) 1 inside each O^p (b), 
0c (b) Pi ^* ^ 0, are poles. The number of poles, counting multiplicity, inside each 

such Oc{b) does not exceed c\k 2 / 3+Sl . The total number of such poles is less then 
c k 2+2s \ 

Proof. The proof of The lemma is analogous to that of Lemma [2.20l up to replacement 
of Lemma 12.161 by Lemma 12.231 ■ 

Definition 2.26. Let us numerate all components 0^\b), 0^\b)r\<&* 7^ 0, by index 
i, i = 1, /. We denote the poles of (H^ (x* (<p) + b) — A; 2 ) 1 in a component (D^\b)i 
by <fi,n t , n>i = l,...,Mj, Mj < c\k 2 / 3+Sl . Let us consider the discs O s (b)i iTH of the 
radius r^ = k-^Ssi-zs arounc [ these poles. This radius is much less then the size of 
discs constituting {b)i- Let Og be the union of all small discs corresponding 
to a component O^t' 

oW^^U^^^n, (2-121) 
Obviously, oi 1] (b)i C O^fyi and, therefore, { a X \b)i n 0^(6)^ = if i + i' . Let 

0«(6) = uf =1 0( 1 )(6) i . (2.122) 
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Lemma 2.27. The size of each connected component of Os(b) does not exceed 

ciA ,_ 3 I_ 5si _3 5 ^ gize j O Wfy does not exceed Co& -2-4 Sl -35 

Proof. Since each disc has the radius = ]^—^— &s i— 3S an j j^j. Cl fc 2 / 3 + s i ; 

size of a connected component of 0^(6) does not exceed ci&; -3 3 -5si-3<5 . The total 

number of discs in 0^\b) does not exceed cok 2+2si . Now the second statement of 
the lemma easily follows. | 

Further we consider only those connected components Ogc (b) of (b) who have 
non-empty intersection with <&\: 0^(b)i n <Pi 7^ 0. 

Lemma 2.28. Function x\{ip) is analytic in every o£)(b) such that oj$(j>)n$i ^ 0. 
T/ie estimates (I2.112P ZioZc? inside 0^(b). 

Proof. By Lemma [2,27l the length of each connected component of O^ 1 ^ (6) is less than 
ciA; _3 3 _5si_ ' 5 . Thus, each connected component of 0^(1)), which has a non-empty 
intersection with is, in fact, in <&\. By Lemma I2.10f 3). x\{^>) is analytic in in 
and, therefore, in O^pib). | 

Lemma 2.29. If lOOsi < ft < 1/12 - 28s x - 145 and ipeS 1 \ { P{b), then 



(H{)i\(ip) + b) — k 1 



<k j{1 \ J« = 5 Cl k 2 / 3+s \ 



(2.123) 

This estimate is stable in the / c - 4 ^ 6s i~ 4<5 neighborhood of Si \ oi 1 \b). The resolvent 
{H(xi{(p) + b) — k 2 ) 1 is an analytic function of (p in every component of O^' '(&), 
whose intersection with (E>\ is not empty. The only singularities of the resolvent are 
poles. The number of poles in each connected component does not exceed c\k 2 ^ +Sl . 
The total number of poles in all components ofOg(b), whose intersection with <Pi is 
not empty, is less then cok 2+2si . 

Proof. Suppose 99 S Oc Using Lemma [2. 251 we obtain 

Mi 

Y[(^-^ n J(H(^(v) + b)-k 2 ) 



m=l 



< cr Clk2/3 k 4 , 



(2.124) 



r being the maximal size of the components 0^\b) 
{ P{b). Considering that \cp - <p i>rH \ > fc-4-6si-3<5 , 



fc-7. Let if G { c\b)i \ 



{H(x*(ip) + b)-k 2 



we obtain 

< ^,5cifc 2 / 3+s i 



Combining the last estimate with (|2.120p . we obtain that 



(H(2*(<p) + b) -k 2 
Hilbert identity yields: 
(#«(*!(</>) + 6) -k 2 )- 



< k 5nk2/3+Sl when<pe£ 1 \0i l \b). 
[H^{x,{ip)+b)-k 2 Y l 



(2.125) 
(2.126) 



(frW^iC^ + b)-*: 2 )- 1 ^ 

(2.127) 
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£={H {^) + b)-H (>i l (v) + b)){H( 1 \^) + b)-k 2 ) \ 
Let us show that ||£|| < ^. Indeed, we write £ in the form £ = £\A, 

£ x = (H (^(<f) + b) - HoiMv) + b)) {H (M<P) + b) + fe 2 )" 1 , 



A = {Ho&fr) + b) + k 2 ) {H^(Mf) +b)-k 2 y\ 

It is easy to show that £\ is a diagonal operator and ||£i|| < | — >c\{(p)\. Using 
(I2.113|) . we obtain: 

||£i|| < c k- l -^- As ^ kV \ 

It easily follows from formula for A that \\A\\ < 1 + 2k 2 \\((H w (x*(ip) + b ) - k 2 )^ 1 ^ 
Using (|2.125p . we get ||A|| < /c 6cifc2 3+31 . Multiplying the estimates for the norms £\ 
and \A% we get ||£|| = o(l). Now (I2.127P yields: 

H^^ip) + b) - k 2 )' 1 <2 (H^iM^ + V-k 2 )- 1 < fc 5cifc2/3+Sl . (2.128) 

Considering that the discs in 0^(b) have the radius / c - 4 - 6s i- 3<5 ) we easily get that 
the estimate (|2.123|) is stable in the h' 4 - 631 " 15 neighborhood of \ { P{b). 

Next, we use determinants to estimate the number of poles of {H^\ii\{^p) + b ) — 

fc 2 ) 1 . We will follow the scheme established in Lemma 12.201 First, we note that 
(|2.127h and (I2.128P hold not only for i?i(c^), but also for all points in the segment 
between x*(</?) and i.e., for x*((p) + a(xi((/j) — a G [0,1]. Second, 

we consider the family of projectors P/y and the operator Hjp = PnH^Pn. We 
use analogs of (I2.125j) . (|2.127p for operators Hn and a E [0,1]. Using a a multi- 
step procedure similar to that in Lemma 12.201 see (|2.105j> — f|2. 107B . we obtain that 
(Hff(2*(ip) + b) - k 2 y l and (H§\ih(<p) + b) - k 2 )' 1 have the same number of 

poles inside each connected component of 0^\b), whose intersection with <Pi is not 
empty. Considering exactly as in Lemma [2,201 see page 1261 we show that all operators 
f-Hjv (?Z*(<p) + b) — A; 2 ) 1 with sufficiently large N have the same number of poles M, 
M < c\ fc 2 / 3 + s i ; inside a component of of 0^\b), whose intersection with <P± is not 
empty. Hence, the same is true for (H$(ji\{}p) + b ) — k 2 } 1 . Considering further as 
in the proof of Lemma 12.201 we obtain that (H^ {k\ ((f) + b) — A; 2 ) 1 has no more 
that c\k 2 ^ poles. ■ 



2.3.2. The set Og '(6) for small bo. Everything we considered so far is valid for b 
obeying the inequality &o > fc— i—ifisi— 12^ ^ere ^ j s distance from b to the nearest 
vertex of K\. However, in the next section and later, &o w ih be taken smaller, since the 
reciprocal lattice is getting finer with each step. To prepare for this, let us consider 
b being close to a vertex of K\\ 



(2.129) 
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We show that for such b the resolvent (if^ 1 ) (y(<p)) — k 2 ) 1 , y(<p) = %i((p) + b has no 
more than two poles tp^ in <P\. We surround these poles by two contours •yt 1 ^ and 
obtain estimate (|2. 134|) for (y (</?)) — k 2 ) 1 when ip is outside r y^ ± . 

Suppose \b\ = bo, i.e., the closest vertex of K\ for b is (0,0). The perturbation 
series (|2.58p converge for both A^ (i?i (<p)) and A^ (y (</?)) when ip E Si, and both 
functions are holomorphic in Si (Lemma I2.10|) . Note that {oki ((f)) = k 2 for all 
cp E Si. We base our further considerations on these perturbation series expansions. 
For b being close to a vertex e other than (0, 0), we take y(ip) = %i(ip) + b — e. 

We define y>& 6 [0, 27r) by the formula b = &o(cos (pb, sin (/9b) when |6| = 6o, and by 
the analogous formula b — e = 6o( cos y>6j sm Pb) when b is close to a vertex e other 
than (0,0). 

Lemma 2.30. Ifb satisfies h2.129\) and \eq\ < 6o^~ 1_16,Sl_12<5 , then the equation 

A« (jftp)) = k 2 + e (2.130) 
has no more than two solutions, ipf^, in Si. They satisfy the inequality 

\<pf - & b ± n/2)\ < Ia-^-ih (2131) 

Proof. Suppose Wi = and |6| = bo, i.e., the closest vertex of Ki for b is (0,0). 
Then the equation ([2.130P has the form \ku + b\ 2 = k 2 + €q, i> = (cos ip, sin ip). It is 
easy to show that it has two solutions ip^ satisfying (|2.13ip . Applying perturbative 
arguments and Rouche's theorem, we prove the lemma for nonzero Wi- A detailed 
proof can be found in Appendix 1. In the case when b is close to a vertex other than 
(0,0), the considerations are the same up to a parallel shift. | 

Lemma 2.31. Suppose b satisfies \2.129}) . <p E Si and obeys the inequality analogous 
to ( fraij) : \ip - (ip b ± tt/2)\ < jfe-a-Wn-U*. Then, 

^-A« (jfa>)) =^00 ±2b Q k(l + o(l)) . (2.132) 

Proof. Let Wi = and |6| = b Q . Then X^(y((p)) = \kv + b\l and 

d —* — > t^z^ 

— |A;z7+ 61? = 2{ku + b,kp)*, a = — = (— sin </?, cos 
dip dip 

For <£> close to ^ ± 7r/2, we have {b, ft) = ±60 (1 + o(l)). Considering also that 
(jl,v) = 0, we obtain ^|/cz?+6| 2 = ±2&oA;(l+o(l)). Applying perturbative arguments, 
we get a similar formula for nonzero W\. For a detailed proof see Appendix 2. In the 
case when b is close to a vertex other than (0,0), the considerations are the same up 
to a parallel shift. ■ 

Definition 2.32. Let T^ib) be two open disks centered at (p Q E Si |j with the 
radius r (1) = k~ 4 - 6si ~ 3S ; j^(b) be their boundary circles and (b) = U 

rW-. 



V is vr f° r e o = 0. 
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Lemma 2.33. For any ip in Si \ Og (b), 

\\ {1) {y(f)) - k 2 \ > b kr {l \ (2.133) 

This estimate is stable in the k~^~ &Sl ~ iS neighborhood of Si \ 0^(5). 

Proof. Suppose (I2.133P does not hold for some ip in S 1 \ oi 1] (b). This means that 
f satisfies equation (|2.130p with some Eq: < b kr^ . By Lemma T2.301 <p obeys 
(|2.13ip . Thus if could be either ipf or tp~ Without loss of generality, assume 



Obviously, the r^-neighborhood of ipf Q also satisfies conditions of Lemma 



12,311 since is small considering with the size 

k -2-W Sl -llS of the 

neighborhood 



Si. Using (I2.132|) . we obtain that \\ w (y(tp)) - k 2 



2/c6 r (1) (l + o(l)) on the 



boundary of this neighborhood. Applying Rouche's Theorem, we obtain that there 
is a point in this neighborhood. It immediately follows that \(p — <p^\ < r^\ i.e., 

if £ C O s (b), which contradicts the assumption (f € \ O^ib). | 

Lemma 2.34. LetO < b < fc-i-l6«i-i2«5 > For any ip inS^O^ib) or its k~ A - 6si - 45 
neighborhood 



< 



4^.3+6si+4<5 

T ' 



(2.134) 



The resolvent is an analytic function of ip in every component of 0^\b), whose 
intersection with &1 is not empty. The only singularities of the resolvent are poles. 
The number of poles in each connected component does not exceed two. 

Proof. Suppose we have proven that 

-l 



A«(^))-fc 2 )(ffW(^)) 



< 2 



(2.135) 



for all ip in S x . Then, using (I2.133|) with r-W = k~ 4 - 6si ~ 35 , we eagily get f|^.134h . 
Let us prove (|2.135p . By Lemma 12.101 i?i(v?) is a holomorphic function in <F X and 
estimates (p7T0|) hold. Therefore, y(<p) = k(ip) + 0(fc" 1 " 16si " 12<5 ). Using (fXM|) with 
2/3 = 1 — 15si — 95, we obtain 

-l 



< k 



16si+10<5 



z eCi, a £ [0, 1], 



(2.136) 



Ci being given by (|2.8|) . Using considerations similar to those in Lemma I2.2UI (see 
(|2.1U5|) and further), we prove that (^H^ (y (</?)) — z^j has at most the same num- 
ber of poles z inside Ci as (y(ip)} — z^j , i.e., at most one pole. The pole 
obviously exists and is located at the point z = \^(y(ip)). Hence, ( A^^tp)) — 



zj (y(ip)) — zj is a holomorphic function of z inside C x for a fixed ip. From 

the definition of Ci and estimate (I2T29]) it follows that |A (1) (y(c/?)) - z\ < 2fc" 16s ' 1 " 10 ' 5 . 
Multiplying (|2. 136j) and the last estimate, we get 

'\( 1 \y^))-z)(H( 1 \y(v))-zY 1 < 2, z € C x . 
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Using the maximum principle, we obtain the same estimate inside the circle and, 
therefore, for z = k 2 . Thus, we have proved (pX34l) for cp £ $ x \ { ^{b). It is easy 
to see that estimate (|2.134[) is stable with respect to a perturbation of tp of order 

£-4-6si-4<5 

Suppose { P{b) n #1 ^ 0. This means that \b) C S 1 . Considering (12.1351) . we 
see that the only poles of the resolvent are the zeros of the function X^'(y((p)) — k 2 . 
By Lemma 12.301 the number of zeros does not exceed two. | 

3. The Second Approximation 

Let us start with establishing a lower bound for k. Let 7] > 3 • 10 . Since rjsi > 
2 + 4si, there is a number > e such that 

a(l + si)k 2+4si lnk < W 81 , C* = 400(c +ci + l) 2 , c = 32did 2 . (3.1) 

for any k > k*. Assume also, that k* is sufficiently large to ensure validity of all 
estimates in the first step for any k > k*. In particular we assume that all o(l) in 
the first step satisfy the estimate |o(l)| < 10~ 2 when k > k*. 

(2) (2) 

3.1. Operator H a . Choosing S2 = 2s\, we define the second operator H a by the 
formula: 

M 2 

= (1) + aW 2 , (0 < a < 1), W 2 = ^ V r , (3.1) 

r=Mi+l 

where is defined by (|2.ip . M2 is chosen in such a way that 2 M2 ~ k S2 . Ob- 
viously, the periods of W% are 2 M " 2 ~ 1 {d\, 0) and 2 M2_1 (0, c^)- We will write them 
in the form: Ni(a\,0) and iVi(0, 02), where ai,a2 are the periods of W\ and A^i = 
2 a/ 2 -m 1) i k s 2 - Sl <Nl< Note that 

M 2 M 2 

||W2||oo< E H^||oo< exp(-2^)<exp(-F sl ). (3.2) 

r=iWi + l r=iW"i+l 

3.1.1. Multiple Periods of W\(x). The operator = Hq + H / "i(x) has the periods 
ai,a>2- The corresponding family of operators, {H^ (t)}t£Ki, acts in L2(Qi), where 
Qi = [0, a{\ x [0,02] and K\ = [0, 2n/ai) x [0,27r/a 2 ). Eigenvalues of H^\t) are 

denoted by Xn\t), n G N, and its spectrum by A^(i). Now let us consider the same 
Wi(a;) as a periodic function with the periods Nia\, Nia>2- Obviously, the definition of 
the operator does not depend on the way how we define the periods of W±. How- 
ever, the family of operators {H^ 1 ' (i)}teifi does change, when we replace the periods 
a±,a2 by ^101,^102. The family of operators {H^ (t)}t^K 1 has to be replaced by a 
family of operators {H^ 1 \t)} t ^k 2 acting in L2(£?2), where Q2 = [0, Niai] x [0,^102] 
and K 2 = [0, 2ir/N iai ) x [0, 27r/iVio 2 ). We denote eigenvalues of ffW(r) by \£\t), 
n G N and its spectrum by A^ 1 ) (r). The next lemma establishes a connection between 
spectra of operators H^(t) and HW(t). It easily follows from Bloch theory (see e.g. 
[22]). 
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Figure 5. Relation between r (x) and i p 



Lemma 3.1. For any r G -fT 2j 

A (1) M=UA (1) W, (3-3) 

tu/iere 

f = {P = (PuPz) G : < pi < iVi - 1, < p 2 < Ni - 1} (3.4) 
and t p = (t Py i,t Py2 ) = (n + 2irpi/Niai, r 2 + 2-np 2 /Nia 2 ) G i^i, see Fig. £. 

We denned isoenergetic set Si (A) C FTi of by formula (|2.57|) . Obviously, this 
definition is directly associated with the family of operators H^\t) and, therefore, 
with periods ai, a 2 , which we assigned to W\(x). Now, assuming that the periods are 
equal to N\ai, N\a 2 , we give an analogous definition of the isoenergetic set Si(A) in 
K 2 : 

Si (A) := {t £ K 2 : 3n e N : A^(r) = A}. (3.5) 
By Lemma EIH Si (A) can be expressed as follows: 

5i(A) = |r G K 2 : 3n G N, p € P : A^ (r + 2irp/N ia ) = a}, (3.6) 

' \N 1 a 1 , N 1 a 2 J 

The relation between Si (A) and Si (A) can be easily understood from the geometric 
point of view as 

Si (A) =/C 2 Si(A), (3.7) 

where /C 2 is the parallel shift into K 2 , i.e., 

K 2 ■ M 2 -»• if 2 , /C 2 (r + 2vrm/iVia) = r, to G Z 2 , r G FT 2 . (3.8) 

Thus, Si (A) is obtained from Si (A) by cutting Si (A) into pieces of the size K 2 and 
shifting them together in K 2 . 

Definition 3.2. We say that r is a point of self-intersection of Si (A), if there is a 
pair ra, to G N, to ^ rh such that Am (r) = \ } (r) = A. 
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Remark 3.3. By Lemma 13.11 r is a point of self- intersection of Si (A), if there is a 
pair p,j)GP and a pair n,n G N such that \p— p\ + \n — h\ ^ and An'' (r+2irp/Nia) = 
Aj ) (r + 2vrp/iV 1 a) = A. 

Now let us recall that the isoenergetic set Si (A) consists of two parts: Si (A) = 
X|(A) U (Si(A) \ Xi(A)), where Xi(A) is the first non-resonance set given by (I2.73p . 
Obviously /C2Xi(A) C K>2S\(X) = Si (A) and can be described by the formula: 

/C 2 Xi(A) = {r G K 2 : 3p G P : r + 2wp/N ia G tf(A)} . 

Let us consider only those self-intersections of Si which belong to JC 2 Xi(A), i.e., 
we consider the points of intersection of /C2Xi(A) both with itself and with Si (A) \ 

K 2 x!(A). 

Lemma 3.4. A self-intersection r of Si (A) belongs to /C 2 Xi(A) if and only if there 
are a pair p,p G P , p ^ p and a pair n,n G N such that r + 2np/N\a G Xi(A) and 

Xn\r + 2irp/Nia) = X^\t + 2np/Nia) = A, the eigenvalue Xn\r + 2-Kp/Nia) being 
given by the series \2. SJ| ) wrai/i i = r + 2-Kp/Nia and j uniquely defined by t from the 
relation p 2 {t) G E\, 

£l = (k 2 - 3k- Wsi - 115 , k 2 + 3/fc- 16si - 11<5 ). (3.9) 

Proof. Suppose r is a point of self-intersection of Si (A) belonging to /C 2 x*(A). Since, 
r G /C 2 X*(A), there is a p G P such that r + 2-KpjN\a G Xi(A). By Lemma |2~TTI 
there is a single eigenvalue Xn\r + 2irp/N\a) = A of ifW (r + 2-Kp/Nia) in si. It 
is given by the series (|2.22p with t = r + 2np/N\a and j uniquely defined by t from 

the relation p 2 {t) G £1. Uniqueness means: Xn\r + 2-Kp/Nia) 7^ a| (t + 27rp/N\a) 
when h n. Since r is a point of self-intersection of Si (A), A^(r + 2-Kp/N\a) = 
A£ ) (r + 27r]5/iVia) = A for some p ^ p. The converse part of the lemma is trivial. ■ 

To obtain a new non-resonance set X2(A) we remove from IC2X1W a neighborhood 
of its intersections (quasi-intersections) with the whole isoenergetic surface Si (A) 
given by (I3.5p ^ (l3.7p . More precisely, we remove from K,2Xi(X) the following set: 

^i(A) = {rG/C 2 Xi(A):3n,nGN, p,peP, p^p: A«(r + 2-K P /N x a) = X, 

t + 27T P /N ia G Xi(A), |A«(r + 27rp/N ia ) - A«(r + 2irp/N ia )\ < ej, (3.10) 

where H 

e 1 = e-i kVS1 . (3.11) 

We define X2(A) by the formula: 

X2(A)=/C 2 Xi(A) \fii(A). (3.12) 



'Note that £1 is an interval, while e\ is a number. 
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3.2. Perturbation Formulae. Before proving the main result, we formulate Geo- 
metric Lemma: 

Lemma 3.5 (Geometric Lemma). If X > k%, there exists a non-resonance set X2(A, 5) C 
/C2X1 suc h that: 

(1) For any r G X2, the following conditions hold: 

(a) There exists a unique p G P such that r + 2irp/Nia G Xi- 

(b) The following relation holds: 

\f ) (T + 2irp/N 1 a) = k 2 , 

where A^ 1 (r + 2-Kp/N\a) is given by the perturbation series 112. 22}) with 
a = 1, j being uniquely defined by t = r + 2irp/Nia and the relation 

(c) The eigenvalue A^ (r + 2irp/Nia) is a simple eigenvalue of H^(t) and 
its distance from all other eigenvalues \v\t+ 2irp/Nia), h G N of H\{t) 
is greater than e± = e 4 

\xf\r + 2ttp/JVio) - A^(r + 27rp/iVio)| > ei. (3.13) 

(2) For any r in i/ie rea/ (eik^ 1 ^ 6 ) -neighborhood of X2, there exists a unique 
p G -P swc/i i/iai r + 2-Kp/N\a is in the (eifc _1_<5 ) -neighborhood of Xi an d 

\xf\r + 2vrp/A^ia) - k 2 \ < 2eifc- <5 , (3.14) 

j being uniquely defined by t = r ' -\-2irp/N\a j and the relation p 2 (t) G E\. An 



estimate analogous to (|3. 13|) holds: 



2\xf{r + 2irp/N ia ) - \ { }\t + 2irp/N ia )\ > e x . (3.15) 

(3) The second non-resonance set X2 has an asymptotically full measure in x* in 
the following sense: 

L(^2X*i \ X2)) ^ „ u -2-2 Sl 



L{X\) 



<cok' 2 - 2S1 . (3.16) 



Corollary 3.6. If r belongs to the complex (e\k 1 <5 )— neighborhood of the second 
non-resonance set X2(X,5), then for any z G C lying on the circle 

C 2 = {z:\z-k 2 \ = e 1 /2} 1 (3.17) 



the following inequalities hold: 



{hV[t)-z)- x \\ <-, (3.18) 



Corollary is proven in Appendix 3. 



From geometric point of view this means that does not have self-intersections. 
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Figure 6. The set ^3. 



Proof. Let us consider 



e>( 2) = u p / e p\{ }O s 



2tt P ' 



(3.19) 



where O s (-) is defined by (|2.122l) in Definition E2S1 b = Note, that the definition 
make sense, since \ > d^ ax k" S2 > fc- 1 - 16s i- 125 ) d max = max{di, ^2}. Let 

# 2 = # 1 \OW, 6 2 = <£> 2 n [0,2tt). (3.20) 

By Lemma ESS 



< A; J(1) ,J« = Cl A; 2 / 3+Sl . (3.21) 



for all p' e P \ {0} and 99 G #2- Estimate (j3Tj) yields < e x \ We consider 

^l,nonres G ^1- 

T>i,nanrea = {^l(<P), ¥ G ©2} • (3.22) 
We define %2 by the formula: 

X2 = K-2^l,nonres- (3.23) 

By definition of /C 2 , for every r in %2 > there are p <E P and j G Z 2 such that 

2tt P 2irj 



T + 



./Via ^ a 



.1(071 



(3.24) 



Considering (|3.2ip . the estimate k jW < e x 1 and the definition of T>i^ nonres , we obtain: 



ffW r + 



2irp 



<e 1 1 , p = p + p'. 



(3.25) 



Note that the index j does not play a role, since it just produces the shift j of 
indices of the matrix elements of the resolvent. Considering that p' can be any but 
zero, we obtain that (|3.25p holds for all p G P \ {p}. Taking into account that 

Xj (t + \ = k 2 and inequality ([335]) . we arrive at (|3.13p for all p ^ p. It remains 



to check (I3T3D for p = p. Let t = r+J^L. By (IBTm t G K\D X = X i- By Theorem^ 
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Xj yr + is a holomorphic function of r in (eik 1 5 )— neighborhood of X2(A, 

and 

\\f\r + 27T P /N ia ) - \ { }\t + 27vp/N ia )\ > k 2 ^ 2 ' 2 ^ > e x . 

Part 2 follows from stability of all estimates with respect to perturbation of r 
smaller then (eifc -1 " 5 ). Indeed, suppose r is in the real (eifc -1-5 )— neighborhood 
of X2(A, 5). Then, there is a to G X2(A, 5), such that |r — tq\ < eik^ 1 ^ 5 . Let p,j 
be defined by ro as in Part 1. Obviously, perturbation series for A^ (t + 2irp/N\a) 
converges and (I3,14h holds. Using (I3.13h . we easily obtain A3. 15j) . Therefore, p, j are 
defined uniquely by (|3.14p . 

Let us estimate the size of 02- According to Lemma 12.271 the total size of each 

O s (^jj^J is less then cok~ 2 ~ 4si . Considering that the number of p does not exceed 

4k 2si , we obtain that the total size of O2 is less then cok~ 2 ~ 2si . Using (|2.70p . we 
arrive at ()3.16p . ■ 

Remark 3.7. Note that every point (m G Z 2 ) of a dual lattice corresponding 
to the larger periods N\a\,N\a2 can be uniquely represented in the form = 

A^a' wnere m = N\j + p and ^ is a point of a dual lattice for periods a%, a 2 , 
while p G P is responsible for refining the lattice. 

Let us consider a normalized eigenfunction ip n (t, x) of (i) in L2(Qi)- We extend 
it quasiperiodically to Q2, renormalize in L2{Q2) and denote the new function by 
ip n (r,x), t = K,2t. The Fourier representations of ip n (t,x) in L2(Qi) and ip n {T,x) in 
£2(^2) are simply related. If we denote Fourier coefficients of ip n (t,x) with respect 

to the basis of exponential functions |Q7p72 e ^* H a~' x \ j G 1? , in L2(Qi) by C n j, 
then, the Fourier coefficients C nm of ip n (T, x) with respect to the basis of exponential 

C nj , if m = jNx+p; 



functions . |l/2 e "i"' , m G Z , in Li2{Q2) are given by the formula: 



nm 



0, otherwise, 

p being defined from the relation t = r + ■ 2 ^, p £ P. Hence, matrices of the 
projections on ip n (t, x) and t/> n (T, x) with respect to the above bases are simply related: 

(En)mrh, if 771 = jiVl + p, 771 = jiVi + p; 



'33 



0, otherwise, 



i? n and E n being projections in L2{Q2) and L2(Qi), respectively. 

Let us denote by £^ (r + j^) the spectral projection Ej(a, t) (see (|2.23j) ) with 
a = landt = r + |£, "extended" from L 2 (Ql) to L 2 (Q 2 ). 

By analogy with (12.18p . f|2.19p . we define functions #r (&,t) and operator-valued 
functions G> (fc, r), r = 1, 2, • • • , as follows: 

gl 2) (k,r) = ^Tr£ ((#W(r) - *)~V 2 ) r ds, (3.26) 
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Gi 2 Hk,r) = t^ll£ ((^(rJ-^V^PW-*)- 1 ^ (3.27) 

#W(r) being denned at the beginning of Section T3.1.H C 2 being given by (|3.17p . We 

consider the operators H a 2 ^ = H^ + aW 2 and the family H a (r), t € K 2 , acting in 
L 2 {Q 2 ). By (JI2]) and fHI} 

||W 2 ||<ef, (3.28) 
HW2II being the norm of the operator here. 

Theorem 3.8. Suppose r belongs to the (eik~ l ~ 5 ) -neighborhood in K 2 of the second 

non-resonance set Xi{\$)> < S < si, e± = e~± kv 1 . Then, for sufficiently large X, 
A > k 2 and for all a, < a < 1, there exists a unique eigenvalue of the operator 

Ha\r) in the interval e 2 (k) := (k 2 — ei/2, k 2 + ei/2). It is given by the series: 

00 

Al 2) (a, t) = xf ] (r + 27rp/N ia ) + ^ a r gW (k, r), j=j+ p/N u (3.29) 



r=l 

7 2 



converging absolutely in the disk \a\ < 1, where p £ P and j G Z are described as 



( 2 ) 

in Geometric Lemma [375[ The spectral projection corresponding to X-. (a, r) is given 
by the series: 



E?\a,T) = Ef ) (T + 2TTp/N l a)+Y,a r G { ? ) {k,T), (3.30) 

r=l 

which converges in the trace class Si uniformly with respect to a in the disk \a\ < 1. 
The following estimates hold for coefficients gi 2 \k,T), G^\k,r), r > 1: 

\g¥\k,T)\ < |i(4e?r, WG^ik,^ < 6r(4 e ?) r . (3.31) 

Corollary 3.9. The following estimates hold for the perturbed eigenvalue and its 
spectral projection: 



X\ >(a,T) - X) ' (r + 2np/N 1 a) < 12aef, (3.32) 

E? ] (a, r) - Ef ] (r + 2^p/N ia ) ^ < 48ae?. (3.33) 

(2) 

Remark 3.10. The theorem states that X\ (a, r) is a single eigenvalue in the interval 

e 2 (k,5). This means that \X~ 2 \a,r) — k 2 \ <e\/2. Formula (|3.32p provides a stronger 

(2) 

estimate on the location of A~ (a,r). 

3 

Proof. The proof of the theorem is based on expanding the resolvent (Ha\r) — .z) 1 

in a perturbation series for z £ C 2 . Integrating the resolvent yields the formulae for 

(2) 

an eigenvalue of H a and its spectral projection. In fact, it is obvious that 
{H a 2 \ T )-z)- 1 = (H^(T)-z)- l (I-aA 2 )-\ A 2 := -W 2 (H^ (r) - z)' 1 . (3.34) 
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Suppose z E C%. Using Corollary 13,61 and estimate (I3.28|) . we obtain: 



ei 



I A, 



<«< 4e ?<l. 
ei 



(3.35) 



The last inequality makes it possible to expand (/ — a^) -1 in the series in powers 
of aA2- Integrating the series for the resolvent and considering as in the proof of 
Theorem 12.31 we obtain formulae (|3.29p , (|3,30p . Estimates (|3.3ip follow from the 
estimates (|3.35p . ■ 

Next, we show that the series (I3.29p . (|3.30p can be extended as holomorphic func- 
tions of r in a complex neighborhood of \2\ they can be differentiated any number 
of times with respect to r and retain their asymptotic character. 

(2) (2) 

Lemma 3.11. The following estimates hold for the coefficients gi- (k, r) and G)- (k, r) 



in the complex {\e\k 1 d ) -neighborhood of the non-resonance set \2 : 
\T(m)gi 2 Hk,r)\ < ml ■ 3 • 2 2 - 1 +M e ^+ 1 -M fc M(i+^ 
||r(m)G( 2 )(it,r)|| < ml ■ 3r • 2 2 '+ 1 +M < >-M A .M(i+,$) - 



(3.36) 
(3.37) 



Proof. Since ()3.18|) is valid in the complex (e±k' 



-i-S\ 



neighborhood of the second 



(2) (2) 

non-resonance set, it is not hard to see that the coefficients gi- (k,r) and G r (k,T) 
can be continued from the real (ei&; _1_5 )-neighborhood of r to the complex (eiA; -1-5 )- 
neighborhood as holomorphic functions of two variables and inequalities (|3.3ip are 
hereby preserved. Estimating, by means of the Cauchy integral formula, the value 
of the derivative with respect to r in terms of the value of the function itself on the 
boundary of the (^ei/c _1_5 )-neighborhood of r (formulas (|3.3ip ). we obtain (|3.36p 
and (1337) . " ■ 

From this lemma the following theorem easily follows. 

Theorem 3.12. Under the conditions of Theorem ] 3.8\ the series |2 

lesf 

5 ) -neighborhood and the following 



, (TPffJ) can be 

continued as holomorphic functions of two variables from the real (ei/c~ 1_<5 ) -neighborhood 
of the non-resonance set \2 to its complex {e\k~ 
estimates hold in the complex neighborhood: 



\m) (\f\a,T) - \f\r + 2irp/N 1 a) > j 
T{m) (e£\ol,t) - Ef ] {T + 2vrp/iVia) 



<aC m e A - lml k^ l+5 \ 
<aC7 m e?- |m| A ; H(i+5) ) 



(3.38) 
(3.39) 



here and below C„ 
Corollary 3.13. 



48m!2l m l. 



2k 



T{m)\f\a,T) < 2 + 2C(W 1 )k- 1 - 2l3+21s ^ 155 , if 



< 2C(Wi)k- 1 - 2 > 3+15si+115 , k = pj{r + 27rp/N ia ) 



777 



(3.40) 
(3.41) 
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The next lemma will be used in the third step of approximation. The operator 
H^(t) is Ha\r) with a = 1. It will play a role of the initial (unperturbed) operator 
in the third step. 

Lemma 3.14. For any z on the circle C2 and r in the complex (e\k~ l ~ 5 )— neigh- 
borhood of X2, 

|| (i? (2) (r) _^ r i||< A. (3. 42) 
ei 

Proof. Considering the Hilbert relation 

(H^{t) - z)- 1 = (hW(t) - z)- 1 + (#W(r) - z)-\-W 2 ){H^ 2 \t) - z)~\ 
and the estimate f|3. 18|) . together with the estimate (|3.28|) . we obtain: 

ll^ (2) (-) - < i TrnuVVL, * 2|l( ^ (1)(r) " z) ~^ * f • (3 - 43) 

1 — \\(HW(t) — z)- 1 W 2 \\ £ i 

■ 

(2) 

3.3. Non-resonant part of the isoenergetic set of H a • Let S^A) be an isoen- 

(2) (2) 
ergetic set of the operator H a : 52(A) = {r G K2 : 3ra G N : A„ (a, r) = A}, 

f2) (2) 

here {An (a, r)}^L 1 is the spectrum of H a (r). Now we construct a non-resonance 

(2 s ) 

subset X2W °f ^2 (A). It corresponds to non-resonance eigenvalues A^ (r) given 
by the perturbation series (|3.29[) . Recall that T>i(\) nonres and x 2 (A) are defined by 
the formulae (13.221) and (13.23R respectively. Recall also that \2 C IC2X1W (see the 
Geometric Lemma) and xf(A) = K.iD\{X), see (|2.73p . Hence, \2 C £2^1 (A). 

Lemma 3.15. T/ie formula IC2T) i(A) nonres = X2 establishes one-to-one correspon- 
dence between V\{X) nonres and %2- 

Proof. Suppose there is a pair x^xi* G T>x(X) nonres such that JC 2 xi = K. 2 ><u = t, 
r G X2- We introduce also ii = and ii* = /C1X1*. The definition (|2.73|> of 

X*(A) implies that G X*W> since #l>^l* G 2?i (A) nonres C Pi (A). Clearly 

^2*1 = ^2*1* = t and, hence, t\ = r + 2-irpi/Nia, tu = r + 2irp2/Nia for some 
Pi,P2 G P. Now, by Part la of Geometric Lemma [3.51 p\ = P2, and, therefore, 
t\ = ti*. Next, by Lemma T2.121 xi = xi*. | 

We define 62(A) as the set of directions corresponding to the set G2 given by (|3.20|) : 

■52(A) = {veS 1 (A):^Ge 2 }. 

Note that B 2 (A) is a unit circle with holes, centered at the origin, and 62(A) C 61(A). 

Let x G £>i(A) nonres . By (|3.23p . r = /C2X G X2(A). According to Theorem 13. 8\ 

(2) 

for sufficiently large A, there exists an eigenvalue of the operator H a (t), given by 
(|3.29p . It is convenient here to denote A^ (a,r) by \( 2 \a, x). We can do this, since, 
by Lemma f3. lb\ there is one-to-one correspondence between x G T>i(X) nonres and the 
pair (r, j), x = 2nj/a + r. We rewrite (|3.29p in the form: 

00 

A^(a,x) = A«(x) + / 2 (a,x), / 2 (a,x) = ^a^ 2 )(x), (3.44) 

r=l 
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here g). (x) is given by (|3.26p . The function f 2 (a, x) satisfies the estimates: 

|/ 2 (a } x)|<12aef, (3.45) 

|V/ 2 (a,x)| < 96ae?£; 1+5 . (3.46) 

By Theorems l3.8l and l3.121 the formulas (|3.44[) - (|3.46p hold even in the real (eiA; -1-5 )- 
neighborhood of T>i(X) nonres , i.e., for any x = xz7 such that v G Z3 2 (A) and |x — 
xi(A, z7)| < eifc~ 1_<5 . We define T> 2 (\) as a level set for A^ 2 ^(a, i?) in this neighborhood: 

X> 2 (A) := jx = xi?: z7 G H 2 , |x - xi(A, < eiAT 1 ^, A^(a, it) = a} . (3.47) 

Next two lemmas are to prove that T> 2 {\) is a distorted circle with holes. 

Lemma 3.16. For every v G B 2 and every a, < a < 1, t/iere is a unique x = 
x 2 (A, z7) in i/ie interval I 2 := [xi(A, z?) — eifc -1- * 5 , xi(A, z?) + ei&;~ 1 ~ <5 ] suc/i t/iat 

A (2) (a,x 2 ^) = A. (3.48) 

Furthermore, 

|x 2 (A, j?) - xi (A, z?)| < 12aef AT 1 . (3.49) 

The proof is based on (|3.44[) . (|3.45p . (j3.40p and completely analogous to that of 
Lemma 3.41 in [23], set 1 = 1. 
Further, we use the notations: 

x 2 (y>) = x 2 (A,z?), h 2 (<p) = x 2 (ip) - xi (y?), x 2 (^) = x 2 (v?)z7. (3.50) 

Lemma 3.17. Lei 10 -4 < si < 10~ 3 , 77 > 10 4 . 27ien t/ie following statements hold 
for A > k 2 : 

(1) T/ie set V 2 {\) is a distorted circle with holes: it can be described by the 
formula: 

V 2 {\) = {xGM 2 : Z=ic 2 (<p), ^G9 2 (A)}, (3.51) 

where x 2 {ip) = xi((/?) + h 2 (99), xi(y>) is t/ie "radius" of V\{X) and h 2 (ip) 
satisfies the estimates 

dh 2 



\h 2 \ < Uaefk- 1 , 



< 96aefk 1+s . (3.52) 



dip 

(2) T/ie foiaZ length of B 2 {\) satisfies the estimate: 

L (£>i \ B 2 ) < c k~ 2 ~ 2si . (3.53) 

(3) The function h 2 (ip) can be extended as a holomorphic function of <p into the 
complex non-resonance set <P 2 and its £- 4 -6si-4<5 ne ighborhood <^ 2j estimates 
13. 52\) being preserved. 

(4) The curve T) 2 {\) has a length which is asymptotically close to that ofT>\{\) 
in the following sense: 

l(v 2 {\)) a = oo L(p 1 (A))(i + 0(A:- 2 - 2si )), (3.54) 
where 0(Ar 2 - 2si ) = (1 + o(l))c AT 2 ~ 2si ; |o(1)| < 10 -2 when k > k*. 



AA 
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Proof. The proof is completely analogous to that of Lemma 3.42 in |24j . set I = 1. 
Here we give just a short version. Indeed, the first inequality in (13.52H is equivalent 
to (|3.49|) . Differentiating the identity A^ 2 ) (K 2 {f)) = A^ (xi(</?)) = ^ 2 with respect 
to ip and using (13.40p . (13,44p . (13.46p . (13.49p . we easily obtain the second estimate in 
(|3.52p . Estimate (|3.53p valid, since the total size of (2) is less than c £:~ 2 ~ 2si . To 
prove the analyticity of /i2( ( / J ) i n ^2 we check the convergence of the perturbation 
series for A^ 2 ) (xi (</?)), ip G <£ 2 - It is enough to show that 



This inequality immediately follows from 



< — , z G C 2 . 



H 



2np \ 



-1 



< 



(3.55) 



(3.56) 



95. 



proven for all p G P. Let p = 0. By Lemma 12.81 

||(i?W(xi(^)) - z)- 1 !! < 2£T 2 ' 3+1+Sl+5 , whenzGCi, 20 = 1-15*1 

It is not difficult to show that the resolvent has a single pole inside C\ at the point 
z = k 2 . The circle C2 has the same centrum and a smaller radius ei/2. Hence 

4 



(# (1) ("iM) 



< 



ei 



when z G C2. 



(3.57) 



Next we use (|2.123p with b = . The right-hand side of (|2.123p is less that e 1 1 , since 
s\T] > 1. Hence, (|3.56p holds for p ^ too. Convergence of perturbation series follows 
from (j335j) and (pHBJ) . Note that the smallest circle in has the size k~ 4 - 6si ~ 35 . 
Since /(;- 4 -6si-4<5 j g j-^^ smaller than the radius of circles, all estimates are stable in 
the / c -4-6s 1 -4<5 of ^ NoW) 

using Rouche's and Implicit function theorems we easily 
show the equation X^ 2 \xu) = k 2 has a solution x = ^{f) which is holomorphic in 
^2 and coincides with x\(}p) + h%{^p) for real cp. Estimate (|3.54p follows from (|3.52p 
and (I333D . 

Let us record a remark for the sequel. Convergence of the series for the resol- 
vent (if( 2 ) (j?i ((f)) — z) , z G Cii means that the resolvent has a single pole z = 
\( 2 \xi{<p)v) inside Ci- Similar result holds when we replace ><\{<p) by ^(f), since 
#i(y) and ^(y?) are close: \>(2((p)—>(i((p)\ = o(e\). Considering that A^ 2 ) (i?2 (<£>)) = \ 

we obtain that (z — A) (ff^ 2 ^ (x 2 (</?)) — z) 1 is holomorphic inside C2 and the estimate 
following holds: 

" < 32. (3.58) 



(z-A) (h®(M<p)) 



Now define the non-resonance set X2W i n S2W by 

X* 2 (\) := /C 2 P 2 (A). (3.59) 

Lemma 3.18. The set X2W belongs to the (l2ae^/c~ 1 ) -neighborhood of X2W * n 

(2) 

K2- If t G xl(^)> ^ en ^ e operator H a (r) /ias a simple eigenvalue equal to A. T/us 
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eigenvalue is given by the perturbation series 113. 29\) . where p € P,j € Z are uniquely 
defined by r as it is described in Geometric Lemma \3.hA, part 2. 

Proof. By Lemma 13,171 T>2(\) is in the (\2aefk~ 1 ) neighborhood of T>i(X) nonres . 
Considering that X2W = (A) and Xa(A) = K.2'Di(X) nonres (see (I3.23P ). we im- 
mediately obtain that X2W ^ s m 

the (12aef )-neighborhood of X2(A). The size 
of this neighborhood is less than ei/c -1 ^, hence Theorem 13.81 holds in it, i.e., for 
any r £ X2IA) there is a single eigenvalue of H a '{t) in the interval E2(k,5). Since 
X2W C 5*2 (A), this eigenvalue is equal to A. By the theorem, the eigenvalue is 
given by the series (|3.29p , where p £ P,j G Z 2 are are uniquely defined by r as it is 
described in Geometric Lemma |3.5| part 2. ■ 

Lemma 3.19. Formula \3. 59\) establishes one-to-one correspondence between X2W 
and T>2(\). 

Remark 3.20. From geometric point of view this means that X2(A) does not have 
self- intersections. 

Proof. Suppose there is a pair x, x* £ T>2 (A) such that IC2X = IC2X* = r, r G xti^)- 

By the definition flSg7]) of V 2 (X), we have X i - 2 \a 1 x) = X^(a,>i*) = A, i.e., the 

(2) 

eigenvalue A of H a (r) is not simple. This contradicts to the previous lemma. ■ 

3.4. Preparation for the Next Approximation. Let b^ G K2 and b^ be the 

distance of the point b^ to the nearest corner of K2: 

b ( n ] = min \b - 2irm/Nia\. (3.60) 

m=(0,0),(0,l),(l,0),(l,l) 

We assume 6q ^ — l^*- 2 ^. In the case when b is closer to a vertex other than (0,0), 

(2) 

the considerations are the same up to a parallel shift. We consider two cases: 6q > 
eifc- 1 - 25 and < b^ < eiAT 1 " 25 . Let 



y 



V(<p) = x 1 (<p) + bP\ (3.61) 



3.4.1. The case b^ > eiAT 1 " 25 . 

Definition 3.21. We define the set 0^ 2 \b^) by the formula: 

0( 2 )(&( 2 )) = U p6 pO( 1 )(2vrp/iV 1 a + b®). (3.62) 
In the above formula we assume 

Lemma 3.22. If<p€& 1 \ 0( 2 )($ 2 )) or its fc-4-6si-4<$ ne i g hborhood, then 

/ \ — 1 l,4+6si+5(5 

|(^)(^)M)-fc 2 ) ||<^— (3-63) 

T/ie resolvent is an analytic function of (p in every connected component (b^ ) 
of C( 2 )(fe( 2 )) J whose intersection with <P\ is not empty. The only singularities of the 
resolvent in such a component are poles. The number of poles (counting multiplicity) 

of the resolvent inside Oc (b^) is less that cok 1+2s2 . The size of each connected 
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component is less than k 3 4si 25 . The total number of poles (counting multiplicity) 
of the resolvent inside 0^(b^) is less than cok 2+2s2 . The total size of 02(6^) is 
less then cok~ 2 ~ 2si . 

Proof. Combining Lemmas 12.291 and 12,341 for = 2Trp/N±a + 

6(2) 

, we obtain: 



4£,3+6si+45 



+2Trp/N 1 a + b (2) ) - k 2 j <maxj£; J , j, (3.64) 

J« =5 Cl k 2 / 3 ^ <c k 2+2s \ p£P, 
for any ip G Si \ 02 (fr*- 2 **)- Using (|3. 1 j) . we easily obtain k jW < e^ 1 . Considering the 
condition b^ > eiAT 1 ^ and the definition of , we obtain (ET63]) in S 1 \ C> (2) (fr (2) ) 
The estimate is stable with respect of the perturbation of order / c - 4 - 6s i- 4<5 ) since the 
size of the discs in e> 2 (£ (2) ) is Ar 4 ~ 6si - 35 '. 

Now we prove the second part of the lemma. Let A* be a rectangle in C: A* = 
{(p : |3ffy? — < k^ 1 , < A;~ 5 } for some y?* G [0, 27r). It is shown in the proof of 
Lemma 18.101 (see the estimate for the number of points in I*) that the total number 
of points V'm(^) i n ^* does not exceed 5co& 1+2si for any b G K\. Therefore the 
number of discs of 0^\b) which intersect A* does not exceed 5cofc 1+2si . Therefore, 
the number of discs of Os(p) which intersect A* does not exceed 5cok l+2si for any 
b G K\. Hence, the number of discs of 0^ 2 \b^) which intersect A* does not exceed 
5cok 1+2s2 . Considering that the size of each disc is less then / c - 4 - 6s i- 35 ; we obtain 
that the total size of {2 Xb {2 ">) n A* is less than ^.-3-4 Sl -2<5 _ It ig obvious 

now that 

the size of each connected component of C*- 2 ^(6^ 2 ^) does not exceed /c _3_4si ~ 25 . If 
Oc 2) (& (2) ) n / 0, than of\b {2) ) C S 1 . Therefore, the resolvent is an analytic 
function in this O^ib^). By construction, the poles of the resolvent are in the 
centra of the discs in 0( 2 )(6( 2 )). Therefore, the number of poles in each connected 
component does not exceed 6coA; 1+2s2 . 

By Lemmas 12.291 and 12.341 the number of poles of the resolvent (^H^ ^Ki{ip) + 

2irp/N 1 a + bW) -k 2 ^j * inside O s (2iTp/N 1 a + b^) is less than cok 2+2si for any p G P. 

Considering that there are no poles in outside O s (2irp/Nia + b^), we obtain that 

the number of poles of the resolvent (^H^ (xi(c^) + 2irp/N\a + b( 2 ^ — k 2 ^j inside 

0( 2 )(£>( 2 )) is less than cok 2+2si . Taking into account that the number of p G P does 
not exceed ck 2 ^ S2 ~ Sl \ S2 = 2s\, we obtain the estimate for the number of poles for 

the resolvent LffM^ 1 ) ((</?)) - A; 2 ) . 

We estimate the size of 0( 2 )(6( 2 )) the same way as we estimated the size of 0^ 2 \ 

Indeed, according to Lemma 12.271 the total size of each O s + U 2 ^ is less then 

cofc~ 2-4si . Considering that the number of p does not exceed Ak 2si , we obtain that 
the total size of 2 (6( 2 - ) ) is less then cok~ 2 ~ 2si . 
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Definition 3.23. We denote the poles of the resolvent (yW(^)) — k 2 ^j in 

OW(bW) as p { n\ n = 1, M^ 2 \ < Co k 2+2s2 . Let us consider the discs {2) (b^) 
of the radius = ^-2- 4s 2-5 r (i) = £,-6-7s 2 -4<5 cen t erec j a t these poles. Let 

Oi 2 \b^) = U^Oi 2 \b^). (3.65) 

Lemma 3.24. The total size of of ] {b^>) is less then c k- 2s2 - s r^ = c /c~ 4 ~ 5s2 ~ 4<5 . 

Proof. The lemma easily follows from the formula = fc -2 — 4s 2-<* r (i) _ £,-6-7s 2 -4<5 
and the estimate M (2) < c k 2+2s2 . | 

Lemma 3.25. If (p E &i \ of\b^), then 



< 



(3.66) 



T/ie estimate is stable in the {r^k ^-neighborhood of 0\\oi 2 \b^). The resolvent 

is an analytic function of ip in every component of oi 2 \b^), whose intersection with 
is not empty. The only singularities of the resolvent are poles. The number of 
poles (counting multiplicity) of the resolvent inside oi 2 \b^) is less than ck 2+2s2 . 

Proof. By the definition of oi 2 \b^), the number of poles (counting multiplicity) of 
the resolvent inside this set is less than ck 2+2s2 . Considering as in Lemma 12.291 we 
obtain 



< v 



-M( 2 ) 



k 4+5 



El 



(3.67) 



where v is the coefficient of contraction, when we reduce (2 \b (2 ^ to of\b^). 
Namely v is the ratio of A 2 ^ to the maximal size of OW(bW). By Lemma E22] 
and the definition of A 2 \ u = k~ 3 ~ 5s2 ~ 2S . It is not difficult to show that that 



-Af< 2 > 



k i+s < e 1 x , when k > k*, see (|3.ip . The estimate (|3.66p easily follows. 



Obviously the total size of oi 2 \b^) is much smaller the smallest disc in 0^ 2 \ 
Therefore, the function ^(v 9 ) is holomorphic inside each connected component of 
{2 \b^) which has a non-empty intersection with (f> 2 - Let 



yt 2 \p) = Z 2 (p) + b^ . 
Lemma 3.26. If<pe£ 2 \ {2 \b^), then 

tf( 2 )(^V))-fc 2 



2 



(3.68) 



(3.69) 



The estimate is stable in the [r^k 5 ) -neighborhood of <P 2 \ oi 2 \b^). The resolvent 

is an analytic function of cp in every component of 0^\b^), whose intersection with 
<P 2 is not empty. The only singularities of the resolvent are poles. The number of 
poles (counting multiplicity) of the resolvent inside oi 2 \b^) is less than ck 2+2s2 . 
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Proof. We use the Hilbert identity 



where y (2) are given by (|3lUjl . (^68]l . £ = £ x + £ 2 £ 3 , 

fl = - fc 2 )~V 2 , £ 2 = (#«(^)) - k^)\ii^{f-)) + fc 2 

Using the estimates (13.28P and (|3.66|) . we see that ||£i|| < ef. Using (I3.66p . we get 
||£ 2 || < 1 + || Wi - 2k 2 \\e^ 2 < 3k 2 e^ 2 . It is easy to show that ||£ 3 || < k~ x \x\ - x 2 |. 
Considering (13.491) . we get ||£ 3 || < 12efk~ 2 . Therefore, ||£|| < 4ef < 1/2 when k > k*. 
It is easy to see now that 

2 



H^U 2 ))-k 2 







r 


< 2 



HW(^)-k' 



i 



Thus, we have obtained (|3.69p . Introducing operators Pat and using the same tech- 
nique as in the proof of Lemma 12.201 (see (|2.105p and further) we prove that the 

number of poles of ( (y^ 2 ^^) — k 2 ) obeys the same estimate as the number of 



poles of (jl^(y^ (</?)) — k 2 ^ , i.e., it does not exceed cok 2+2s2 . | 

3.4.2. The set of ] (fo (2) ) for small b^. Everything we considered so far is valid if 

6q 2 ^ > eik^ 1 ^ 25 . However, in the next section and later, is taken smaller, since the 
reciprocal lattice is getting finer with each step. To prepare for this, let us consider 
&( 2 ) being close to a vertex of K2' < b^ < e\k~ l ~ 2& . We show that for such 

&( 2 ), the resolvent (^H^(y< 2 \ip)) - k 2 ^j has no more than two poles in ^ 2 . We 

surround these poles by two contours 7 =t ^ and prove an estimate for the norm of 

^£f( 2 ) (y( 2 )((^)) — k 2 ^j when cp is outside these contours. 

Suppose \b^\ = 6q 2 \ i.e., the closest vertex of K2 for b^ is (0,0). We consider 
the functions A^^ 1 ) ((/?)) and A^^y*- 2 ) ((£>)) defined by perturbation series (|2.58p and 
(|3.44p for ip G # 2 - The convergence of these series can be easily justified. In fact, 
by Lemmas 12.101 and 13.171 >?i(y) and i? 2 ((/?) are holomorphic functions of ip in <£ 2 . 
The perturbation series (I2.58P and (I3.44|) converge for A( 1 )(i?i(^)) and A( 2 )(x 2 (^)), 
respectively. Since the estimates involved are stable with respect to a change of h\ 2 
not exceeding e\k~ l ~ & , the perturbation series for A^ (y^ ((</?)) and A^ 2 ^ (jf 2 " 1 (cp)) also 
converge, both functions being holomorphic in ^ 2 . We base our further considerations 
on these perturbation series expansions. For b^ being close to a vertex e other than 
(0, 0), we take y@\(p) = Z 2 (<p) + b^ - e. 

From now on, we denote the solutions tpf Q of the equation A*- 1 ) (y^ 1 ) (((£>)) = k 2 + Eq, 
introduced in Lemma EMU by ip^. We set r^ = r (l)fc-2-4s 2 -5_ 
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Lemma 3.27. If < 6 Q 2) < e x k 



,-1-25 



and |eo| < b^k 1 s r^ 2 \ then the equation 



A (2) (^ 2 V))=A: 2 + 



eo 



(3.70) 



has no more than two solutions ip^^ in &2- For any y>eo ^ ^2 there is ^ 
such that 



(1)± 



< r (2)/ 4 



(3.71) 



here and below (f^^ is ip^^ for eo = 0. 

Proof. First, we expand \W{y< 2 \ip)) - A^HyW^)) near the point b^ = and 
consider that A( 2 )(x 2 (^)) = AW(xi (</?)) = A. Then, using (pM) . (l27o3l) . Q374g]) and 
(j3"~4Tj) . we check that 



A( 2 )( y -< 2 )M)-A«(^)M) <^ ei 



(3.72) 



in #2 and even in its r (^-neighborhood, the neighborhood being a subset of Si. 
Suppose (ETTOD holds for some <p 6 <£ 2 - By ([37721) . A^yW^)) = A: 2 + e' , e < 
eo + &o^ e i < 2&g 2 ^/c 1 ~'M 2 ) when k > k*. Hence, 99 satisfies conditions of Lemmas 
12.301 and 12.311 Surrounding tp by a circle C of the radius r( 2 )/4 and using (|2.132p . 
we see that (if- 1 ^ (ip)) — k 2 — e' \ « |fc& Q 2 M 2 ) » |e | on this circle. Applying 
Rouche's theorem, we obtain that there is a solution of AM(yW (</?)) = k 2 inside this 
circle. Thus, any solution of (|3.70p is in the circle of the radius r^ 2 ^/4 surrounding 

<£q 1 the point being in the (r( 2 )/4)-neighborhood of <f>2- It remains to check 

that (|3.70p has no more than two solutions and no more than one in a vicinity of 

each </?o . We construct the disk of the radius r^ 2 ^/4 centered at <Pq^ and note 
that |AW(yW (^)) - k 2 \ > \kbf\^ outside the circle. Using ([3772]) and Rouche's 
Theorem, we obtain that there is only one solution of (|3.70p in the disk. | 

- 1-L-2S an d (p £ <p 2 obeys the inequality analo- 



Lemma 3.28. Suppose < b < e\k 

< r^ 2 ) . Then 



gous to HTW: 



±(i) 
<P-<Po 







dip' 



A( 2 ) U 2 \ip) 



=^00 ±26^^(1 + 0(1)), (3.73) 

where \o(l)\ < 10~ 2 + e\ when k > k*. 

The proof is completely analogous to that of Lemma 12.311 

Definition 3.29. Let (b&A be the open disks centered at ip^^ with radius 

r (2). 7 (2)± 



(b^) be their boundary circles and of ] (b^) = U T 



,(2)- 



Lemma 3.30. For any ip> in S 2 \ of\b {2 ">), 

\\M$ 2 \<p))-tf\>b®&- B rV>. (3.74) 
The proofs of this and the next lemma are analogous to those of Lemma [2. 331 12.341 
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Lemma 3.31. For any ip£& 2 \ of 1 (& (2) V 



Hm(f>^-^l <w ^-, (,75) 

The estimate is stable in the (r^kr 5 } -neighborhood of <P 2 \ (b^ ) . The resolvent 

is an analytic function of ip in every component of 0$ (b^), whose intersection with 
^>2 is not empty. The only singularities of the resolvent are poles. The resolvent has 
at most two poles inside 0^\b^). 

4. The Third Approximation 



Squaring both sides of (|3.ip . we easily obtain the relation which coincide with (|3.ip 
up to the substitution of s 2 instead of s±: 

C*{l + s 2 )k 2+As nnk < kP s \ C* = 400(c + l) 2 , c = 32d 1 d 2 , s 2 = 2 Sl . (4.1) 

for any k > fc*. We will use (|4.ip in the next step. 

(3) • (3) 

4.1. Operator H a • Choosing S3 = 2s 2 , we define the third operator H a by the 
formula: 

M 3 

= H® + aW 3 , (0 < a < 1), W 3 = ^ V r , 

r=M 2 +l 

where M3 is chosen in such a way that 2 Mi ~ A; S3 . Obviously, the periods of W3 
are 2 Af3-1 (di,0) and 2 M3_1 (0, d 2 ). We write them in the form: N 2 Ni(a 1 ,0) and 
N 2 N 1 (0,a 2 ), here N 2 = 2 M3 ~ M *, \k s *~ S2 < N 2 < Ak s ^~ S2 . Note that 

M 3 M 3 

IIW3II00 < £ IIKII00 < £ ex P (-2^) < exp(-fc" S2 ). 

r=M 2 + l r=M 2 +l 

4.2. Multiple Periods of W2(x). The operator, = H\+W 2 (x), has the periods 
N\ai, N\a 2 . The corresponding family of operators, {H^ (t)} T £k 2 , ac ts hi L 2 (Q 2 ), 
where Q 2 = [0,^%] x [0,-/Via 2 ] and K 2 = [0,2^/^^) x [0, 27r/Aqa 2 ). Since now 
on we denote quasimomentum t from the first step by quasimomentum r from 
the second step by t^ . Correspondingly, quasimomentum for H a 3 ^ we denote by 
t(3). Eigenvalues of .ff( 2 ) (i^) are denoted by A^ (i^), n E N and its spectrum by 
A( 2 )(t( 2 )). 

Next, let us consider W 2 (x) as a periodic function with the periods N 2 N\a\, 
N 2 N\a 2 . When changing the periods, the family of operators [t^ 2 ') } t m eK 

is replaced by the family of operators |if( 2 ) (t )} t (3) 6 ^ ) acting in L 2 (Q3)j where 
Q 3 = [0,N 2 N ia i] x [0,Ar 2 Aqa 2 ] and K 3 = [0,27r/N 2 N iai ) x [0, 2ir/N 2 Nia 2 ). We 
denote eigenvalues of (t^) by An (i^), n G N, and its spectrum by A^ 2 ) (t®). 
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We denote now by P^ 1 ' the set P, introduced by (|3.4|) . its elements being p^\ By 
Bloch theory (see e.g.|22j), for any i {3) G K 3 , 



A (2) ^(3) j = |J A (2) ^(3) + 27rp (2) /A r 2iVl aj , (4.2) 

p (2) g p(2) 

where P( 2 ) = jj/ 2 ) = (pf } ,p^ 2) ) G Z 2 : < pf ) < N 2 - 1, < < iV 2 - l}, 

\N 2 Niai N 2 N±a 2 J 
An isoenergetic set ^(Ao) C K% of the operator is defined by the formula: 
5 2 (A) = p>^ 3 :3nGN: A^ 3 ))=a} 

= [t® G K 3 : 3n G N, p( 2 ) G P^ : (V 3 > + 2irp^ /N 2 N ia ) = a} . 
Obviously, S 2 = K.^S 2 , where /C3 is the parallel shift into K3, that is, 

/C 3 : R 2 -> K3, /C 3 (V 3) + 2irm/N 2 N 1 a) = i (3) , m G Z 2 , i (3) G if 3 . 



We denote index j, introduced in Part 1 of Lemma 12.11 (Geometric Lemma for the 
First approximation), by and j, introduced in Part 1 of Lemma 13.51 (Geometric 
Lemma for the Second approximation), by j( 2 \ = + p^ 1 ' /N\. 

4.3. Perturbation Formulae. 

Lemma 4.1 (Geometric Lemma). For a sufficiently large X, A > k 2 , there exists a 
non-resonance set X3(A, 8) C /C3X2 suc ^ that: 

(1) For any point t^ G X3> the following conditions hold: 

(a) There exists a unique p^ G P^ such that t^ + 2-Kp^ /N 2 Nia G x 2 - 

(b) The following relation holds: 



X% (t^ + 2«pW/N 2 N ia )=\, 



where A^j (i® + 2irp^ /N 2 N\a) is given by the perturbation series {2 



with a = 1 and = j + p/N\, here j and p are defined by the point 
t = t^ + 27rp( 2 ) /N 2 N\a as it is described in Part 1(b) of the Geometric 
Lemma for the previous step. 
(c) The eigenvalue (t^+2'KpW /N 2 Nia) is a simple eigenvalue of (t^) 

and its distance to all other eigenvalues An (t^ + 2irp^ /N 2 N\aj of 
H@)(tW) is greater than e 2 = e~^ 



\fl (V 3 ) + 2^/N 2 N ia ) - Af (V 3 ) + 27rpW/N 2 N ia ) |> e 2 . (4.3) 
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(2) For any in the real (e 2 k 1 5 ) -neighborhood of \3> there exists a unique 
p (2) g p(2) suc /j ^ t (3) + 27rp (2) /JV 2 iVl a is in £/je (e 2 k 1 5 ) -neighborhood of 
X 2 an d 



X% (tW + 2^/N 2 N ia )-k< 



< e 2 k~ 



(4.4) 



J 



(2) 



J 



+p/N\, here j andp are defined by the point r = t^ + 2-7rp( 2 ) /N 2 N%a 
as ii is described in Part 2 of the Geometric Lemma for the previous step. An 
estimate analogous to (|4.3[) holds: 



\\f(t<® + 2irpW/N 2 N ia ) - A£V 3) + 2TTp^/N 2 N ia )\ > (l + o(l))e 2 . (4.5) 

(3) The third nonresonance set X3 has an asymptotically full measure on x 2 in 
the following sense: 

L(IC 3X * 2 \X3)) 



< cok 



-4-2si-2s 2 



(4.6) 



L(X* 2 ) 

Proof. The proof of the lemma is analogous to that for Geometric Lemma in the 
second step. Indeed, let us consider 

/ 2vrp( 2 ) \ 



o (3) = Up (2)eP(2)U0} o( 2 ) 



\N 2 N ia J 



(4.7) 



when- of ] ( is defined by (^65]) with b^ = -g^. Note, that the definition 

> d-l x k- s 3 > e 2 k- l ~ 25 . Let 



\N 2 Nu 
make sense, since | NiNia 1 



<£ 3 = # 2 \0( 3 ), G 3 = <^ 3 n [0,2vr). 



By Lemma 13.261 



* <2) (^>+f£)- 



<3 



(4.8) 



(4.9) 



for all p^ G P^ \ {0} and 92 G ^3, here and below ^3 is the (r^k s ) -neighborhood 
of <P 3 . We consider T> 2tnonres C V 2 : 

^2,nonres = {>?2(V)) V ^ ©3} . (4.10) 

We define X3 by the formula: 

X3 = ^3^2,nonres- (4.11) 

By the definition of /C3, for every in X3, there are G p (2) G P (2) and 

j G Z 2 such that 

2vrp (2) 2vrp (1) 2vrj 



t( 3 ) + 



+ 



+ 



N 2 N%a Nia a 
Considering the definition of T> 2 nonres , we obtain: 



K 2 (ip), >c 2 ((p) G T> 2 , 



H (2)( t (3) 



+ 



2irpM 
N 2 N 1 a 



k 



< 



-2 ' 



(4.12) 



(4.13) 
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p( 2 ) = pi 2 ) _|_p( 2 )_ Note that the indices j, p^ do not play a role, since they just 
produce a shift of indices of the matrix elements of the resolvent. Considering that 
p( 2 ) can be any but zero, we obtain that holds for all p^ G P (2) \{p {2) }. Taking 

into account that A^ 2) (^ 3 ) + 2irp^ /N^a) = k 2 and e 2 < e\ , we arrive at (14. 5p for all 
p^ p. It remains to check (g3]) for p = p. Let = + 2irp^ /N 2 N 1 a. By (|4"7l2j) . 
t^ 2 ) G fC 2 V 2 . Using (I339|) . we get *( 2 ) G x|. By Theorem ESI Aj 2) (t( 3 )+27rp( 2 )/iV 2 Af 1 a) 
is the only eigenvalue of ff( 2 ) (t^ + 2-Kp^ /N 2 N\a) in the interval £2- Hence, 

lAfV 3 ) + 2irpW/N 2 N ia ) - 4 2) (* (3) + 2vrp( 2 )/iV 2 Af 1 a)| > £l . 

Part 2 holds, since all estimates are stable with respect to the perturbation of 
less then e 2 k~ 1 ~ s . 

Let us estimate the size of . According to Lemma 13.241 the total size of each 
of ) {2-K P ^/N 2 N 1 a) is less then c k 4 5s2 46 . Considering that the number of p( 2 > 
does not exceed fc 2 ( S3-S2 ) = k 2s2 , s 2 = 2s%, we arrive at (14. 6p . 

■ 

Corollary 4.2. If t@' belongs to the complex {e 2 k~ l ~ 5 ) — neighborhood o/x3(A, 5), 
then for any z lying on the circle C3 = {z : \z — k 2 \ = e 2 /2}, the following inequality 
holds: 

(HW(tW)-z)- 1 <-. (4.14) 
(-2 

Remark 4.3. Every point 2irq/N 2 Nia (q G 1?) of the dual lattice for periods 
N 2 N\a\, N 2 N\a 2 can be uniquely represented in the form: 2nq/N 2 Nia = 2irm/Nia+ 
2irp^ /N 2 Nia, where m G Z 2 , p^ G P^ . Note that 27rm/Nia is a point of a dual 
lattice for periods Nia%, N\a 2 and p^ G P^ is responsible for refining the lat- 
tice. By Remark 13.71 2-Kq/N 2 N\a also can be uniquely represented as 2-Kq/N 2 N\a = 
2irj /a + 2irp^ /N ia + 2irpW /N 2 N ia , here j G Z 2 , pW G Pto, p( 2 ) G P< 2 ). 

Let us consider a normalized eigenfunction if) n (t^ 2 \x^ of i?( 2 ) (^ 2 **) in L 2 (Q 2 ). 
We extend it quasiperiodically to L 2 (Qs), renormalize and denote the new function 
by VVi (t( 3 \%), = ICst^ 2 \ The Fourier representations of ip n (t^ 2 \x) in L 2 (Q 2 ) 
and tpnit^ ,x) in L 2 (Qs) are simply related. If we denote Fourier coefficients of 
V>n (t (2) ,x) with respect to the basis \Q 2 \-i/ 2 e < 2 ™/Nxa+tV) ,*) > m £ ^ in L2 (q 2) by 

C nm , then, obviously, the Fourier coefficients C nq of ij) n (p 3 ',i) with respect to the 
basis |Qg|-i/V(^/* a * iaft(3) ' x ), q £ ^ in £ 2 (q 3 ) are given by the formula: 

if q = mN 2 + p^ ; 
otherwise, 

p( 2 ) being defined from the relation t^ = t^ + 2-np^ /N\N 2 a. Correspondingly, 
matrices of the projections on ip n (T,x) and ^p n (t^ 3 \ x) with respect to the above 
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bases are simply related: 



{E n )qq 



(En) mm, ^q = mN 2 + p( 2 \ q = mN 2 + p^; 
0, otherwise, 



E n and E n being projections in L 2 (Q 3 ) and L 2 (Q 2 ), respectively. 

We define functions gi 3 \k, t^) and operator- valued functions G^ (k,t^), r 
1, 2, • • • , as follows: 



r+l 



o 3 

( t (3)) _ ^"V 3 V (H^(t^) - zY'dz. 



HW(t®)-z) W 3 ) dz, 



Theorem 4.4. Suppose belongs to the (e 2 k~ 1 ~ ~* ') -neighborhood in K% of the third 
nonresonance set X3(^i$)> < 5 < s\, e 2 = e~^ kn 2 . Then, A > k 2 and for all a, 
< a < 1, there exists a unique eigenvalue of the operator Ha ' 1 (i*- 3 * 1 ) in the interval 
e^(k) := (k 2 — e 2 /2, k 2 + e 2 /2). It is given by the series: 

00 

\%(a,tV>)=\% (tW + 2*pW/N 2 N ia ) + £c/#(M (3) ), (4.15) 

r=l 

converging absolutely in the disk \a\ < 1, where := + /N%N\, p^ 2 \ j^' 
being described in Geometric Lemma \4-l\ The spectral projection, corresponding to 
A^(3) (a, t^), is given by the series: 

00 

£$)M 3) ) = E% (tW + 2npW/N 2 N ia ) +Y,crGW(k,tW), (4.16) 

r=l 

which converges in the trace class Si uniformly with respect to a in the disk \a\ < 1. 
The following estimates hold for coefficients gfj?\k,t^), Gl?\k,t^): 

»W(M (8) )| < N)' > IK^M^IL < 6r-(4el) r • (4.17) 

Corollary 4.5. The following estimates hold for the perturbed eigenvalue and its 
spectral projection: 



\%(a,t^)-\%(t^ + 2npV/N 2 N 1 a) 



E 



%{a^)-E%{t^ + 2^/N 2 N l a) 



£-(2) 



< I2aet 



< 48ae 3 ,. 



(4.18) 
(4.19) 



Proof of the theorem is analogous to that of the Theorem l3.81 The series (|4.15p . (|4.16|) 
can be extended as holomorphic functions of t^ in the complex (e2fc _1_<5 )-neighborhood 
of X3; they can be differentiated any number of times with respect to and retain 
their asymptotic character. The results analogous to Lemma 13.11} Theorem I3.12[ 
Corollary 13.131 and Lemma 13.141 hold. 
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4.4. Nonresonant part of the isoenergetic set of Ha \ This section is analogous 
to Section [3731 for the second step. Indeed, let £3 (A) be an isoenergetic set of the op- 
erator Hi 3) : S 3 (A) = {t& G K 3 : 3n G N : Al 3) (a,i( 3 )) = A}, here {Al 3) (cM (3) )}~=i is 
the spectrum of Ha\t^). Now we construct a non-resonance subset X3W °f ^(A). 
It corresponds to non-resonance eig envalues A ( | } (^ 3 )) given by the perturbation se- 
ries (|4.15p . Recall that T>2(X) nonres and %3 are defined by formulas (14,10p and (14. lip . 
Recall also that X3 C IC3X2W ( see Geometric Lemma [47T]) and X2W = ^2^2 (A), see 
(|3.59p . Hence, X3 C £3^2 (A). The following lemma is analogous to Lemma 13.151 in 
the second step. 

Lemma 4.6. The formula fC3T>2(X) nonres = X3 establishes one-to-one correspondence 
between T> 2 {\) nonres andx3- 

Proof. The proof is analogous to that of Lemma 13.151 up to the shift of indices by 
1, i.e., X2 — > X3; X*W ~^ X2W) T = ^ ~~ ^ * > we use formula (13.59P instead of 
(|2.73p . Part la of the Geometric Lemma for the third step instead of Part la of the 
Geometric Lemma for the second step, and Lemma 13.191 instead of 12.121 | 

We define #3 (A) as the set of directions corresponding to 03, 63 being given by 

S 3 (A) = {ve6 2 (A):^e 3 }. 

Note that -63(A) is a unit circle with holes, centered at the origin, and 63(A) C 62(A) C 
61(A). We define ^(A) as a level set for \( 3 \a, 5e) in the (e 2 k- 1 - 6 ) -neighborhood of 

( A)nonres • 

D 3 (A) := {>?= kv: z7g B 3 (X), |x - x 2 (A, v)\ < e 2 k~ l ~\ A (3) (a, i?) = a} . 

Next two lemmas are to prove that T>% (A) is a distorted circle with holes. Their 
formulations and proofs are analogous to those of Lemmas 13.161 and 13.171 

Lemma 4.7. For every v E 63(A) and every a, < a < 1, there is a unique 
x = xs(X, v) in the interval I3 := [2*2 (A, u) — €2k~ 1 ~ s , X2(A, v) + e2&~ 1 ~' 5 ] such that 

A (3) (a,x 3 i7) = A. (4.20) 

Furthermore, 

|x 3 (A, v) - x 2 (A, v)\ < 2aef fc -1 . (4.21) 

Further we use the notations xz{tp) = X3(A,z7), h 3 [ip) = ^3(99) — %z[fp) = 

Lemma 4.8. The following statements hold for A > k%: 
(1) The set V 3 (X) is a distorted circle with holes: it can be described by the 
formula: 



^s(A) = {xGM 2 : K=K Z [if), ^0 3 (A)}, 



(4.22) 
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where >c 3 ((p) = xi{$>) + h 3 (ip), x-ii^p) is the "radius" of T>2(X) and h 3 {ip) 
satisfies the estimates 

dh 3 



\h 3 \ < Uaejk- 1 , 



< 96aelk 1+s . (4.23) 



dip 

(2) The total length of -83(A) satisfies the estimate: 

L (B 2 \ B 3 ) < 4irk- 4 - 2s ^ 2s2 . (4.24) 

(3) The function h 3 (<p) can be extended as a holomorphic function ofip to the com- 
plex non-resonce set <P 3 and its (k~ s A 2 ^ -neighborhood S 3 , estimates {1^.23 ) 
being preserved. 

(4) The curve T> 3 {X) has a length which is asymptotically close to that ofT>2(X) 
in the following sense: 

l(v 3 (X)) a =^l(p 2 (A))(i + 0(£T 4 - 2si - 2s2 )), (4.25) 

where 0(Ar 2 - 2si - 2s2 ) = (1 + o(l))c A;" 2 " 2si " 2s2 , |o(l)| < 1(T 2 when k> k*. 

Proof. The proof is analogous to that of Lemma 13.171 Note only that, in Part 2, 
when proving convergence of the series for the resolvent (if® (#2 (<£>)) — z) , we use 
the estimate 

H {2 \h 2 {^)) - zY 1 <-, zeC 3 , (4.26) 

analogous to ()3.55j) . the operator iT® acting in ^2(^3)- The estimate ()4.26p follows 
from ()4.9p and (j3.42p . As a side result of these considerations, we obtain an estimate 
analogous to (|3.58p for the new resolvent and z being inside C3. ■ 

We define the non-resonance set, X 3 W m S 3 (X) by the formula analogous to f)3.59|) : 

X* 3 (X):=JC 3 V 3 (X). (4.27) 
The following lemmas are analogous to Lemmas 13.181 and 13.191 

Lemma 4.9. The set X 3 W belongs to the (^ae^k^ 1 ^ -neighborhood of X3W i n 

If i® G X 3 W> then the operator if® (i®) has a simple eigenvalue equal to X. This 
eigenvalue is given by the perturbation series fi4.15\ ). 

Lemma 4.10. Formula (f-^.ffTp establishes one-to-one correspondence between X 3 W 
and D 3 (X). 

Remark 4.11. From geometric point of view this means that X 3 W does not have 
self- intersections. 

4.5. Preparation for the Next Approximation. Let 6® e K 3 and b^ be the 

distance of the point 6® to the nearest corner of K 3 : 

6® = min |fe (3) - 2-Km/N 2 N 1 a\. (4.28) 

m=(0,0),(0,l),(l,0),(l,l) 
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We assume = \b^\. In the case when is closer to a vertex other than 
(0,0), the considerations are the same up to a parallel shift. We consider two cases: 
b ( o ] > e 2 k' 1 - 25 and < b^ < e 2 k~ 1 - 25 . Let 

y^) = i? 2 (</?) + £ (3) . (4.29) 



4.5.1. The case bf ] > e^ 1 - 25 . 

Definition 4.12. We define the set 0( 3 )(&( 3 )) by the formula: 

(3) ( g(3)) = Up(2)eP(2) 0( 2 )(2vr^ 2 )/iV 2 Ar ia + ^(3) )) (4.30) 

f 2 \ ( 2 \ 

set Og being defined by formula f)3.65|) . We assume Og fl #2 ^ 0- 

The set 0^ 3 \b^) consists the disks with the radius r^ = r ( 1 )/ c -2- 4s 2-5 centered 
at poles of the resolvent (^H 2 {jf 2 \ip)^ — k 2 

Lemma 4.13. If<p£& 2 \ C (3) (^ (3) )> then 

H^^i^-eY 1 <!Zl^!! (4.31) 

The estimate is stable in the (r^ k~ s ) -neighborhood of S 2 \0^ (b^) . The resolvent is 

an analytic function of \p in every connected component Oc (b^) of (b^), whose 
intersection with <P 2 is not empty. The only singularities of the resolvent in such a 
component are poles. The number of poles (counting multiplicity) of the resolvent 
inside 0^(b^) is less than ck 2+2s3 . The total size of 0^ 3 >(b^ 3 >) does not exceed 

CoA ,_4- 2si -2s 2 _ 

Proof. The proof is analogous to that of Lemma 13.221 To obtain (|4.3ip we combine 
Lemmas I3T261 and [373T1 for b^ = 2npW /N 2 N x a + 6 (3) , p (2) € P (2) , and take into 
account that e 2 > e 2 . By the same lemmas the number of poles of (^H^ {K 2 (tp) + 

2vp {2 ^lN 2 N x a + $ 3 )) - k 2 } 1 inside 0? ] (2tt P ^ /N 2 N ia + b&) is less than ck 2+2s ' 2 . 
Considering that this resolvent does not have poles outside 0^\2i:p^ /N 2 Nia-\-b^), 
we obtain that the total number of poles of (^H^ ( K K 2 {ip)+2-Kp^ /N 2 Nia+b^) — k 2 ^J 
inside 0^(b^) is less than cok 2+2s2 . Taking into account that the number of 

p (2) G 

p( 2 ^ does not exceed ck 2 ( S3 ~ S2 \ S3 = 2s 2 , we obtain the estimate for the number of 

poles for the resolvent (^H^ ((<p)) — k 2 ^j 

We estimate the size of 0^(b^) the same way as we estimated the size of e>( 3 ). 
Indeed, according to Lemma f3.24l the total size of each of ] (2tt P ^ /N 2 N ia ) is less then 
co£;~ 4-5 ' S2-4<5 . Considering that the number of p^> does not exceed /c 2 ( S3_S2 ) = k 2s2 , 
s 2 = 2si, we arrive at the estimate for 0( 3 )(&( 3 )). ■ 
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-1 



Definition 4.14. We denote the poles of the resolvent \^H^ (if 2 } ((</?)) — k 2 

in 0( 3 )(6( 3 )) as ip { n \ n = 1,...,M^~> < cok 2+2s3 . Let us consider the circles 

On Qffl) of the radius = fc-2-4s3-(5 r (2) aroun( j these poles. Let 

0(3)(6(3)) = U ^ 3 1 ) Oi 3 )(6( 3 )). (4.32) 

Lemma 4.15. The total size of Og (b^) is less then cok^ 2s3 ^ s A 2 \ 

Proof. The lemma easily follows from the formula = /^- 2 - 4s 3-<5 r (2) anc [ ^ e esti- 
mate M( 3 ) < cok 2+2s3 . | 

Lemma 4.16. Ify€& 2 \ of\b^), then 

' H^^^)) - k 2 



< 



4 



(4.33) 



The estimate is stable in the (r^k ^-neighborhood of &2\@s 3 \b^) ■ The resolvent 

an analytic function of if in every connected component oi^(b^) ofO^\b^), whose 
intersection with <&2 is not empty. The only singularities of the resolvent in such a 
component are poles. The number of poles (counting multiplicity) of the resolvent 

inside oi 3 \b^) is less than cok 2+2si . 

Proof. By the definition of oi 3 \b^), the number of poles (counting multiplicity) of 
the resolvent inside this set is less than ck 2+2s3 . Considering as in Lemma 12.291 we 
obtain 



< v 



_ M(3) l7r( 2 h 3S 



(2 



(4.34) 



where v is the coefficient of contraction, when we reduce G^(b^) to oi 3) (b^). 
Namely v is the ratio of r^ to the maximal size of e>( 3 )(fc( 3 )). By Lemma EE3] and 
the definition of A 3 \ v = / c - 4 - 4s 3-4s2-2<5_ Considering that v~ M(A ' r^k s < ej 1 , we 
obtain (fPBI) . " | 

By Lemma f4.151 the total size ofCi 3) (fe( 3 )) is less thenr( 2 ). Since the smallest circle 
m 0( 3 ) has the size r' the function ok%{}p) is holomorphic inside each connected 
component of oi 3 \b^) which has non-empty intersection with $3. Let 



Lemma 4.17. If<pe£ 3 \ of ] (b^), then 

V 3 V 3 V))-^ N " 



2 

^72 



(4.35) 



(4.36) 



The estimate is stable in the (r^k s ) -neighborhood of &$\Os (b^). The resolvent 

is an analytic function of ip in every component of ol 3 \b^), whose intersection with 
^3 is not empty. The only singularities of the resolvent are poles. The number of 
poles (counting multiplicity) of the resolvent inside 0( 3 )(&( 3 )) is less than cok 2+2ss . 
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Proof. The proof is analogous to that of Lemma 13.261 up to the shift of indices by 
one. We use Lemma f4.16( || W3H00 < ^2> an d the first estimate in (I4.23p . | 

4.5.2. The set for small . The considerations of this section is analo- 

gous to those of Section [3.4.21 up to the shift of indices by one. The following lemmas 
and the definitions are completely analogous to 13.271 - 13.311 

Lemma 4.18. If < < e 2 £;~ 1_2<5 and |eo| < b^ k 1 ~ 5 A 3 ^ , then the equation 



(4.37) 



i (3) ~ i (3) i(2) 

has no more than two solutions (p eo in<P 3 . For any (p £0 there is <p £ $2 such 
that 



,^ (3) 

here and below y>Q is ipt ^ for eo : 

r(3) 



±(2) 
<Po 



< r^/A 



(4.38) 



0. 



Lemma 4.19. Suppose < 6q < e 2 £; and <p E <F 3 obeys the inequality analo- 

gous to &3§: ^-^l < ^ (3) . Then, J^\® = fc ^oo ±2b®k(l + o(l)) , 

where |o(l)| < 10 -2 + ei + £2 u>/ien k > k*. 

Definition 4.20. Let (b^ be the open disks centered at (f^^ with radius 

r (3). 7 (3)±^( 



(6 (3) ) be their boundary circles and { p (6 (2) ) = r( 3 )+ U 



Lemma 4.21. For any ip inS 3 \ oP(b^), \X^(y< 3 \ip)) - k 2 \ > &{j 3 V - M 3 ) 
Lemma 4.22. For any tp £ S 3 \ (& (3) ), 

V 3) (^m)-^)" 1 



16 



6 (3) r (3) fc l-5 



(4.39) 



The estimate is stable in the (r^k ^-neighborhood of$3\0^\b^). The resolvent 

is an analytic function of ip in every component of ol 3 \b^), whose intersection with 
^3 is not empty. The only singularities of the resolvent are poles. The resolvent has 
at most two poles inside Os (6®). 

5. The b-th Step of Approximation. Swiss Cheese Method. 

5.1. Introduction. On the n-th step, n > 4, we choose s n = 2s n _i and define the 
operator by the formula: 



H (n) =H (n-l) +aWni (0 <«<!), 



W n 



r=M n -x+l 



where M n is chosen in such a way that 2 Mn ~ k Sn . Obviously, the periods of W n are 
2 Mn ~ 1 (d 1 ,0) and 2 M "~ 1 (0,d 2 ). We write them in the form: iV n _i • ... • iVi(ai,0) and 



iV n _! • ... • iVi(0,a 2 ), here 7V n _i = 2 M "" M ™- 1 , < N n ^ < 4A; S »- S »- 1 . Note 
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that ||W n ||oo < Y^f=M -i+i ll^lloo < exp(— /c^™- 1 ). By analogy with the definition 
of ei, 62 we introduse the notation e n = exp(— k VSn ). 

The n-th step is analogous to the second step up to replacement the indices 3 by 
n, 3 by n— 1, the producr N2N1 by N n _\ ■ ... ■ N\, etc. 

We note that k, satisfying (|3.ip . obeys the analogous condition for with any s n 
instead of s\ 

C*(l + s n )k 2+4Sn In k < k T,8n (5.1) 

with the same constant C*. The inequality (|5.ip can be obtained from ()3.ip by 
induction. This is an important fact: it means that the lower bound for k does not 
grow with n, i.e., all steps hold uniformly in k for k > k*, k* being introduced by 
(|3.ip . Further we assume k > k*. 

The formulation of the geometric lemma for n-th step is the same as that for Step 
2 up to a shift of indices, we skip it here for shortness. Note only that in the lemma 
we use the set Xn-iW t° define Xn(A). In fact, we started with the definition of 
Xi(A) and then use it to define Xi(^) (Step 1). Considering X1W1 we constructed 
X2(A) (Lemma I3.5P and later used it to define X2W (Section 13. 3p . Using X2M' we 
introduced X3(^) (Lemma l4.ip and, then x|(A) (Section l4.4p . Thus, the process goes 
like xi - > Xi — ► X2 — > X2 ~ * X3 ~~ > xt- Here we start with the set Xn-iW defined by 
(|4.27p for n = 4 and by (|5,15p for n > 4. @ The estimate (|4.6p for n-th step takes the 
form: 

LiJCnX^Axn)) <fc -5„ ? Sn = 2j2(l + Sl ). (5.2) 

L \Xn-l) i=1 

The formulation of the main results (perturbation formulae) for n-th step is the same 
as for the second and third step: Theorems analogous up to the shift of indices to 
Theorem I3,8[ / |4~4"1 Lemma 13.114 Theorem 13.121 Corollary 13.131 and Lemma 13.141 hold. 

5.2. Proof of Geometric Lemma. The proof of the lemma is analogous to that 
for Geometric Lemma in the second step. Let us consider 

0( " , = U P<"-'»^"-'AM°^ 1> (i|S^)' < 5 ' 3 > 
where is defined by Definition MM with = 2 */*l w h en 



n = 4. When n > 4 we use Definition 15.61 with n — 1 instead of n and take 6( n-1 ) = 
■jt?^ — rf—. The radius r^" -1 ) of Os ^ is defined by the recurrent formula: r*" -1 ) = 

iv n _l...iVia 15 •> 

^.-2-4s„_i-<5 r (n-2)_ jj^g means 

r (n-l) = fc -2n-4EfcI 1 1 ^-2si-25 _ fc -2n-2££ =1 s„-2<5 = k -S n+1 -28 ^ /g^ 



^Strictly speaking we assume that there is a subset Xn-iM °f the isoenergetic surface 5 n -i(A) 
of H^ n ~^ such that perturbation series of the type (|4.15p . 1)4.161) converges for t^ n ~ 1 ' € Xn-i(A) an( i 
X*i_i(A) has properties described in Section f4. 41 up to replacement of 3 by n — 1. In particular, we 
assume that Xn-iW = ^n-iI'n-i(A), where X> n _i(A) satisfies the analog of Lemma 14.81 and that 
the analogs of Lemmas 14. 91 and 14. 101 hold too. Here, IC n -i is the parallel shift into K n -i. Further in 
this section we describe the next set Xn(A) which has analogous properties. 
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> 



S n being denned by (|5.2p . Note, that the definition make sense, since |^ — - — 

dmaxk~ Sn > tn-ik ~ 1-25 , when k > fc*, the estimate (|5.1h has been used. The last 
inequality can be easily proved by induction. Let 

# n = £ n _! \ 0<- n \ 9„ = <2> n n [0, 2vr), (5.5) 

3>3, ©3 are given by (j4~8l) . By Lemmas liTfl |5"U1 

\ -l 



i^ (n - 1) (>?n-lM + 



2vrp( n - 1 ) 
iV„,_i...A r ia 



k 



< 



-n-2 



(5.6) 



for all p( n 1 ) g p( n *) \{0} and <p G S n , here and below is A n ''-neighborhood 
of & n . +We consider V n - 1)nonres C X> n _i: 

T^n-l,nmres = {^h-lOp)) V 3 G @n} • (5.7) 
We define x« by the formula: 

Xn — K-'n'Dii— 1 .nonre.s • (^-^) 

By the definition of IC n , for every in Xm there are G P^,... £ 
and j G Z 2 such that 

2vrp("- 1 ) 27rpW 2ttj _ . . _ . . . . 

"TJ 77 h ••• + "^7 1 = ^n-UV 3 )) Jin-lVP) e Ai-l.nonres- (5.9) 

7V n _1...7Via Wia a 



Considering the definition of £>n_l,nonres an d (|5.6p . we obtain: 

x -l 



H {n-l)Un) + 



27rp 



(n-l) 



N n -!...N ia 



< 



c n-2 



(5.10) 



~(n-l) 



-j-pC™- 1 ), Note that the indices j, p^,...,^" --1 ) do not play a role, 
since they just produce a shift of indices of the matrix elements of the resolvent. Con- 
sidering that p( n_1 ) can be any but zero, we obtain that (|5.10p holds for all p( n ~ 1 ') £ 
p(n-i) \ Taking into account that aJ* -1 ^") + 27rp( n ~ 1 ViV„_i...iVia) = k 2 

and e„_i < e 2 _ 2 . We arrive at the analog of (|4.5p for all p( n_1 ) ^ p( n— -v.This 



also proves that p( n ^ is uniquely defined by (|5.9|) . It remains to check the analog 
of for pt"- 1 ) = p^" 1 ). Let iC"" 1 ) := iW + 2vrp( n - 1 )/iV„_i...iVia. By (Q, 

t(n-i) e K, n _xD n _ x . Using (|43T|) for n = 4 and (f5TT5|) with n-l instead of n for 
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> 4, we get rj(™ _1 ) € Xn— i- By the analog of TheoremlOlfor step n-l. Xy j (rj(™) + 
2 7r p( n - 1 )/JV n _i...iV 1 a) is the only eigenvalue of ijt 71-1 ) (t( n ) + 2vrp("- 1 )/iV ri _i...Afia) 
in the interval £ n -\- Hence, 

| A (— D^(3) + 2vrp^ 1 )/iV„_ 1 ...iV 1 a) - X^ (t^ + 2irp^ /N n _ x ...N x a)\ > e n _ 2 . 

Thus, the analog of holds for all p^ n ~^ G P^" 1 ). Part 2 holds, since all estimates 
are stable with respect to the perturbation of less then e n -\k~ l ~ & . 

Let us estimate the size of O n . According to Lemma 14.151 15.71 the total size of 
each 0i n ~ 1) (27rp(™- 1 )/iV„_i...iVia) is less then c £;- 2s "- 1 ~M n - 2 ). Considering that 
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the number of p( n_1 ) does not exceed k 2 ( s n~ s n-i) _ k 2s n-i ) Sre _ 2s n ._i, we obtain 
that the size of O n is less than k~ s r^ n ~ 2 ^ . Using the formula for A n ~ 2 \ we obtain that 
the total size of O n does not exceed fc~ 2 ("~ 1 )~ 2 ^fc=i Sfc . Using this estimate we easily 
arrive at (|5.2p . It is easy to see that S n = 2(n — 1) + (2™ — 2) si and S n « 2 n si « s n 
for large n. The lemma is proved. 

5.3. Nonresonant part of the isoenergetic set of Now we construct a 

nonresonance subset XnW °f the isoenergetic surface S n (X) of in K n , S n (X) C 
i^ n . It corresponds to nonresonance eigenvalues given by perturbation series. The sets 
Xi(X), X2(^)) XtW are defined in the previous steps as well as the non-resonance sets 
Xi(A), X2(A), X3(A)- Let us recall that we started with the definition of xi(A) and 
then use it to define X*(A)- Considering X*(A), we constructed X2(A) (Step 2). Next, 
we defined x^W- Using X2W1 we introduced X3(A) an d, then X;j(A)- Thus, the 
process looks like xi ~ > X* — ^ X2 - >■ x\ ~ * X3 — * X\- The geometric lemma in this 
section gives us Xi an d every next Xn if Xn-i * s defined. To ensure the reccurent 
procedure we show now how to define x* n W using Xn(A). 
We define B n as the set of directions corresponding to 6 n : 

B„(A) = G 5i : p G e n }. 

Note that B n is a unit circle with holes centered at the origin and B n (X) C B n -±(X). 
We define T> n (X) as a level set for X^ n \a, i?) in the (e n _ifc~ 1_<5 )-neighborhood of 

T^n— 1 ,nonres ( A) • 

^n(A) = I i? = xz7 : £7 g B„, |x - x n _i(A, z7)| < e n _iA; _1_<5 , A (n) (a, x) = a| . 

Considering as in the previous step, we prove the analogs of Lemmas 14.71 and 14.81 
For shortness, we provide here only the second lemma. By analogy with previous 
sections, we shorten notations here: x n -\{ip) = x n _i(A, u), K n —\{^p) = K n —i(X,v)v. 

Lemma 5.1. For A > A; 2 ; 

(1) The set T> n (X) is a distorted circle with holes: it can be described by the 
formula: 

P„(A) = {xGM 2 :H = a n (ip), <p G 6 n (A)}, (5.11) 
where *c n (<p) = >c n -i((p) + h n ((p), and h n {ip) satisfies the estimates 

\h n \ < \2cte n _ik 

(2) The total length of B n (X) satisfies the estimate: 

L (B n -i \ B n ) < Airk~ Sn . (5.13) 

(3) Function K n {ip) can be extended as a holomorphic function of ip to <P n , esti- 
mates \5.12§ being preserved. 

(4) The curve T> n (X) has a length which is asymptotically close to that o/X> n _i(A) 
in the following sense: 

L(v n {X)) A =^L(p n _!(A))(l + 0(£r 5 ")). (5.14) 



dh r 



dip 



l+<5 



(5.12) 
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Now define the nonresonance set, XnW m S n (X) by the formula analogous to 
HHZZH . Indeed, 

X* n {\) := IC n V n (X). (5.15) 
The following lemmas are analogous to Lemmas 14.91 and 14.101 

Lemma 5.2. The set XnW belongs to the (2ae^_j_A; J -neighborhood of Xn(^) * n 

K n . If t^ G XnW> then the operator H^\t^) has a simple eigenvalue equal to X. 
This eigenvalue is given by the perturbation series analogous to \4- 15\ ). 

Lemma 5.3. Formula A5.15\) establishes one-to-one correspondence between XnW 
and D n (X). 

5.4. Preparation for the Next Approximation. Let G K n and b^ be the 



distance of the point &( n ) to the nearest corner of K n . We assume b^ = |&( n )|. We 
consider two cases: > e n -\k~ x ~ 2& and < 6q < e n ^\k~ l ~ 26 . Let 

^ n - 1 \tp) = K n _ 1 {tp) + b^. (5.16) 
5.4.1. The case b^ ] > e n -ik- 1 - 25 . 

Definition 5.4. We define the set 0^ n \b^) by the formula: 

C>W(&(«)) = U p( n-i) e p(„-i)Oi n -^ (2Trp^- 1 )/N n ^ 1 ...N 1 a + , (5.17) 

set oi n ^ being defined by Definition 15.41 for n = 4 and by 15.61 for n > 4. 

The set 0^{b^) consists the disks with the radius r^" 1 ) = r (n-2)£-2-4 Sn _i-5 
centered at poles of the resolvent ^-ff n -l (y*™ -1 ) (tp)) — k 2 ^j 

Lemma 5.5. If tp G $ n -\ \ C (n) (6 (n) ), then 

H^^-'Xp)) - k 2 Y X < 17rin ~ 1)k3S (5 . 18) 

' £rt— 1 

The estimate is stable in the (r( n-1 )/c -<5 ) -neighborhood of $ n —l \ 0^ n \b^ n )). The 
resolvent is an analytic function of tp in every connected component 0^\b^ n )) of 
0( n \b( n )), whose intersection with ^ n -i is n °t empty. The only singularities of the 
resolvent in such a component are poles. The number of poles (counting multiplicity) 
of the resolvent inside is less than ck 2+2s " . The total size of 0( n )(6( n )) 

does not exceed cok~ Sn . 

Proof. The proof is analogous to that of Lemma 13.221 Let n = 4. To obtain (|5.18p 
we combine Lemmas IH3 and H22 for = 2vrp( n " 1) /^-i---A r ia + & (n) taking 

into account e^_ 2 > e n -i and the estimate for 6q . If n > 4, then using the recur- 
rent procedure, we apply Lemmas 15.91 and 15.141 with n — 1 instead of n. Moreover, 

the number of poles of (^H^ 1 ) (x n _i(ip) + &( n_1 )) — A; 2 ^ inside 0( n )(6 ra ) is less 

than ck 2+2,Sn ~ 1 . Considering that the number of p( n_1 ) g p( n_1 ) does not exceed 
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ck 2 ( Sn Sn ~ 1 \ s n = 2s n -i, we obtain the estimate for the number of poles for the 
resolvent (h^- 1 ) (y (n_1) ((¥')) - • We estimate the size of 0^(b^) the same 

way as we estimated the size of 0^ n \ | 



Definition 5.6. We denote the poles of the resolvent yH^ 71 ^ (%f n ^{(pj) — k 2 

in 0( n )(6(™)) as tp$ , m = 1,...,M^ AfW < c k 2+2s ". Let us consider the circles 
Om\b^) of the radius r^ = / c - 2 - 4s ™-< 5 r O-i) around these poles. Let 

0W(feW) = u£ ) 0W(feW). (5.19) 

Lemma 5.7. The total size ofOs(b^) is less then co£r 2Sn ~M n_1 ) . 

Proof. The lemma easily follows from the formula r^ = / C " 2 - 4s «- <5 r ( n - 1 ) anc l the 
estimate JlfW < c k 2+2Sn . | 



Lemma 5.8. If tp G £ n -i \ oi n \b^), then 



< 



c n-l 



(5.20) 



The estimate is stable in the (r^k s ) -neighborhood of & n -i \ o\ ,(b^). The re- 
solvent is an analytic function of tp in every connected component Oj£ (b^) of 

0a (£>( n )), whose intersection with <& n -i is n °t empty. The only singularities of the 
resolvent in such a component are poles. The number of poles (counting multiplicity) 

of the resolvent inside oi n \b^ ) is less than ck 2+2Sn . 

Proof. By the definition of oi n \b^), the number of poles (counting multiplicity) of 
the resolvent inside this set is less than = cok 2+2s " . Considering as in Lemma 

12.291 we obtain 



ff(n-l) (^«-l)(^)) - A; 2 



< v 



17r (n-l) k 3S 



£n-l 



(5.21) 



where v is the coefficient of contraction, when we reduce 0^ n \b^) to oi n) (b^). 
Namely v is the ratio of r^ to the maximal size of 0^(b^). By Lemma[53]and 
formula (15.41) . v = ^- 4 - 4s «- 4s n-i. Considering that r^ n ~^k s < 1 and using (15. ip , 
we obtain v~ M r { n ~ l )k s < when k > k*. Thus, we obtain (|5.20[) . 

■ 

By Lemma 15.71 t ne total size of oi n \b^) is less then r( n_1 ). Therefore, the 

function >c n {p) is holomorphic inside each connected component of oi n \b^) which 
has non-empty intersection with <£> n . Let 

yi n \<p) = x n (<p) + $ n \ (5.22) 
Lemma 5.9. IfpeS n \ Oi n \b^), then 



< 



(5.23) 



-n-l 
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The estimate is stable in the (r^k - ^ -neighborhood ofS n \Og(b^). The resolvent 

is an analytic function of if in every connected component oi^\b^) of Og (b^), 
whose intersection with ^> n is not empty. The only singularities of the resolvent 
in such a component are poles. The number of poles (counting multiplicity) of the 
resolvent inside 0^ n \b^) is less than cok 2+2s " . 

Proof. The proof is analogous to that of Lemma 13.261 up to the shift of indices by 
two. We use Lemma [5\8| ||Wn||oo < e^-ii an d the first estimate in (I5.12D . | 

5.4.2. The set {b^\ for small b^ . The considerations of this section is analo- 
gous to those of By analogy with the previous subsection, we choose = ^- 2 - 2s «^ <5 r ( n ^ 1 ) 
here. The following Lemmas 15.101 and 15.111 are identical to 14.181 and 14.191 up to re- 
placement of indices (3) by (n), (2) by (n — 1), etc. Next, Definition 15. 121 is analogous 
to KM Lemmas 15331 lOH are analogous to ICTl K22\ 

Lemma 5.10. If < b$ < e n _i/c~ 1-2<5 and | e 1 < b^ k l ~ s r^, then the equation 

\ {n) (f n \p))=k 2 + e (5.24) 

has no more than two solutions p>% ^ in S n . For any pt ^ there is Pq^"" ^ £ ^n-i 
such that 

' <t»/4, (5.25) 



here and below ip^ ^ is ^ for eo = 0. 

Lemma 5.11. Suppose < b^ < e n -\k~ l ~ 25 and ip £ S n obeys the inequality anal- 



±(n-l) 
P~Po 



< rW. Then, J^\W (y^M) = k ^ ±2b^k(l + 
< W~ 2 + ei + ... + e n -i when k > k*. 



ogous to 115. 25]) : 
o(l)), where |o(l) 

Definition 5.12. Let (R n A be the open disks centered at <Pq w ^h radius 

r (n). 7 (n)±^(„)^ be thdr boundary drcles and {n) (£(n-l)) = p(n)+ y p(n)-_ 

Lemma 5.13. For any p inS n \ oi n \b^), \\( n) (tp)) - k 2 \ > b^k 1 ' 5 ^. 
Lemma 5.14. For any tp £ # n _i \ d n) (R n A , 

(H^Uy^i^-k 2 )- 1 < 16 (5.26) 
V V J J b^'r^k 1 - 6 

The estimate is stable in the [r^k - ^ -neighborhood of ( P n \oi n \b^) . The resolvent 

is an analytic function of ip in every component of oi n \b^), whose intersection with 
<P n is not empty. The only singularities of the resolvent are poles. The resolvent has 
at most two poles inside oi n \b^). 
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6. LlMIT-ISOENERGETIC SET AND ElGENFUNCTIONS 

6.1. Limit-Isoenergetic Set and Proof of the Bethe-Sommerfeld Conjec- 
ture. At every step n we constructed a subset B n {X) of the unit circle, and a function 
x ra (A, v\ V G B n (X), with the following properties. The sequence B n (X) is decreasing: 
B n {\) C B n -\(\). The set V n {\) of vectors 5k = xr n (A,z?)z?, G -8„(A), is a slightly 
distorted circle with holes, see FigJU Figj2j formula (|1.15j) and Lemmas 12.101 13.171 
14.81 15.11 For any 5k n (X, v) G T> n (X) there is a simple eigenvalue of H (jl \5k n ) equal to 
A and given by a perturbation series. Let B ao (X) = f^ =1 B n (X). Since B n+ \ C B n 
for every n, Boo(X) is a unit circle with the infinite number of holes, more and more 
holes of smaller and smaller size appearing at each step. 

Lemma 6.1. The length of Boo(X) satisfies estimate hi. 9^ with 73 = 5. 

Proof. Using (I277TD . ([333]), (g23| and (EU3l) and consideirng that S n « 2 n si, we 
easily conclude that L(B n ) = (l + 0(/c~ <5 )), = A 1 / 2 uniformly in n. Since B n is a 
decreasing sequence of sets, (|1.9p holds. ■ 

Let us consider Xoo(A, z?) = lim n _ 5>00 x n (A, z?), z? G ,600 (A). 

Lemma 6.2. T/ie /imii ^(A, z?) exists for any V G i3oo(A) and the following estimates 
hold when n > 1: 

|x 00 (A, I 7)-x n (A, l 7)| <14e^- 1 , e n = exp(-^ s "), s n = 2^ 1 s 1 . (6.1) 

Corollary 6.3. For every z? G £>oo (A) estimate il.lO\) holds, where 74 = 3 — 30si — 
205 > 0. 

Proof. The lemma easily follows from the estimates (|3.52p . (j4.23p and (j5. 12|) . To 
obtain corollary we use (]2,70p , ■ 

Estimates (|3"32l . (1031) and (I5TT21) justify convergence of the series ^2m=i ^dtp ' 
hence, of the sequence We denote the limit of this sequence by ^§ 2S -- 

Lemma 6.4. The estimate \1.18\) with 75 = 1 — 33si — 226 > holds for any 
ugBoo^X). 

Proof. The lemma easily follows from (|2T70]) . (13321) . (|4T23|) and ([5A2]) . | 

We define P^A) by (|1.8p . Clearly, Poo(A) is a slightly distorted circle of radius 
k with the infinite number of holes. We can assign a tangent vector ^f^-v + Xqo/2, 
/2 = (— sin (p, cos to the curve 1500(A), this tangent vector being the limit of corre- 
sponding tangent vectors for curves V n {X) at points 5k n (X, v) as n — > 00. 

Remark 6.5. We easily see from (|6,ip . that any 5k G T> !X} {X) belongs to the (l4e^ -1 )- 
neighborhood of T> n (X). Applying perturbation formulae for n-th step, we eas- 
ily obtain that there is an eigenvalue X( n '(5k) of H^ n \5k) satisfying the estimate 

^The operator H < - n \>c) is defined for every % £ R 2 . The perturbation series is given by a formula 
analogous to (I3.44|) . which coincides with (|3.29|l up to a shift of indices corresponding to the parallel 
shift of x into K n . 
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A(")(x) = A + 5 n , 5 n = O (4), the eigenvalue A^(x) being given by a perturba- 
tion series of the type (|3.44|) . Hence, for every x e Poo(^) there is a limit: 

lim A^(x) = A. (6.2) 

n— >oo 

Theorem 6.6 (Bethe-Sommerfeld Conjecture). T/ie spectrum of operator H contains 
a semi- axis. 

Proof. By Remark [631 there is a point of the spectrum of H n in the ^-neighborhood 
of A for every A > k 2 , being introduced by (|3,ip . Since ||iT n — < e„, there is 
a point of the spectrum of H in the ^-neighborhood of A, <5* = 5 n + e*. Since it is 
true for every n and the spectrum of H is closed, A is in the spectrum of H. | 

6.2. Generalized Eigenfunctions of H. A plane wave is usually denoted by e i ^ k ' x ', 
keM 2 . Here we use x instead of k to comply with our previous notations. We show 
that for every it in a set 

Soc = U A > At P 00 (A), A* = A; 2 , (6.3) 
there is a solution ^(^(x, cc) of the equation for eigenfunctions: 

-A¥ 0O (x,x) + V(x)V oo (5i,x) = A 00 (i<)^ 00 (i<,rE), (6.4) 
which can be represented in the form 

^> OQ {x,x) = e i{ ^ x) (l + u OQ {5i,x)^ |hoo(^,^))|| Loo(K 2) <c(F)|x|- 71 , (6.5) 

where Uoc(x,x) is a limit-periodic function, 71 = 1/2 — 15si — 85; the eigenvalue 
Aoo(x) satisfies the asymptotic formula: 

A oc (i?) = |x| 2 + 0(|i?r 72 ), 72 = 2 - 30si - 205. (6.6) 

We also show that the set satisfies (\l.7\) . 

In fact, by (|6.1|) . any x G T> 00 (\) belongs to the (e n /c~ 1_5 )-neighborhood of V n (X). 
Applying (p39l) . ([230]) with 2/3 = 2/3 = 1 - 15si - 95 and the perturbation formulae 
proved for next steps, we obtain the following inequalities: 

H^W^)-^ )^)^ < c||Wi||AT 1 / 2+12si+8<5 , WE^+^iayE^iH)^ < A8e 3 n , n > 1, 

(6.7) 

|A (1) (x)-|i<| 2 | < C(iyi)£r 2+30si+20 ' 5 , |A (n+1) (x)-A (n) (x)| < 12e 4 n , n > 1, (6.8) 

where E^ n+l \ are one-dimensional spectral projectors in L2(Q n +i) correspond- 
ing to potentials W n+ i and W n , respectively; \^ n+l \>c) is the eigenvalue correspond- 
ing to E^ n+l \>i), E^\k) corresponds to V = and the periods 01, a%. This means 
that for properly chosen eigenfunctions \l/ n -|-i(x, x): 

" *o|U 2 ( Ql ) < c[| W x [| ^-1/2+12^+85 1 Qi |a/ 2 ^ ^ o{x) = e i(2, x )^ ( g Q) 

ll*n+l " *nllL 2 (Q n+1 ) < WOef, \Q n+1 \ l ' 2 , (6.10) 
where is \& n extended quasi-periodically from Q n to Q n +i- Eigenfunctions ^> n , 
n > 1, are chosen to obey two conditions. First, H^n II L 2 (Q„) = |Qn| 1//2 ; second 

10 The condition ||*»|U 2 (q„) = \Q„\ 1/2 implies ||*„||l 2 (q„ +i) = \Q n +i\ 1/2 ■ 



68 Y. KARPESHINA AND Y.-R. LEE 

(# n , #„_i) n > 0, here (•,•) n is an inner product in L2(Q n )- These two conditions, 
obviously, provide a unique choice of each ty n . Considering that ^n+i and *ff n sat- 
isfy equations for eigenfunctions and taking into account (|6.8p . (|6.10p we obtain: 
||*n+i - ^n\\wi(Q n+1 ) < 2fc 2 e^|Qn+i| 1/2 , n > 1, and, hence, - ^n\\ Loa (Q n+1 ) < 

6k 2 e 3 n \Q n+l \ 1 / 2 . Since ^n+i an d ^ n obey the same quasiperiodic conditions, the same 
inequality holds in the whole space M 2 : 

Pth-i - Vnh^m < Qk 2 el\Q n+1 \ 1 ' 2 , n>l, (6.11) 

where ^n+ij are quasiperiodically extended to M 2 . Obviously, we have a Cauchy 
sequence in L oc (M 2 ). Let ^' 00 (x, x) = lim n _ Sl00 ^ n (x, x). This limit is defined point- 
wise uniformly in x and in W 2 loc (M?). 

Theorem 6.7. For every sufficiently large X, A > AoQ^Li ll^ll)^)* an< ^ ™ ^ ^oo(A) 
i/ie sequence of functions ^ n (x,x) converges in Loo(M 2 ) and W" 2 2 , oc (M 2 ). The limit 
function ^^(x^x), * 00 (i?, x) = lirrin-^oo ^f n (>e, x), satisfies the equation 

- A* 00 (i?, x) + y(x)^oo(i?, x) = A^oo(x, x). (6.12) 

It can be represented in the form 

^oo(x,x) = e i( ^ ) (l + u 00 (i?,x)), (6.13) 

where Uoo(>i,x) is a limit-periodic function: 

oo 

u^k.x) = ^£t n (x,x), (6.14) 

n=l 

u n (H,x) being a periodic function with the periods 2 Mn ~ 1 d\,2 Mn ~ 1 d2 with 2 Mn ~ 

WMl^) < c(V)k-^, 7l = 1/2 - 15 Sl - 85, (6.15) 
IMiooCR') <Qk 2 el_ 1 \Q n \ 1 ' 2 , n>2. (6.16) 
Corollary 6.8. Function n^x, x) obeys the estimate W.5\) . 

Proof. Let us show that is a limit-periodic function. Obviously, = ^0 + 
J2^=o(^ n+i — ^n), the series converging in L 0O (lR 2 ) by (|6.1ip . Introducing the no- 
tation u n+ i = e~ % ^ ,x \^ n+ i — ^ n ), we arrive at (|6.13|) . (|6.14p . Note that u n is 
periodic with the periods 2 Mn ~ 1 d\, 2 Mn ~ 1 d2- Estimate (|6.16p follows from ([6. lip . 
We check (|6.15p . Indeed, by (|6.9I) . Fourier coefficients (u\) m , m € Z 2 , satisfy the 
estimate |(«i) m | < c(F)A;- 1 / 2+12si+8 ' 5 |Qi| 1 /2 < c (y) fc -i/2+i3 S1 +85_ We win uge thig 

estimate for m : \m — j\ < k Sl . Next, we obtain a stronger estimate for other m- 
s. Indeed, the inequality \E^(x) mj \ < ( c (y) / fc-( 1 /2-i2s 1 -75)y m -J'l /? o 1 follows from 
(|2.25p and (|2.28p . Hence, similar estimates holds for Fourier coefficients of ^i: 

|(%)m| < (c(y)A;^( 1 / 2 ~ 12si_7<5 )) ' m . Summarizing the inequalities, we obtain 
that (|6.15p holds. It remains to prove (|6.12p . Indeed, \P n (i£,x), n > 1, satisfy equa- 
tions for eigenfunctions: H n ^ n = A( n )(x)\I> n . Considering that \I> n (x, x) converges 
to x) in W 2 loc and relation (|6.2p . we arrive at (|6.12p . | 
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Theorem 6.9. Formulae (uO), $6. 5|) and $6. 6\) hold for every x G C/oq. TTie se£ Q, 



is Lebesgue measurable and satisfies U. 7|) wnt/i 73 = 5. 



Proof. By Theorem 16.71 (|6.4p . (|6.5p hold, where Aoo(x) = A for x G Poo (A). Using 
(jl.lOP , which is proven in Corollary 16.31 with = |x| , we easily obtain (|6.6p . It 
remains to prove (jl.7p . Let us consider a small region C/ n (Ao) around an isoenergetic 
surface £>„(A ), A > k%. Namely, U n (X ) = U| A _ Ao | <Pn P n (A), p n = e n ^ik~ 2& ', /c = 

Aq^ 2 . Considering an estimate of the type (|2.62p for A( n )(x), which holds in the 
(e n _iA; _1_2 ' 5 )-neighborhood of P n (Ao), we see that £/„(Ao) is an open set (a distorted 
ring with holes) and the width of the ring is of order e n -i /c~ 1 ~ 2 ' 5 . By Part 4 of Lemmas 
I2TU1 I3T71 S3 EH the length of £>„(A ) is 2nk (l + 0(Ar 5 )). Hence, |C4(A )| = 
2irkp n (l + 0(/c~ 5 )). It easily follows from Lemma 15. II that U n +\ C U n . Definition of 
^oo(Ao) yield: Poo(A ) = n~ =1 f7„(A ). Hence, = C\™ =1 Q n , where 

Gn = Ua >a,^(A ) = U A>A ,_ Pn(M P„(A). (6.17) 

Considering that U n+ \ C U n for every Ao > A*, we obtain Q n +i C Q n . Hence, 
\Goo n Br I = lim n _ ) . 00 |£/ n n Br|. Summarizing volumes of the regions U n , we easily 
conclude \Q n n BrJ = |Br| (l + 0(i?~ 5 )) uniformly in n. Thus, we have obtained 
([LTD with 73 = 5. ■ 



7. Proof of Absolute Continuity of the Spectrum 

The proof is completely analogous to that for the case I > 6. Here we give only 
the list of lemmas and the main result. For details, see |25j . 

7.1. Projections E n {Q' n ), Q' n C Q n . Let us consider the open sets Q n given by (|6.17p . 
There is a family of Bloch eigenfunctions ^(x, x), X G £7 n , of the operator H^ n \ 
which are described by the perturbation formulas (|1.12p . Let Q' n be a Lebesgue 
measurable subset of Q n . We consider the spectral projection E n (Q' n ) of H^ n \ corre- 
sponding to functions *„(x,x), x G ^. By [26], £ n (Q' n ) : L 2 (M 2 ) -> L 2 (M 2 ) can be 
presented by the formula: 

E n {Q'n) F = ^ f (F,* n (x))*„(x)dx (7.1) 

for any F G C^°(M 2 ), here and below (•, •) is the canonical scalar product in L 2 (M 2 ), 
i.e., 

(F, = f F(x)* n (Z,x)dx. 

The above formula can be rewritten in the form: 

En {Q'n) = S n {Q' n ) T n {Q' n ) , (7.2) 

T n : C7 °°(M 2 ) L 2 (^) , S n : L M (&) -> L 2 (M 2 ), 

T n F=^-(F,* n (x)) for any F G C -(M 2 ), (7.3) 
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T n F being in (g' n ), and, 

S n f = ^-f x)dx for any /eL M (G' n ). (7.4) 

By [26], ||T n .F[| L2 (g/j < ||F|| i2(R 2) and \\S n f\\ L2 ^ < \\f\\ L2 (g> n )- Hence, T n , S n can 
be extended by continuity from C^°(M. 2 ), (Q' n ) to L 2 (M 2 ) and L 2 (Q' n ), respectively. 
Thus, the operator E n (Q' n ) is described by (|7.2p in the whole space Z^O^ 2 )- 

Obviously, for every x in Qm there exists a pair (A n ,^) such that X n — A^ n ^(j<) 
and that (cos (p, sin 99) = j||r. Let us introduce new coordinates (X n ,ip) in A ra = 

X (n \x), (cos 99, sin ip) = ^. 

Lemma 7.1. Every point x in Q n is represented by a unique pair (X n ,ip), X n > A*, 
if £ [0,2ir), such that 

x(\ n ,ip) = x n (\ n ,u)v, v = (cos v?, sin v?), (7.5) 
x n {Xn,v) being the "radius" of the isoenergetic curve'D n {X n ) in the direction v. where 

For any function f(x) integrable on Q n , we use the new coordinates and write 

/ f{x)dx= I x(Gn,>c) f(x)dx 
Jg n Jr 2 

x(Gn, y)) /(>?(A n , ¥>)) ^ gx n —d<p d\ n , 

where \ (Ci n ,x) is the characteristic function on Q n , x(X n ,(p) is given by (17.5P and 
d£ = (yxW(x),v)\z =MXntff) . Let 

G n ,x = {ZeG n :X^(x)<X}. (7.6) 

This set is Lebesgue measurable, since Q n is open and X^ n \x) is continuous on Q n . 

Lemma 7.2. \G n ,x+e \ Qn,\\ < 27re when < s < 1. 

By ([7TD, E n {G n ,\ + e) - E n (Q n ,x) = E n {Q n ,\+s \ Qn,\)- Let us obtain an estimate 
for this projection. 

Lemma 7.3. For any F G Cg°(R 2 ) and < e < 1, 

\\{E n (g n>x+£ ) -E n (g n>x ))F\\ 2 L2m < C(F)e, (7.7) 
where C(F) is uniform with respect to n and A. 

7.2. Sets t/oo and £/oo,A- The sets g^, g n are given by (|6.3|) . (|6.17|) . As it was shown 
in the proof of Theorem 16. V\ g n +i C g n , £?oo = fl^Li Qn- Therefore, the perturbation 
formulas for X^ n \x) and ^ n (x) hold in g&o for all n. Moreover, coordinates (A n , y?) 
can be used in g^ for every re. Let 

Goo,x = G ^00 : Aoo(x) < A} . (7.8) 

The function Aoo(i?) is a Lebesgue measurable function, since it is a limit of the 
sequence of measurable functions. Hence, the set g o0; x is measurable. 
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Lemma 7.4. The measure of the symmetric difference of two sets Goo,x and G n ,x 
converges to zero as n —> oo uniformly in A in every bounded interval: 

lim \Goo,x^G n ,x\ = 0. 

n— >oo 

7.3. Spectral Projections E{Q O0 \). In this section, we show that spectral projec- 
tions E n (Goa,\) have a strong limit E O0 (G O o, a) i n -^(R 2 ) as n tends to infinity. The 
operator E 00 (G o,x) is a spectral projection of H. It can be represented in the form 
Eoo(Qoo,x) = SooToo, where Soo and are strong limits of S n (Goo,x) and T n (Gaa,\), 
respectively. For any F G Cq°(R 2 ), we show: 

Eoo {Qoo,\) F = ^ I (F, (i<)) l-oo (x)dx, (7.9) 
HE 00 (G 00 ,x)F = ^ f A 00 (x)(F,^ 00 (i<))* 00 (i<)dx. (7.10) 

Woo, A 

Using properties of £"00 (£/oo, a) > we prove absolute continuity of the branch of the 
spectrum corresponding to functions ^^(x). 

Now we consider the sequence of operators T n (Goo,x) which are given by (|7.3p and 
act from L2OR 2 ) to ^2(^00, a)- We prove that the sequence has a strong limit and 
describe its properties. 

Lemma 7.5. The sequence T n (Goo,x) has a strong limit T OQ (Goo,x)- The opera- 
tor T 00 (Q 00) x) satisfies \\Too\\ < 1 and can be described by the formula TooF = 
g^tfooC*)) for any F G C^°(M. 2 ). The convergence of T n {Q^ x )F to T^G^F 
is uniform in A for every F G L<2(M?). 



Now we consider the sequence of operators S n (Goo,x) which are given by (|7.4p with 
Q'n = ^00, a: 

S n (aoo,A) : L 2 (&o,a) -> ^2(M 2 ). 
We prove that the sequence has a strong limit and describe its properties. 



Lemma 7.6. The sequence of operators S n (Goo,x) has a strong limit S 00 (Goo,x)- The 
operator Soo {Goo, x) satisfies ||Soo|| < 1 and can be described by the formula 



(Soof)(x) = J- / /(i^^ax (7.11) 



/or any / G (£oo,a)- convergence of S n (G 00, \)f to S 00 (G o,x)f is uniform in 

A /or every / G L 2 (£00) ■ 

Lemma 7.7. Spectral projections E n (Goo,x) have a strong limit E 00 (G 00 ,x) in L2(M 2 ), 
i/ie convergence being uniform in A for every element. The operator E OQ (G C xi,x) is a 
projection. For any F G Cq°(]R 2 ) it is given by ( |7.9[ ) and formula ( |7. i0[ ) holds. 

Lemma 7.8. There is a strong limit E oc (Goo) 0/ £/ie projections E^ (Goo, x) a $ A goes 
io infinity. 

Corollary 7.9. T/ie operator E^Gco) o projection. 

Lemma 7.10. Projections E 00 (Goo,x}} A G R, and E 00 (G o) reduce the operator H. 
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Lemma 7.11. The family of projections -Eoc(£/oo> A) is the resolution of the identity 
of the operator HE 00 (Q 00 ) acting in E oo (Q O0 )L2(M. 2 ) . 

7.4. Proof of Absolute Continuity. Now we show that the branch of spectrum 
(semi-axis) corresponding to is absolutely continuous. 

Theorem 7.12. For any F £ (M 2 ) and < e < 1, 

\{Eoo(Goo,x+e)F,F) - (E 0O (Q OOt x)F, F)\ < C F e. (7.12) 

Corollary 7.13. The spectrum of the operator HE 00 {Q 00 ) is absolutely continuous. 

Proof. By formula (j7.9j) . 

I (E 00 (G 00j x+ £ )F, F) - (E(Q 00 ^)F,F) \<Cf\Q oo,A+e \ foo,A| • 

Applying Lemmas 17.21 and 17.41 we immediately get (|7. 12|) . 

■ 

8. Geometrical Lemmas. 

8.1. Proof of Lemma 12.11 (Geometric Lemma, step I). Corollaries 18.51 and 18. 7| 

obtained in Sections 18.1.21 and 18.1.31 together give Lemma 12. II and its corollary. 

8.1.1. Set xi(W,si,5). The set xi(A,/3,si,(5) is given by formulas (gnity-i^Tfl). It 
remains to define radii r m of the discs O m . First, we need more notations. Angle tp m 
is defined as a polar coordinate of p m : 

Pm ■= Pm(0) = p m (cosip m ,sm(p m ), m = (mi,m 2 ) G Z 2 . 

In addition, let H m)qi ± be the distance in C between \k(<p) + p m \ 2 and \k{ip) + p m+q \ 2 
when ip = (fmi i- e -) when ip solves \k{<p) + p m \ 2 = k 2 . Obviously, 

n m , g ,± := \ (k(ip m ) +Pm,Pg)*+p 2 q \- (8.1) 
We will be interested in values of H m ,q,± for q : < p q < k Sl . Let 

Kn, q ,± ■= \(k<Pm) +Pm,P q }*\ + k 2si . (8.2) 
Obviously, U' m ^ ± = n mj9)± + O (k 2s i), U' m ^ ± > k 2s \ Next, 

Um := o< P ™ln , ± U ' m ^ = 0<P £& , ± 1 + ^' I + ^ ■ (8 - 3) 

The value II m characterizes the smallest distance between \k((p) + p m \ 2 and \k(<p) + 
Pm+q\ 2 when \k(ip) + p m \ 2 = k 2 , < p q < k Sl . It coincides with this distance up to 
the value 0(k 2si ). Further we use that II m > k 2si . 



SPECTRAL PROPERTIES OF A LIMIT-PERIODIC SCHRODINGER OPERATOR 73 

Definition 8.1. The radius r m of the two open disks O m in C centered at ip m is 
defined by 



k 



-3s\ 



p m Tl m 



if < p m < 4k Sl , 
if p rn > 4k Sl and 



Prn 

TP 



P7„ 
lk- 



> 



(8.4) 



Pr, 



if 



1 



l>7„ 
IP 



< 



Remark 8.2. If p m > 4k, then the imaginary part of (p m is too large to be in <Pq 
and so we do not need to remove corresponding disks from <Po to construct #i- 

Lemma 8.3. For every m with < p m < 4k, r m = o(l) as k — > oo. 

Proof. Let us consider the following four cases: 

If < p m < 4k Sl , then r m < k rzr^ = k~ 2si ~ s = o(l), since the periods a\, 02 of 
are of order k Sl . 



If 4k°i < Pm <k, then r 2 m < fc 4 fc Ji 3 / 4 < 2k^~ 2 ~^ = o(l), since U m > k 2s K 

< k 2/3 ~ 2 ~ 2si = o(l). 



If k < p m < 4k and 



Pm 

Lk- 



> 



k 2 P- 2 

n m 



, then r 2 ^ < 



If 



Pm 

TP 



< 



then of course p n 



k 2 m 

> k and r m < 



fc2,9-2 



m n m 
JkP 
k 1+a i 



0(1). 



8.1.2. Proof of the statements (i), (ii) in Lemma \2.1\ 

Lemma 8.4. Let 4s\ < 2/3 < 1 - 15si - 85, ip is in @i or its (k 2 P~ 3 ~ sl ~ 25 )- 
neighborhood, m,q G Z 2 , m,q,m + q ^ 0, < p q < k Sl . Then 



2\\k((p) +p m \l ~ k 2 \\\k{if)+p 



|2 



fc 2 | > A^, 



(8.5) 

\\Kv)+P q \l-k 2 \ >k 1 -^- 5 . (8.6) 

Corollary 8.5 (Statements (i), (ii) in Lemma l2.1| Corollary [22]). Ift G Xi(A, j3, s\, 5), 
then there is a unique j G Tl? such that Pj(t) = k and (|2,4p , (|2.5p /10W. For any t in 
the (/j2/3-2-si-25j .neighborhood of the non-resonance set in C 2 , there exists a unique 
j G Z 2 suc/i that (USD and (|2Ti|> . (|23]1 . fllgD (12401) ftoZe?. 



Proof of Corollary \8. 5[ If i G xi(A, /?, si, 5), then, by the definition of Xl> there is 
<p G 0i such that t = K\k(<p). Hence, there exists a j such that Pj(t) = k. Let us 
show that j is unique. Suppose there is f = j + m, m ^ such that pj +m {t) = k. 
Then pj +m {t) = k{ip) + p m and \k{ip) + p m \ = k. The last relation contradicts to 
([83]) . Therefore, j' = j. 

Substituting fc(c^) = p^-(i) into (|8.5p . (|8.6p and using the notation j + m = i, 
we obtain (|2.4p and (|2.5p . Similar arguments work when t is in the (k 2 P~ 2 ~ Sl ~ 2 ^)- 
neighborhood of xi(A, (3, si, 5) in C 2 . 

The inequalities (|2.9p . (|2.10p are obvious, when pi > 4k. We assume now pi < 
4k. Let us prove ([2J|. When i = j, §ZM) follows directly from §ZM and 
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Suppose i / j. Clearly, \p 2 (t) - p 2 i+q (t)\ = 0(k 1+Sl ). This together with {23J) yields: 
> cfe 2 ^~ 1_si . Considering ([%5|> . we obtain 



2|j>?(t)-z| > \pf{t) - p 2 At)\. 



(8.7) 



The estimate (|2.9jl easily follows. Let us prove (|2.10p . In the case of i = j, using 
(|2.5p and the definition of C\, we get 



\v] +q {t) -z\> \ v ){t) - v ) +q (t)\ - k 2 ^ 1 -^- 6 > h 



l-3si-<5 



Using this together with the definition of C\, we get (I2.10P in the case % = j. In the 
case of i ^ j (|2.4p and (|8.7p yield (|2.10p . It is easy to see that all the estimates are 
stable under a perturbation of t of order ^ 2 ^- 2 - s i- 2<5 . Therefore, the estimate (|2.10p 
can be extended to the complex (fc 2/3_2_,Sl_2<5 )-neighborhood of Xl(&; P,si,5). | 



Proof of Lemma \8.4\ If Pm > 4 A: then (I8.5P is obvious, since both factors are greater 
than k 2 . Assume < p m < 4k. Noting that 



\H<P) +Pm\l 



k' 



2kp m cos((p - (p m ) +p 



and recalling that tp^ are the solutions of \k(ip) + p r , 



k 2 , we see: 



cos{ip m - ip m ) = -—, I sin(p r 



<Pr. 



Pin 

Ak 2 



(8.9) 



UO) 



Now let tp be on the boundary of 0+ or O m . Expanding (|8.9p around tp^, we get: 

\k((p)+p m \l-k 2 = 2kp m sm(tp m -<p m )r m (l + 0(r 2 m )) -kp m cos(^-99 m )r 2 a (l + O(r^)) . 

(8.11) 

Next, we prove 

\\k(p) + p m \l ~ k 2 \ > k 2 ^ /U m . (8.12) 

in three cases as in Definition 18. lf i). In special, we get (I8.6P in the first case. 

Case (i) : r m = k — when < p m < 4k Sl . The modulus of the first term 



Pm 



111 



3D is 2/c 1 ~ 3si - ,5 a/i 



+ o(l)) > |A; 1 ~ 3si ~ <5 and that of the second term 



is 0(fc 6si 2<5 ), which is much smaller than the first term. Thus, we obtain 
Considering that /c 1-3si ~ 5 ^ - A *~ 



>f^,we obtain ([8J2D . 



Case (ii) : r n 



k 2 ?- 





V 2 




1 



when 4k Sl < p m < 4A; and 



Pm 

1A-- 



Substituting r m into (|8.1ip . we get that the modulus of the first term is 2 ^- (1 + o(l)) 



> 



and that of the second term is 



k 2 ? 



k^~ 







n 2 


^ Pm 


±1 m 


4fc 2 



(1 + o(l)). Using the condition 



/<7„ 

TP 



> 



— , one can easily see that the former is at least twice greater than the latter. 
Thus, we get (^~12l) . 
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Case (iii) : r„ 
second term is 



when 



1 



TP 



< 



-. This time the modulus of the 

16k 2 ~ 



1 + o(l)) and that of the first is smaller than -^-(1 + o(l)). 
Therefore we again have (|8.12p . 

Now we prove f)8.5[) . If < p m ,p m +q < 4k Sl , then (|8.5p easily follows from (|8.6p 
proven in the Case (i) above for such m,m + q. Next, if both factors in the left 
hand side of (|8.5p is greater than k 13 , then (|8.5p is obvious. Therefore, without loss 
of generality, we can assume 



\\k&) + Pm\t ~ k 2 \ <kP, Pm >4k S \ 

Suppose we have proved that 

\k(ip) +p m + q \l ~ k 2 



(8.13) 



U4) 



for every m : II m > k 13 and ip satisfying ((HUHJ. Then, notice that II m > A;' 3 follows 
from (l8TT2j) . Considering (j87T2j) and dgTHjl together, we obtain (j8"3j) . 
It remains to prove (|8.14|) . First, we check that 



\\H(p)+p m \l-\k(<p)+p m+q \l\ > -U m + 0(k 2si 



1.15) 



for every m : Tl m > k$ and ip satisfying (|8.13p . Indeed, taking into account (|8.1ip 
and considering as in the proof of (I8.12p . we conclude that the set of ip satisfying 
(|8.13p is inside the circle of the radius r' m around <£>^, where 



2 f~> 

--, when 1 - > k 13 



wk^ 2 -\ 



otherwise. 

It is easy to show that for any ip inside those circles: 



\k(ip) +Pm\l ~ \k{ip) +p m + 



2 

q I * 



\k{<P>ti) +Pm\t ~ |fc(<Pm) +Pm+ 



2 

q I * 



< 



(8.16) 



3.17) 



2kp q r' m < ~k? + o(l) < hl m + o(l). 
Considering the definition ()8.3p of II m , we easily get: 

{fat) +Pm\l ~ \k(ipt) +Pm+ q \l > m m + 0(k 2 ^ 



(8.18) 



Combining (fSTTTT) and (J87T8]), we arrive at (f87T5|) . Now inequalities <KW\ and (pTL5j) 
yield (l8TT4"l) . Thus, <^E) is proven. 

Note that all estimates are stable with respect the perturbation of ip of order of 
fc 2/3-3-si-2«5_ Therefore, (pT5l) . §8M hold not only in but also in its (Z^- 3 " 51 " 25 ) -neighborhood. 



The only conditions required for (|8.14[1 are n m > fc 8 and (|H7T3|) . We will use (f8T4)l under this 
condition in Lemma 18.91 



76 



Y. KARPESHINA AND Y.-R. LEE 



8.1.3. Proof of the statement (Hi) in Lemma [Ql 
Lemma 8.6. 

£ r m = 0(k- 5 ) 



m££" 

0<p m <4fc 3 l 



£ 



4fc s i<p m <4fc 



70^ 



3-19) 
3.20) 



where 



7o 



1 - 2/3 - 15si - 76, if p > 1/6; 
5/6 - /3 - 15si - 7<5 if p< 1/6, 



the sum ^ n <|8.20[) includes only such m that O m <t Uo<p 9 <4fc s i,±C^- 

Corollary 8.7. 7/2/3 < 1 - 15si - 8(5, f/ien f/ie totaZ /en#f/i o/ £>W n [0, 2vr) does noi 
exceed 0(k~ s ). Estimate (|2.7p holds. 

Proof. We split m G Z 2 with < p m < 4fc into five sets: 

Ji := {m G Z 2 : < p m < ik Sl } , 

jm G Z 2 : 4A; S1 < p m < fc l-4 sl -4«5 j ^ 

72 ■- >fc 1 - 4si - 45 , 2/fc-p m > k 1 - 8 ^ 46 }, 



h 
h 

h 



•fm G : p r 
m G Z 2 :2k-p m < fc 1-8si_4<5 , 



Pm 

Ak 2 



> 



4LL 



'5 " 



m G Z 



Pm 

Ak 2 



< 



k 2 ?~ 



4il 



and define Ej := X^me/ rm > J = ^> ••' ^ an< ^ an u PP er bound for each. 

<i) By ea, 



rnGh 



-3si—5 



Pn 



.21) 



Considering that the size of the lattice formed by p m is ck Sl , we obtain: 



5 Pm 



2.si 



rneh 



2vr /-fc s l 



JO 



(8.22) 



Substituting the last estimate into (|8.2ip . we obtain f|8. 19[) . 
(ii) Now we estimate £2. Let us check that 

n m > ^k l ~ 3si " 38 , when m G h- 

Indeed, suppose H m < 1 y/i; 1 ~ 3si ~ 3 ' 5 . Then for some q G Z 2 with < p q < k Sl : 
\(k{<£m)+Prn-,P q )*\+k 2si < ifc 1_3si_35 , where ± means + or -. Hence, |(«(<Pm)>jPg)*| < 
k 1 ' 331 ' 36 , since \{p m ,p q )*\ < k 1 ' 481 ' 46 ■ k Sl = k l - 3si ~ 4& when m G h- Therefore, 

\\kv m )+p q 



Z |2 



2(fc(</>± + p 2 \ < 2k 1 - 331 - 35 + /c 2si < SA; 1 - 351 ^, which 



,1— 3s\— 35 
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means ip^ G Of, see (|8.6p in Lemma [8,41 Moreover, </?^j is relatively close to the 
centrum of Of, its distance to the centrum of Of is less than r„k~ 2S . Using (JH 



one can easily check that r m « r q k 25 when m & I-i- Thus O^ is completely con- 
tained in Of . We do not consider such an m in the sum (18.20p . Thus, we can use 

n m > ik 1 - 3 * 1 - 35 . By (jsaj 

fc 2/?-l 

/ — - /,! 



me/5 ri TT 



1 Pm 

- 1 IF 



M— 3si— 35 / - ^ 
me/2 



We estimate the sum on the r.h.s. by the integral k 1 J J 4fcsi --rdrdip. Com- 
puting the integral, we obtain : 

(iii) Let m E I3. First we note that 



3.23) 



since 2A; — p m > A; 1-851-45 . For the moment we assume to have 
Then, 



3.24) 



.25) 



me/3 



E 



me/3 -n TT 



< 



1 Pm 

- 1 IF 



kl-4si-45 . i.— 4si— S 

me/3 



si-25 X/ TT 

_ -LJ-m. 



where we used the definition of ^3 and (|8.24|) . Considering also (|8.25p , we obtain: 

S3 < c^- 1+15si+7<5 . (8.26) 

Now we need to show (|8.25p . Obviously, ip^ are real, when m G Is- Therefore, all 
(•, ■)* and I • | 2 in the following formulas are real too. We note that |£?(92 m )| 2 = k 2 , 
\Kvti)+Prn\l = k 2 and \{k{tp±)+p m ,i> qm )*\ +k 2si = U m for some q m , < p Qm < k Sl 
and yj+ or /and <^~. Assume for definiteness that the last equality holds for 99+. 
Denoting fc(y TO ) + p m by k m gives us: 



I k m p n 



I ^m I * ^ 



I {^miPq m )* \ ~\~ k 1 — II 



(8.27) 
(8.28) 



Obviously, k 2si < U m < k 1+Sl . Define 0,- as the set of p m in I3 satisfying the 
inequality : k 2s ^ 5 < U m < fc2si+(j+i)5 , j = 0, • • • , J, J = L^J • Here |rj = 
max{z £Z:z< r}. Obviously, Uj C Uo<p 9 <fc s i where Qj >g consists of p m in I3 
satisfying (|8,27p and 



k 2s ^ 5 < \(k m ,p q )*\ + k 2s > < fc2»i+0-+D*. 



(8.29) 
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Let p m £ &j,q- By (I8.27D . ()8,29p . p m belongs to a circle of radius k centered at a 
point k m , where k m belongs to a circle of the radius k centered at the origin and 
satisfies the inequality k 2si+ i s < \(k m ,Pq)*\ + k 2si < / c 2s i+(i+ 1 )< 5 ) i. e . f~ m belongs to 
one of two arcs. Considering that p q > k~ Sl , we easily see that all such points p m 
belong to a couple of rings of the same radii k — k 3si+ ^ +1 ^ s , k + / c 3s i+0'+ 1 )< 5 centered 
at two points k q,ziz satisfying the conditions \k q,ziz \ = k, (k q)± ,p q ) = 0. Estimating the 
number of points p m at this region by its area, we obtain, that the number of points 
in Q,g does not exceed cfc 1+5si+ ^ +1 ^. Therefore, the number of points in f2j does 
not exceed cA; 1+9si+ ^ +1 ^. Thus, we obtain the estimate 

y -L<y y -L< y ck i+9s 1+ u+i)s . < ck w Sl+ s 
^ n m ~ ^ ^ n m ~ ^ k 2s ^ s ~ 

which finishes this case. 

(iv) Let m £ I4. For the moment, we assume the following two estimates, El 

n m > -k 1 - 9 ^ 45 , (8.30) 

v 1 f ck^ 2 +^+ s , if p > 1/6; 

- m % 4 ^2k^\ I ck-W^+ s , if P < 1/6. (8 - 31) 



Then, 



"16/4 n TT 



1 Pm 

IF 



/ ^-1+2/3+12^+5^ if ^ > 1/6; 
" \ cA; -5/6+/3+12 Sl+ 55 ! tf ^ < 1/6 

Now we show (|8.30p . We consider two cases: \2k — p m 

_£,l-8si-45_ s ^ ar ^ w ith the former. Suppose, (|8,30p is not true, i.e., H m < 
i^,i-9si-45 fphen, there is a q such that < p q < k si and 

|(^)+p m ,^)*|<^ 1 " 9S1 " 45 

for or/and Note that 93^ has a non-zero imaginary part when p m > 2k. 

Denoting k((fm) + p m by fc^, fc^ G C 2 , we get (|8.27p for both k^. It follows from 
(|8.27p that k^ = k((f_ m ), we mean here that the pairs of vectors are the same. 
Therefore, there is at least one (p^ m , such that \{k(ip^ m ),p q )*\ < \k 1 ~^ Sl ~ A5 . We 
denote it by <#L m . The last inequality means that tpt m is in O q , O q = C + U O q , 
and even relatively close to the center of a disc O^: its distance to the centrum is 
less than r q k~ 2si , the radius r q of being given by r q = k~^ Sl ~ 5 p q ~ l , see (|8.4p . 
The distance between centers of 0+ and 0~ is tt + 0(p q k^ 1 ), It is easy to see also 



12 To prove (|8.30[) we use only the following conditions: 2k — p m < —fc 1 8si 4,5 or \2k — p m \ < 
j c i-aa 1 -4S an£ j at j eagt Qne Q £ ^± - g j n ^/ _ ^ Q ^ Uo<p ? <4fc=i±C'* or its 2fc~ 4si ~ 2,5 -neighborhood. 

Under these conditions, we use estimate (|8.30p also in Lemmas 18.91 18.101 
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that (p^ m = + 7r modulo 2tt, we mean here that the pairs of angles are the same. 
Considering the last two relations, we obtain that y>+ or/and is also in O q . Using 
\2k - Vm\ < k l ~ 8si ~ 45 , we get from (l83U|) : \cp+ - ip m \ < 2k~ 4si ~ 2S . It is much 
less than r g . Hence, both and 99" are in either 0+ or 0~, relatively close to a 
center. It is easy to show that r m = o(k^~ 1 ) = o(r q ). This means O rn C O q . We 
don't include include such r m in the sum (I8.20|) . Thus, (I8.30P is proven for the case 
\2k- Pm \ < k 1 - 831 ' 45 . Let 2k- Pm < -k^ 831 ' 46 . By (l8T0|) . 



\(k(<pt)+p m ,p g )*\ = \(k(<ptj,p q )*\ = fep,|cos(^ m -^)| > -k Pq k- 8s ^ s > -k 1 -^-* 6 . 



It is easy to see that from the definition of ip^ that + p m = k(<p_ m ). Hence, 

h Pq k- 8si ' 45 > - 

Therefore, U m > i/c 1 " 951 " 45 . Thus, we proved (IQjjl . 

To show (|8.3ip , we split I 4 into three subsets Z4 = -Z4 + ^4 + ^4" : 

4 = { m G / 4 : |2fc - p m | < AT 1 / 3 }, J 4 ' = {m E / 4 : A;" 1 / 3 < \2k - p m \ < k 5 }, 

Jf = {m<El 4 : \2k - p m \ > k S }. 

Correspondingly, oj = oj' + oj" + uj'" . Let us estimate oj' . We use a well known 
estimate (see e.g. [27]) for the number N(k) of a rectangular lattice points in the 
circle of radius k: N(k) = irv^ik 2 + 0(k 2 ^ 3 ), where v is the area of the elementary 
cell of the lattice and the implicit constant depend on the periods. Considering 
that our lattice has a size of order k~ Sl , we rescale the last estimate, so it becomes: 
N(k) = ■Kv~ d 1 k 2+2si + 0(k 2 / 3+2si / 3 ), v d = A-k 2 ^ 1 ^ 1 . This means that I' A contains 
less than ck 2 ^ 3+2si points. Note also that \2k — p m \ > \k 2/3 ~ 2 ~ Sl ~ & follows from 
|1 - |§H > ^ and n m < k 1+s \ Hence, 

J < cv^k 2 / 3 ^ ■ k -^i+s l/2+ s/2 = o(yr /3+5/3+3 Si +^ 

Next, considering that the region I'l contains no more than ct/^ 1 /c 1+<5+2si of points 
and \2k — p m \ > k~ 1 / 3 , we obtain 

oj" < cv d l k l+5+2s ^ < cv d l k 7 / &+2 ^ +& . 

Further, we estimate oj'" by an integral: 

f 2n f 1 3 

oj'" < ck 2si v d ~ 1 / • / =r drdip < cv^k^ 

Jo J\2k-r\<2k y/\2k — r\ 

Combining the last three estimates, we obtain (|8.3ip . 

(v) As in the previous case, we have U m > |fc 1_9si_4<5 . The number of points in 1$ 
admits the estimate: #(^5) < cv ( ^ 1 A;3 +2si . Therefore, 

_ v ^ 4k? 8k? v ^ 

5 ~ ~ k • £;(l-9si-4<5)/2 " 



so 
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. Adding (1Q31) . (|Q6]l . (IQ2D . and the last estimate, 



Hence, S 5 < c^ 1 /t /3 -t +7si+2 ' 5 
we obtain (|H720|) . | 

8.2. Proof of Geometric Lemma 12.161 The set 0^(b) is defined by formula 
(12331) . where 0±(6) are open disks around tpm{b), the points being the two 

zeros in C of \k(tp) + p m {b)\ 2 = k 2 . Radius r m {b) of 0^{b) is defined in Section [8.2.1i 
Lemma 12.161 is the combination of Lemma 18.91 proven in Section 18.2.21 and Lemma 
18. 101 proven in Section [8.2.31 

8.2.1. Definition of the set 0^\b) . Let us recall that r m are the radii of the discs 
of the first non-resonant set 0^ l \ They are defined by formula (|8.4p . The radii 
f m (b) of the discs in 0^(b) are, roughly speaking, defined as follows: r m (b) = 
£.-2+2/31+2251+15^ when m = 0; r m (b) w r m , when p m (b) > 3k Sl ,\2k - p m (b)\ > 
fc i-8si-4<5. rm fy << when \2k-p m (b)\ < k 1 - 831 -^ or p m (b) < 3/c Sl . We have to 
reduce the r m (b) for smaller p m (b) and for p m (b) close to 2k in order to ensure that 
each component of 0^'(b) is sufficiently small (by (|8.4|) . r m becomes bigger, when 
p m is small or tends to 2k). 

Definition 8.8. The radius r m (b) of the two open disks 0^(b) in C centered at 
<p m (b) are defined by 

r fc-2+2/3i+22si+155 if m = 0, 



r m (b) := < 



£,-2+2ft+3si+25 

Pm{b) 

fc -l+2/3i 



Pm (b)U m (b) 



£,-4+2ft + 10si+45 



P 2 m (b) 
4fc 2 



if m ^ 0,p m (b) < A; 1 ~ 8s i- 4 ' 5 , 

if Pm(&) > k 1 - 8 ^- 46 , \2k-p m (b)\ > A; 1 - 881 " 45 , 



P 2 m(b) 

Ak 2 



4^,-2+/3i+5si+2<5 



if |2/c -p m {b)\ < k 



l-8si-45 



1 _ P 2 m (b) 



Ik? 



> l^.-4+2/3i+10si+45 
— 4"' ) 



if 



1 _ P 2 m {b) 



"IF" 



< l£,-4+2/3i + 10si+45 ; 



3.33) 



where n m (6) is defined by (|8.3p . p m (b) being substituted instead of p m . 

Considering as in Lemma 18.31 we easily show r m (b) = o(l) as k — > oo. 
8.2.2. Lemma\£M 

Lemma 8.9. Let lOOsi < A < 1/12 - 28si - 145. Suppose if G S' \ 0^\b). Then 
for every m E Z 2 such that 

\\k&) + p m (b)\l - k 2 \ < (8.34) 
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the following inequalities hold: 



min \\k(<p)+p m+a (b)\l-k 2 \ > 



, ■ 'in— <i 1 " ' u - '< | '■ (8.35) 

0<p q <k s l 

| l^(c^) ^^T^(^)|2 — A: 2 | | ^^^^ _ fc 2 1 1 ^^^^^ (g*)|2 _ fc 2| > k 2 Pl) (8 . 36) 
when < p qi (0),p q2 (0) < k Sl . This property is preserved in the / c - 4 + 2 / 3 i^ 2s i^ 5 neigh- 
borhood o/#' \OW(6). 

Proof. First, we consider the case m = 0. It is easy to see that |fc(y) + b\ 2 — k 2 \ > 
kboro, when ip Oq (6). Using the estimate 6o > ^ _1_16,S:L_125 and the formula for 
ro, we get 

||jfc + 6 |2_jfe2| > fc -2+2 / 3 1 +6 S1 +3<5_ (8>37) 

Next, 



|l^)+p gi WI*-^| = ||fcM+p 9l 



fcM,P«(0)) + 0(F 1 ), z = 1,2. 



Considering that ^ S #'o, we obtain 



l-3si-<5 



Hence, 



\\k( V )+p qi (b)\l-k 2 \>-k 



lSsi-5 



1,2. 



(8.38) 



(8.39) 



Thus, (fOST) holds for m = 0. Multiplying (15371) and (jOl for i = 1,2, we arrive at 
dHJSD for m = 0. 

Let p m (b) < Sk 1 - 831 - 46 , m ^ 0. It easily follows from the definition of 0±{b), 
+Pm(b)\l - k 2 > kp m (b)r m > fc -i+2/Ji+3»i+M 



Next, 

\k(<f) +p m +q t (b)\l ~ \H<f)+Pr. 



k{if),p q M) +0(k 



.1-781-45% 



(8.40) 



1,2. 



Considering that ip £ <P' Q , we again obtain (I8.38j) . Hence, 

\k(ip) +p m+qx {b)\l ~ \H<P) +Pm(b)\ 2 

Taking into account (I8.34p . we obtain 

\k&) + p m+q M\l ~ k 2 



i 



>-k 



3.41) 



Thus, (JO?l) holds for the case p m (b) < 3k^ Ssi - /iS . Multiplying (1CTD and (IH^Tj) for 
i=l,2, we obtain ([OS]) . 

If p m (b) > /c 1_8si ~ 4<5 , |2A; — p m (&)| > A; 1-8 * 1-45 , we use the same considerations as 
in Lemma 18.41 to show that 



|&(¥>) +Pm(ft)|* - A; 2 |fc(ip)+p m + g 



The estimates (I8.35p . (I8.36P easily follow. 



> k 2fSl , for i = 1,2. 
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Let \2k — p m (b)\ < k 1 ~ 8si ~ 4S . Using (|8.1ip and considering as in Lemma 18.41 we 
obtain that 

\\k(<p) + p m (b)\l - k 2 \ > k - 2+2 ^ +Ws ' +4S (8.42) 

when p OW(6). Suppose at least one of Pm(b) is in the (fc _4si_25 )-neighborhood 
of S' . Then, IL m (b) > \k l ~ 9si ~ 4S , see formula flOOj) and the footnote there. Using 
the inequality (|8.14p . see the footnote there, we obtain: 

\\k(cp)+p m+gi (b)\ 2 -k 2 \ > in m (6), i = l,2. (8.43) 

Estimates (|8.35p . (|8.36p easily follow. To finish the proof it remains to show that at 
least one p^(b) is in the k~ 4si ~ 2S neighborhood of &' . If both f^(b) are not in the 
^,-4si-2(5 ne igh Dor h ooc i f then each is inside a disc O^, < p q < k Sl and further 
then 2k~ 4si ~ 2S from its boundary. Therefore, the distance of p to p^ is greater than 
k- 4si ~ 25 . However, (18341 and (KIM , yield \p-p^\ < lOfc^/ 2-1 . This contradiction 
proves that this case is not possible. Since all estimates are stable with respect to 
the perturbation of cp of order ^- 4 + 2 / 3 i^ 2s i^ <5 5 the statement of the lemma holds in 
the A;-4+2/3i-2 Sl -5 neighborhood of S' \ 0^\b) too. ■ 

8.2.3. Lemma\8lM- 

Lemma 8.10. Suppose lOOsi < Pi < 1/12 - 28si - 14<5 and { c\b) is a connected 
component of 0^(b). Then, the size of 0^\b) does not exceed ck~^ , 7 = 11/6 — 
Pi — 12si — 45. A component 0^\b) contains no more than c\k 2 ^ +2si discs. The 
total size of 0^\b) does not exceed 2ircik~ 5 ^ + ^ 1+12si+4S . The total number of discs 
in 0^\b) is less than ck 2+2si . 

Proof. Let A* be a rectangle in C: A* = {cp : \3lp — p*\ < k~ 1 , \$scp\ < k~ s } for some 
<p* G [0,27r). Clearly, <Pq is the union of such rectangles. Let 

/* = {m G Z 2 , p m (b) < 4k, p±(b) G A,}, 

where p^(p) € A* means that either ^(6) G A* or p^b) G A*. We will show that 
the number of points in i* does not exceed 6co& 1+2si and 

£* := Yl rm $) < ck-~< (8.44) 

This means that the total size of (6)n A* is less than k^" 1 '. Considering that 7 > 1, 
we easily obtain that the size of each (b) is less then cfc -7 . Since ^0 consists of 
2nk rectangles A*, the total size of 0^\b) does not exceed fc-5/6+/3 1 +i2 Sl +4<S_ Let us 
prove (|8.44p . In accordance with the Definition 18.81 and by analogy with the proof of 
Lemma 18.61 we split i* into six sets: 

/,„ : = { m = 0}, I*i := {m G J* \ {0} : p m (b) < k 1 - 8 ^' 45 }, 

J, 2 : = { m G 7* : p m (b) > k 1 " 8 ^ 45 , 2k - p m (b) > k 1 - 8 ^ 45 }, 

1*3 := { m G I* ■ Pm (b) -2k> k 1 - 8 ^- 46 , }, 
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'*4 



m £ h : \2k - p m (b)\ < k 1 - 8 ^- 45 , 



1 Pm(b) 

4k 2 



> l /t -4+2/3i + 10s 1 +45 



1*5 ■= \ m G I* : 



1 



Pm(&) 



4/c 2 



1 



< - fc -4+2/3i + 10si+4<5 



We define E*j := ^ me/(j r m (b), j = 0, ..,5, and get an upper bound for each sum. 
Formula 

E„ = k - 2 + 2 ^+ 22 ^s (8.45) 
immediately follows from Definition 18.81 Obviously, 

E*i = k -2+2h+3s 1+ 28 J2 — ^. (8.46) 

We estimate the sum in the right-hand side by an integral: 

V — U- < k 2si + k 2si [ -rdrd^ < k 3s \ 

where 1^ is a (ck~ Sl ) neighborhood of Substituting the last estimate into ()8.46|) . 
we obtain: 

E*i < ck - 2 + 2 ^+^+ 2S . (8.47) 

since 2k — p m > 

(8.48) 



Let us estimate E* 2 - First we note that \J I - > jk- 4si - 2S 

^i-8si-4<5_ p or m0 ment we assume that 

v — U- < ck 7s i+ 35 . 

Then, 

k -l+2 h C k- l+2 ^ ^ 

/ ^ / — t.1— 8sj— 45 . I.— 4si— 25 / / 



J *2 



ie/»2 p m II m (&Wl 



1 



l>T„ 
Ik' 1 



n m (6) 



Thus, 



E* 2 < dfe- 2+a *+ 19 *+ M . 
Let us prove dSaSD. Definition of p±(6) yields: 



(8.49) 

fc. (8.50) 

k, \k(tp*) + p m (6)| = 



Considering that \<Pm(fi) ~ f*\ < we obtain: 

A; + O(l), |0(1)| < 1. This means that p m (b) belongs to the ring 1Z of radii k ± 1 
centered at (</?*)• We split 7£ into several components depending on the value of 
n m (6). Indeed, let 

fi*i = {m G Z 2 : G 71, £; 2si+ ^ < II m (&) < jfe 2 *^^ 1 )*}, 



j = 0, • • • , J, J = L^x 1 ] ■ Let us show that 



s m6 o n- 1 ^) < 4fe 7si 



7s!+28 



3.51) 



8-1 
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We estimate the number of points in 0*,-. Indeed, by the definition of LT m for every 
m G r2*j there is a p qm , such that < p Qm < k Sl and 

+ k 2s > < k 2s ^ +1 ^ s . 



< 



k (ymQ>)) +Pm(b),p qrl 



> 



Considering that \<f m (b) — <p*\ < k \ we obtain that 

< k{^*)+Pm{b),P q > < 2k 2si+{j+l)& for a q : < p q < k s \ and p m G TZ. 

The number of such points in TZ, obviously, does not exceed 4/c 5si+ ^ +1 ^ for a fixed q 
and 4^ 9s i+0+ 1 ) 5 for all q. Using now the estimate U m > k 2si+ ^ , we arrive at (I8.5ip . 
The estimate (|8.48p easily follows. 

Let us estimate E* 3 . Using ([QUI) , we get r m < 2k- 3+2fSl+13si+(i5 . By (|OJ|) . 
Sf{^ + 7T = ip m (mod 27r) when m G 1*3. Therefore, for any m G 1*3, we have 
Wm — ~k — V 9 *! < & an( i 2A; < p m < Ak. Now it is obvious that the number of points 
in 1*3 does not exceed ck. The estimate 

x* 3 < C £r 2+2 ^ +13si + 65 



easily follows. Next, we estimate E*4. Suppose we have checked that 



E 



< a 



^5/3-/31+25! 



(8.52) 
(8.53) 



Then, 



11*4 — 



,e/*4 y |2fc -p m (6)| 

^,-4+2/3 1 + 10si+45 



1 



/'7» 



< cA;~5 +2/3l+10 - Sl+4<5 



Therefore, 

E* 4 < c £r n / 6+ ^ +12si+45 . (8.54) 

To show (1533"j) . we note that \2k - p m \ > lk~ 3+2l3l+l0si+AS follows from 1 - |p > 
iA;~ 4+2 ^ 1+10si+4<5 . We split 1*4 into three regions: 

w*i := {m G 1,4, ^ < |2& - p m | < A: 1 - 881 " 45 }, 

w*2 := {m G 7*4, < \2k — p m \ < k S }, 

u,* 3 := {m G 1*4, h- 3+2 ^ +10s ^ 45 < \2k-p m \ < k- 1 }. 

The corresponding sums we denote as er,-, j = 1, 2, 3. It is easy to estimate <ti by an 
integral: 

,,£,1-86! -4i 



E 



< ck 



k 6 



2n 



1 

Vt 



d<pdt < ck^ 



-4s!-25 



Now we estimate 02- It is easy to see that the number of points in does not 
exceed cofc 1+<5+2si . Using the estimate 2k — p m > A; -1 , we obtain: 

a 2 < ck 3 ' 2+2s - +5 . 
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Let us estimate 03 . Definition of f^ib) yields k ^V'm(^)) = ^ 2 > & ^m(^)) + Pmib) 
k 2 . By (|QIj> . |9<p±| < AT 1 when m G w*3. Therefore, \ip^ — (p*\ < 2k 1 . Hence, 
we obtain: |fc(y*)| = A;, + p m (&)| = k + 0(l), |0(1)| < 1. This means that 

Pmip) belongs to the ring 1Z of radii k ± 1 centered at &(<£>*). We also introduce the 
ring 1Z* = {k : \k\ — 2k < k^ 1 }. Obviously w*3 C 7*!.*. Thus, u>*3 C 7£ n 7£*. It is 

not difficult to show that the ring 7Z cuts from 7Z* an "arc" of the length 0(k 1 ^ 2 ). 
Using well-known results from the theory of lattices, see [27], we obtain that IZnTZ* 

contains no more than ck l ^ +Sl points p m (b) for any b. Considering the estimate 
\2k — p m \ > 1/^-3+2/31+1051+4(5^ we arr i ve at fa e inequality: 

a 3 < ck^-P^- 25 . 

Adding the estimates for a\, 02 and (73 and considering that fti < 1/6 — 10si — 115, 
we obtain (|8.53|) . 

Let us estimate £*5. The number of points in 7*5 can be estimated the same way 
as the number of points in w*3, i.e., it is less than cfc 1//6+Sl . Hence, 

£* 5 < dfe-ii/e+ft+e-i+M. (8.55) 

Adding the estimates ([5~15|> . ([8171) . ([532]) . (18^91) . ([531]) and (1835]) and considering 
that /9i < 1/6 - 10si - 115, we obtain ([513]) . 

Let us show that the number of points in I* does not exceed co/c 1+2si . First, 
we count all real (p^(b). Obviously, \tp* — (p^(b)\ < k^ 1 . Using (|8.50p . we obtain 
+ Pm(b)\ = k + 0(1), I O < 1. The number of points p m (b) in this ring 
does not exceed cok 1+2si , s\ appearing since the size of our lattice is of order k~ Sl . If 
9(^(6) / 0, then, by (I8.10p . cp m = SR^(6) + tt modulo 2tt. Hence, \cp m — tt — ip*\ < 
k^ 1 . Considering that p m {b) < 4/c, we obtain that the number of such points does not 
exceed 4cok 1+2si . Hence, the total number of points in I* does not exceed 6cofc 1+2si . 

It remains to show that a component O^p (b) contains no more than c\k 2 ^ +Sl 
discs. Indeed, the length of Oc (b) is less then c/c~ 7 . This means that all points 
<fm(b) belonging to 0^\b), are, in fact, in a square of the size c/c~ 7 centered at 
a point y?** € A*. Let us consider all real ip^i") m t ne component. Obviously 
|%**-^±(6)| < ck-y. Using (1530]) . we obtain +p m (b)\ = k + Oik 1 ^). By 

[27] , the number of points p m (b) in this ring does not exceed c\k 2 / 3+Sl , s± appearing 
since the size of our lattice is of order k~ Sl . If Qc/?^(6) / 0, then, by (|8.10p . (p m = 
?RiPm(b) + tt modulo 2tt. and cosh Sy>m (b) = p m /2k. Hence, \tp m — it — 3J<^**| < c/c~ 7 
and \p m — < c/c 1-7 , 2** = k cosh 3^**. Obviously the number of points p m (b) 
is such a region is 0(1). Hence the total number of points (f^ip) in 0^\b) is less 
than c\k 2 l^ +Sl . Considering that the number of points p m (b) satisfying the inequality 
Pm(b) < 4/c is less than ck 2+2si , we obtain that the total number of discs in 0^'(b) 
is less than ck 2+2si . ■ 
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Appendix 1. Proof of Lemma l2.30L By Lemma l2,lCH Part 2. the function xj (ip) is 
holomorphic in Si and AW(j?i (<p)) = k 2 . Hence, the equation (|2,130p is equivalent to 
AW(yW (</?)) = AW(yW(»-6) + eo. We use perturbation formula ([238]): \y^(<p)\l + 
/lfy^Hv)) = ly^H 1 / 9 ) ~~ ^1* + fi(y^ 1 \ i P) ~ fy + e o- This equation can be rewritten as 

2<y< 1 V), 6), - |6| 2 + - /i^^M - S) - e = 0. (9.1) 

Using the notation 6 = 6o(cos sin^ft), dividing both sides of the equation (|9.ip by 
2bok, and considering that = xi(<£>) + & = (fc + + we obtain: 

cos(ip - if b ) - e gi(cp) + g 2 (<p) = 0, (9.2) 

where gi{<p) = (2bok)^ 1 and 

{hi(tp),b}* b 



92(<f) 



b k 



Let us estimate #2 (</?)• Using the inequality (12.70P for hi, and considering that 60 < 
k-l-l6si -12<^ we eas iiy obtain: 



(hi(<p),b), 



bnk 



< 



2/n 



-4+30si+20<5 



), 



^0 < 1 ju-2-16si-125 
2fe - 2 



therefore 



< sup|V/i|6o = O(k- 2+33si+228 b ), and 

(/i(/^))-/l(^%) -&))<?! (¥>)| = 0(£r 3+33si + 225 ). Thus, we have 
= 0(£r 2 - 16s i~ 125 ). Using e < 6 ^ 1_16si " 12<5 , we obtain e 5i(<^) < A;- 2 - 16 * 1 " 125 . 

Thus, 

32M-e giM = O(A;- 2 - 16si - 125 ). (9.3) 



7T 



Suppose that ^ ± — is in the 2 16si -neighborhood of Si. We draw two 
circles C± centered at ipi,± — with the radius |/c~ 2 ~ 16si ~ 11<5 . They are both inside the 

complex 2/c~ 2_16si_11<5 -neighborhood of <Pi, the perturbation series converging and 
the estimate ([93]) holds. For any tp on C±, \<p - {(p b ± w/2)\ = ^k~ 2 ~ 16ai ~ ns and, 
therefore, |cos(<^ - <p b )\ > ^fc" 2 " 1651 " 115 > \g 2 ((p) - e gi(ip)\ for any (p G C±. By 
Rouche's Theorem, there is only one solution of the equation (|9.2p inside each C±. 
If ipb + vr/2 is not in the (|/c _2 ~ 16si ~ 115 )— neighborhood of Si, then | cos((/9 — ifb)\ > 
Ygfc~ 2 ~ 16si ~ 11<5 in Si and, hence, equation (|9.2|) has no solution. Thus, there are at 
most two solutions in Si and \tpf - (ipb ± vr/2)| < |fc _2 ~ 16,Sl ~ 115 . 



Appendix 2. Proof of Lemma 12.311 Using the perturbation formula (|2.58p . we 
obtain: 



d_ 

dip 
d 



V^ 1 ) (y« fo,)) - V f A« (y« ( V ) - 6) , (?) 



SPECTRAL PROPERTIES OF A LIMIT-PERIODIC SCHRODINGER OPERATOR 87 

n2 



V|^V)|: - V\y^(<p) - b\l, (k + h 1 )fl + h' 1 u^ + 

(v/ x (y(%)) - V/i^fa) -b),(k + h l )jl+K^)^ , (9-4) 
where V = (cos 93, sin (p) and /2 = V = (— sin 99, cos 93). Note that 

V|y<%)|* - V|^%) - b\l = yW{y) - 2(tfW{<p) -b)=2b (9.5) 
Substituting (£5} into ([93]), we get ^A^ (y^M) = Zi + T 2 , 
Ti = 2(b,(k + h 1 )fl+h' 1 P^ , 

T 2 = ( V/i (y^ (*>)) - V/i (y« (p) - S) , {k + h x )ji + h! x V 



Considering that if is close to ± ir/2, we readily obtain: (b, z?)* = o(bo), (6,/!)* = 
±6 (1 + o(l)). Using also estimates (p770j) for /ii, we get Ti = ±2b k(l + o(l)). By 

(M), |v/i - V/i^fcp) -b) = O(b k- 2+36s ^ s ). Hence, T 2 = o{b k). 

Adding the estimates for Ti,T 2 , we get (|2.132p . 

Appendix 3. Proof of Corollary 13.61 Let To G % 2 . Taking into account the 
relation X^\tq + 2irp/N\a) = k 2 and the definition of C 2 , we see that \\f\ro + 
2-Kp/N\a) — z\ = ei/2. Using (|3. 13|) and the last equality, we easily obtain: \\n \i~o + 
2-Kp/Nia) -z\> ei/2 for A&^To + 2-Kp/N ia ) ^ A^ro + 2Trp/N ia ). Therefore, for 
any 2 G C 2 , 

\\(m i Hr )-zr 1 \\<2/e 1 , (9.6) 

i.e., ()3. 18f) is proved for to G X2- Now we consider r in the complex (ei/c -1-5 )— neighborhood 
of X2- By Hilbert relation, 

(tfW(r) - z)- 1 = (#«(r ) - z)" 1 + T x T 2 {hM(t) - z)~\ 
T\ = (H(t ) - z)- 1 (H q (t q ) + k 2 ), T 2 = (H (t ) + k 2 )- l (H Q ( TQ ) - H (t)). 



Suppose we have checked that HT1T2H < k 5 . Then, using (|9.6p . we easily arrive 
at (j3. 18|) . The estimate ||Ti|| < 4/c 2 /ei, easily follows from (|9.6j) . The estimate 
||T 2 || < 2eik- 2 - 5 easily follows from |r - r | < e^ 1 ^ . Thus, ||riT 2 [| < 8k- 5 and, 
hence, (13. 18H is proved. 
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